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ABSTRACT 


It  is  often  desirable  to  simplify  multivariable  feedback 
controllers  designed  using  modem  control  theory  so  that  they  will  be 
easier  and  more  practical  to  apply  to  the  control  of  industrial  pro¬ 
cesses.  Typical  simplifications  include  linearization,  assumption 
of  time- invariant  models,  and  reduction  of  the  order  of  the  govern¬ 
ing  equations.  This  thesis  is  concerned  primarily  with: 

1.  The  simplification  of  linear  time- invariant  state  space 
models  using  modal  analysis  techniques  to  eliminate  selected 
state  variables  followed  by  the  design  of  multivariable  con¬ 
trollers  based  on  these  reduced  models. 

2.  The  use  of  modal  analysis  techniques  to  eliminate  selected 
state  variables  from  a  "high  order"  control  law  designed 
using  the  original  process  model. 

3.  The  use  of  least  squares  techniques,  instead  of  the  modal 
analysis  mentioned  above,  to  fit  the  low  order  system  (model 
or  control  law)  to  data  generated  by  exciting  the  high  order 
system  with  a  sequence  of  random  numbers. 

4.  The  development  of  computer  programs  to  implement  the  above 
techniques  plus  examination  of  related  theoretical  points 
such  as  the  equivalence  of  different  reduction  techniques. 

All  of  the  techniques  developed  in  this  work  were  applied 
to  the  computer  controlled  pilot  plant  evaporator  at  the  University 
of  Alberta  that  was  used  in  previous  multivariable  control  studies. 
Simulated  and/or  experimental  data  are  presented  to  show  the  relative 


iv  - 


■ 


. 


. 


performance  of  different  process  models  and  controllers.  Specific 
contributions  and  conclusions  are: 

1.  Existing  modal  approaches  to  the  reduction  of  continuous¬ 
time  models  are  extended  for  use  with  discrete-time  models. 
Several  of  the  basic  modal  results  presented  in  the  litera¬ 
ture  are  shown  to  be  equivalent. 

2.  In  calculating  a  reduced  order  discrete-time  model  from  a  high 
order  continuous-time  model  using  the  modal  approach,  the 
relative  order  of  the  model  reduction  and  the  model  discret¬ 
ization  steps  is  unimportant,  since  the  same  reduced  order 
model  results.  However,  using  the  least  squares  approach, 

and  the  same  sequence  of  random  numbers,  different  reduced 
order  discrete-time  models  result  when  the  order  of  the 
reduction  step  and  the  discretization  step  is  interchanged. 

3.  Control  laws  designed  by  eliminating  selected  state  variables 
from  a  high  order  control  law  using  the  modal  analysis  were 
the  only  controllers  to  give  consistently  good  control  of  the 
evaporator. 

4.  When  the  least  squares  approach  with  random  data  is  used  to 
simplify  control  laws,  the  result  obtained  by  minimizing 
the  sum  of  the  squares  of  the  difference  between  the  orig¬ 
inal  control  vector  and  the  desired  reduced  order  control 
vector  is  the  same  as  that  obtained  by  minimizing  the  sum  of 
the  squares  of  the  difference  between  the  closed-loop  responses 
of  the  high  order  model,  when  controlled  by  the  high  order 
controller  and  by  the  desired  reduced  order  controller. 
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In  general,  it  was  concluded  that  a  multivariable  control 
law,  such  as  the  tenth  order  optimal  controller  designed  for  the  pilot 
plant  evaporator,  can  be  simplified  by  the  elimination  of  selected 
state  variables  and  that  the  resulting  control  schemes  can  be  easily 
and  successfully  implemented  on  actual  process  units. 
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CHAPTER  ONE 


INTRODUCTION 

1.1  GENERAL 

According  to  the  Encyclopedia  Britannica,  a  control  system 
is  any  ineans;  natural  or  artificial,  by  which  a  variable  quantity  or 
set  of  variable  quantities  is  caused  to  conform,  more  or  less,  to 
some  prescribed  norm  [29].  The  most  common  type  of  control  system 
presently  implemented  in  the  chemical  and  petroleum  industry  is  a 
single  variable  feedback  system  applied  through  the  use  of  electronic 
and  pneumatic  instrumentation.  The  theory  and  application  of  these 
single  loop  controllers  has  been  discussed  in  detail  in  many  text¬ 
books.  [5, 6, 11, 22 ] .  They  are  relatively  easy  to  apply  to  the  process 
and  can  provide  good  control  of  the  variable  of  interest.  However, 
difficulties  can  develop  when  process  variables  interact  or  when  the 
control  variable  for  one  section  of  the  plant  is  a  feed  variable  for 
another  section.  Some  modifications,  such  as  feedforward  control, 
help  to  overcome  this  problem.  A  more  recent  approach  to  eliminating 
these  interactions  is  the  use  of  a  multivariable  control  scheme. 

Optimal  multivariable  control  and  its  theoretical  potential 
has  been  discussed  at  length  [3,13,14,21].  However,  its  application 
to  the  control  of  a  physical  process  is  more  difficult  than  the 
application  of  single  variable  controllers.  This  is  because  each 
control  value  is  calculated  as  a  function  of  many  process  values 
(states).  This  practical  application  has  become  feasible  with  the 
development  of  the  control  computer  capable  of  performing  the  required 
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on- line  calculations.  Also,  in  the  application  of  optimal  multivariable 
control,  values  must  be  obtained  for  each  of  the  process  variables  that 
are  in  the  state  vector  of  the  process  model.  Usually,  these  cannot 
all  be  measured  so  that  estimated  values  must  be  obtained.  This,  in 
turn,  requires  more  on-line  computer  calculation  time  and  so  increases 
the  loading  on  the  computer.  An  alternate  approach  for  dealing  with 
unmeasurable  states  is  to  eliminate  them  from  the  control  law.  It  is 
this  approach  which  is  treated  in  this  thesis,  using  two  general  tech¬ 
niques.  In  the  first,  the  unmeasurable  states  are  eliminated  from 
the  process  model  and  the  resulting  model  is  used  to  calculate  the 
optimal  multivariable  control  law  (model  reduction)^  in  the  second, 
the  unmeasurable  states  are  eliminated  from  the  optimal  control  law 
itself  (control  law  reduction) . 

Previous  control  studies  at  the  University  of  Alberta  include 
the  following.  Andre  [2]  built  a  pilot  plant  evaporator  and  controlled 
it  with  conventional  single  loop  feedback  electronic  instruments. 

Wilson  [25]  controlled  the  pilot  plant  with  electronic  feedforward 
plus  feedback  controllers.  Jacobson  [12]  and  Fehr  [7]  performed 
similar  studies  using  direct  digital  control  on  the  newly  acquired 
IBM  1800  Data  Acquisition  and  Control  Computer  [8,28].  They  also 
studied  inferential  control  of  the  product  concentration;  controlling 
the  concentration  using  calculated  values  instead  of  measured  values. 
Newell  [17]  applied  optimal  multivariable  control  with  proportional 
feedback,  integral  feedback,  feedforward  and  setpoint  control  modes. 
Additional  advanced  control  techniques  have  been  applied  to  the 
evaporator  by  Nieman  [ 18]  (optimal  state  driving  control),  Oliver 
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[19]  (model  reference  adaptive  control),  Alevisakis  [1]  (time  delay 
compensation)  and  Hamilton  [10]  (Kalman  filters  and  Luenberger 
observers) . 

Control  studies  at  the  University  of  Alberta  using  a  pilot 
plant  distillation  column  [4,16,20,23]  and  identification  studies  on 
a  heat  exchanger  [15,24]  have  also  been  done. 

1.2  SCOPE  AND  STRUCTURE  OF  THE  THESIS 

This  thesis  treats  the  general  problem  of  controlling  a  pro¬ 
cess  which  has  unmeasurable  states  in  the  process  model.  The  basic 
model  for  this  study  is  a  tenth  order  lineat  model  of  a  pilot  plant 
evaporator  which  was  obtained  by  a  numerical  linearization  of  the 
nonlinear  model  developed  by  Newell  [17] .  The  general  form  of  the 
model  is  in  one  of  two  forms;  the  continuous- time  form  as 

x  =  Ax  +  Bu  +  Dd_  (1.1) 

JZ  =  C  x  (1.2) 

or  the  discrete-time  form  as 

x(j  +  1)  -  0  x(j)  +  A  u(j)  +  0  d(j)  (1.3) 

Z(j)=Cx(j)  (1.4) 

where  x  is  the  n-dimensional  state  vector 

u  is  the  m-dimensional  control  vector 
d  is  the  q-dimensional  disturbance  vector 
X.  is  the  p-dimens ional  output  vector 


=  >  =>  =’  E->  <k)  A,  §  are  constant  coefficient  matrices  of  approp¬ 
riate  dimensions 

j  is  a  counter  such  that  x(j)  =  x(t)  when  t  =  jT  and 
T  is  the  discrete-time  interval 

The  basic  control  law  used  is  a  linear  proportional  feedback, 
integral  feedback,  feedforward  plus  setpoint  controller  which  results 
from  calculations  using  dynamic  programming  techniques  to  minimize 
a  quadratic  performance  criterion.  Unless  otherwise  stated,  the  con¬ 
trol  laws  are  discrete-time  controllers  so  that  they  can  be  applied 
to  the  control  of  the  pilot  plant  evaporator  using  a  digital  computer. 

The  topics  discussed  in  this  thesis  can  be  outlined  as  part 
of  the  general  scheme  presented  in  Figure  1.1.  The  model  reduction 
steps  represented  by  path  one  for  a  continuous- time  model  and  by 
path  four  for  a  discrete-time  model  are  treated  in  Chapters  Three, 

Four,  and  Six.  Procedures  for  control  law  reduction  represented 
by  path  eight  are  developed  in  Chapters  Five  and  Six.  The  other 
topic  discussed  in  the  thesis  is  the  calculation  of  a  reduced  order 
control  law  as  the  optimal  controller  for  a  reduced  order  model  which 
had  been  obtained  by  model  reduction  techniques  (Part  of  Chapter  6). 

This  is  shown  in  Figure  1.1  by  path  four  followed  by  path  seven. 

The  thesis  is  structured  in  the  general  form  of  a  series 
of  publications  dealing  with  the  methods  outlined  in  Figure  1.1. 

Earlier  versions  of  several  chapters  have  been  published,  as  noted, 
in  the  discussion  which  follows. 
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FIGURE  1.1  PATHS  FOR  THE  CALCULATION  OF  LOW?:  ORDER  MODELS 
AND  LOW  ORDER  CONTROL  LAWS  FROM 
HIGH  ORDER  MODELS 
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In  Chapter  Two  a  GEneral  Multipurpose,  Simulation  and 
COntrol  PackagE  (GEMSCOPE)  is  presented.  For  ease  in  simulation 
and  design  calculations,  a  collection  of  computer  programs  was 
assembled  into  one  integrated  package.  This  package  allows  the  con¬ 
trol  system  designer  to  specify  a  model,  perform  various  calculations 
on  his  model  (model  reduction,  calculate  the  discrete-time  model, 
calculate  a  discrete-time  model  with  a  different  sampling  interval), 
design  or  specify  a  control  law  and  calculate  a  time  domain  response. 
An  earlier  version  of  this  chapter  has  appeared  in  the  literature  [9]. 

In  Chapter  Three,  several  existing  approaches  to  model  reduc¬ 
tion  for  continuous- time  models  are  extended  for  use  in  reducing  the 
order  of  discrete-time  models.  The  resulting  reduced  order  discrete¬ 
time  models  calculated  using  the  modal  approach  from  a  high  order 
continuous- time  model  along  two  paths  of  Figure  1.1  (Path  one  followed 
by  path  three  and  path  two  followed  by  path  four)  are  compared.  An 
earlier  version  of  this  paper  appeared  in  reference  [  26]  . 

The  subject  of  model  reduction  using  least  squares  to  fit 
the  response  of  the  reduced  order  model  to  that  of  the  high  order 
model  is  discussed  in  Chapter  Four,  where  an  alternate  method  of 
obtaining  the  data  used  by  least  squares  is  presented. 

A  modal  approach  to  control  law  reduction  (path  eight  of 
Figure  1.1)  is  presented  in  Chapter  Five.  An  earlier  version  of  this 
paper  appeared  in  reference  [27]. 

Chapter  Six  includes  several  aspects  of  both  model  reduction 
and  control  law  reduction.  The  modal  approach  to  model  reduction  is 
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discussed  in  more  detail  and  several  of  the  results  presented  in  the 
literature  are  compared.  Three  approaches  to  reduced  order  control 
law  design  are  discussed  and  applied  to  the  control  of  a  pilot  plant 
evaporator . 

Finally,  some  overall  conclusions  and  recommendations  for 
future  work  are  presented  in  Chapter  Seven. 
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CHAPTER  TWO 


DESCRIPTION  OF  A  COMPUTER  PACKAGE 
FOR  CONTROL  SYSTEM  DESIGN 


ABSTRACT 


A  GEneral,  Multipurpose  Simulation  and  COntrol  PackagE 
(GEMSCOPE)  is  described  which  assists  with  the  design,  analysis,  and 
^Sital  simulation  of  linear,  time-invariant  dynamic  systems.  The  pro 
cess  model  may  be  defined  as  transfer  functions  or  in  the  form  of  a 
standard  state  space  model  in  either  the  continuous -time  form  or  the 
discrete** time  form.  Options  are  included  for  generating  the  optimal 
control  matrices,  determining  the  state  feedback  matrix  for  non-inter¬ 
action,  reducing  the  order  of  the  model  and  of  the  control  law,  and 
calculating  the  open  or  closed- loop  system  response.  GEMSCOPE  was 
used  for  all  the  model  and  control  law  calculations  used  in  this 
thesis . 
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2.1  INTRODUCTION 


GEMSCOPE  (a  GEneral,  Multipurpose,  Simulation  and  COntrol 
PackagE)  is  an  integrated  series  of  computer  programs  to  aid  in  the 
simulation  of  linear  models  and  in  the  design  of  control  systems  for 
these  models.  It  has  been  designed  to  run,  either  as  a  batch  job  or 
from  a  time  sharing  terminal,  on  an  IBM  360  model  67  computer  operating 
under  the  University  of  Michigan,  time- sharing  terminal  oriented  system 
(MTS).  It  has  been  set  up  in  order  to  make  use  of  some  of  the  file  man¬ 
ipulation  programs  available  with  the  MTS  system.  Some  of  the  programs 
in  GEMSCOPE  were  used  in  recent  studies  by  Newell  [20,21]  and  most  of 
the  programs  were  used  in  the  design  and  simulation  calculations  to  be 
outlined  in  this  thesis. 

The  basic  data  used  in  GEMSCOPE  as  well  as  some  of  the  basic 
programs  (data  input  and  output,  calculation  of  a  state  space  model 
from  a  transfer  function  model,  calculation  of  a  discrete-time  model 
from  a  continuous- time  model,  calculation  of  a  time  domain  response) 
are  those  which  made  up  the  Continuous  Systems  Analysis  Program  (CSAP) 
described  by  Agostinis  and  Fisher  [,1,2,3].  GEMSCOPE  itself  is  des¬ 
cribed  in  a  series  of  research  reports.  The  GEMSCOPE  User's  Manual 
[29  ] describes  the  use  of  the  existing  programs  as  well  as  the  required 
data.  Other  reports  summarize  the  file  structure  and  the  subroutines 
used  in  GEMSCOPE  [30]  ,  indicate  the  required  steps  to  add  a  new  pro¬ 
gram  to  GEMSCOPE  [31]  ,  and  outline  a  sample  problem  using  GEMSCOPE  [29]. 

Many  other  digital  simulation  and  design  programs  exist.  One 
of  the  first  was  developed  by  Kalman  and  Englar  [15]  and  others  are  des¬ 
cribed  in  textbooks  by  Franks  [ 9]  and  Melsa  [  19]  .  Several  companies 


■ 


and  organizations  also  have  their  own  programs  for  internal  use.  This 
chapter,  therefore,  describes  an  example  of  such  programs,  not  a  unique 
program,  by  discussing  the  basic  philosphy  used  to  set  up  GEMSCOPE  and 
briefly  outlining  the  main  programs  in  the  package. 

2.2  GEMSCOPE  PROGRAMS 

GEMSCOPE  is  a  series  of  separate  programs  which  share  common 
disk  storage  files  and  can  be  executed  either  singly  or  sequentially  to 
perform  operations  such  as  model  definition,  control  calculations  or 
simulation  of  time  domain  responses.  From  the  point  of  view  of  simula¬ 
tion  the  important  programs  are: 

1.  Specification  of  the  programs  to  be  executed 

2.  Definition  of  the  state  space  model 

3.  Calculation  of  the  state  difference  equation 

4.  Calculation  of  the  time  domain  response. 

These  programs  are  shown  in  the  center  column  of  Figure  2.1.  However, 
the  usefulness  of  GEMSCOPE  is  increased  significantly  by  the  functions 
represented  by  the  blocks  in  the  right  most  column  of  Figure  2.1  which 
includes : 

5.  Generation  of  an  equivalent  state  space  form  for  models  entered 
as  transfer  functions 

6.  Reduction  of  the  order  of  the  state  space  model 

7.  Generation  of  state  feedback  required  for  non-interaction 

8.  Calculation  of  the  control  matrices  for  optimal  regulatory 
control 

9.  Simplification  of  high  order  control  laws 


i  r. 


. 


FIGURE  2.1  SIMPLIFIED  FLOW  DIAGRAM  OF  GEMSCOPE 


- 
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10.  Utilization  of  a  discrete  control  vector  for  open-loop  or 
optimal  state-driving. 

These  methods  are  discussed  in  the  following  sections.  Some  of  the 
"utility  functions"  for  performing  tasks  like  listing  output  data  are 
shown  in  Figure  2.1  but  are  not  discussed  further. 

The  GEMSCOPE  programs  in  their  present  form  consist  of  about 
7700  FORTRAN  source  cards  and  in  object  form  occupy  about  140,000  words 
of  storage.  Core  memory  requirements  vary  because  the  IBM  360/67  uses 
dynamic  memory  paging  but  a  typical  average  would  be  30,000  words. 
Central  Processing  Unit  time  requirements  vary  with  each  run  but 
average  about  40  seconds. 

2.2.1  Job  Specification 

It  should  be  emphasized  that  each  program  in  GEMSCOPE  is  a 
separate  mainline  program  and  can  be  executed  individually  or  in  any 
sequence  that  the  user  specifies.  The  mainline  sequencing  executive 
permits  the  user  to  specify  which  program  blocks  he  wants  executed  by 
simply  entering  a  series  of  integer  code  numbers  either  via  punched 
cards  (if  he  is  running  GEMSCOPE  as  a  batch  program)  or  via  a  key¬ 
board  on  a  time  sharing  terminal.  The  executive  then  generates  the 
appropriate  series  of  job  control  statements  and  stores  them  in  a  file. 
Special  files  of  job  control  statements  can  be  retained  for  frequently 
used  combinations  of  programs.  The  computer  monitor  system  will  then 
execute  the  GEMSCOPE  programs  in  accordance  with  the  job  control  state¬ 
ments  in  this  file  and  then  exit,  or  return  to  the  executive.  This 
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approach  maximizes  flexibility,  facilitates  reassignment  of  input- 
output  devices  and  files  for  each  program  block,  and  makes  it  a  simple 
matter  (since  each  block  is  a  self  contained  mainline  program)  to  add, 
modify  or  delete  programs.  This  approach  also  allowed  the  full  use  of 
the  computer  system  support  programs,  such  as  file  manipulation  rou¬ 
tines,  which  are  specified  by  the  use  of  job  control  statements. 

2.2.2  State  Space  Model  Specification 

The  system  to  be  simulated  can  be  specified  in  standard  state 
space  form  (order  <30)  or,  as  described  later,  in  the  classical  transfer 
function  and  block  diagram  form.  For  matrix  input  the  coefficient 
matrices  are  entered  by  the  user  or  obtained  from  previously  established 
input  files  in  the  form  of  a  continuous -time  system  as, 


x  =  Ax  +  Bu 


(2.1) 


(2,2) 


or  in  the  form  of  a  discrete- time  system  as  Equation  (2.2)  and 


x(j  +  1)  =  £x(j)  +  £u(j) 


(2.3) 


where:  x  =  state  vector 


u  =  input  vector 


=  output  vector 


w  =  input  plus  delayed  vector 

A  ^  A 

A,  C,  E  >  ij  A  =  constant  coefficient  matrices  of  approp- 


tiy  tt  y 


riate  dimension 


j  is  a  counter  for  time  intervals  such  that  x(j)  denotes 


x(t)  at  t  =  jT,  where  T  is  the  discrete  time  interval. 


■ 
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2.2.3  State  Difference  Equation 

When  a  state  space  model  in  the  continuous- time  form  of 
Equation  (2.1)  is  available,  the  equivalent  discrete- time  model  of 
Equation  (2.3)  may  be  needed  for  some  applications.  Thus,  a  program 
is  available  to  calculate  Equation  (2.3)  from  Equation  (2.1). 

Equation  (2.3)  agrees  exactly  with  Equation  (2.1)  at  each 

A 

sampling  point  if  the  input  vector,  u,  is  constant  over  the  sampling 
interval  (which  is  true  when  ^  originates  from  a  digital  control  com¬ 
puter),  otherwise  an  approximation  is  involved.  This  discrete  form 
is  particularly  convenient  for  digital  calculation  and  <t>  and  A  are 
available  for  use  by  other  programs. 

The  fundamental  matrix,  0,  can  be  evaluated  using  either  a 
power  series  expansion  or  the  Cayley -Hamilton  technique  as  described 
by  Agostinis  [1]  or  in  textbooks  such  as  Ogata  [25]  and  Saucedo  and 
Schiring  [28].  Implementation  of  the  latter  technique  includes 

algorithms  for  determining  real,  complex  and/or  repeated  eigenvalues  < 

A 

of  non- symmetric  A  matrices.  The  A  matrix  is  obtained  by  numeric 

A 

integration  of  0B  over  the  interval  T. 

This  program  to  calculate  the  discrete-time  model  was  orig¬ 
inally  part  of  CSAP  [1,2,3]. 

The  specification  of  the  time  interval,  T,  has  a  pronounced 
effect  on  the  numerical  stability  of  some  of  the  algorithms,  is  usually 
the  dominant  factor  affecting  simulation  time  and  determines  the  resol¬ 
ution  or  precision  of  the  time  domain  response.  Therefore,  to  provide 
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additional  flexibility,  a  program  has  been  included  to  modify  the  pre¬ 
viously  calculated  coefficient  matrices  in  Equation  (2.3)  so  that  they 
apply  to  an  interval  mT  (where  m  is  an  integer)  instead  of  T  and 

0(mT)  =  0m(T)  (2.4) 

A(mT)  =  j^0m_1(T)  +  .  .  .  +  0(T)  +  ij  A(T)  (2.5) 

In  general,  since  use  of  Equations  (2.4)  and  (2.5)  is  consid¬ 
erably  faster  than  recalculating  0  and  A,  it  is  most  expedient  to  specify 
a  small  interval  T  for  the  first  evaluation  of  0  and  A.  Relatively 
large  multiples,  mT,  can  be  used  for  exploratory  studies  and  then  m 
can  be  reduced  until  satisfactory  results  are  obtained. 

2.2.4  Time  Domain  Response 

The  time  domain  response  is  evaluated  directly  using 
Equations  (2.2)  and  (2.3).  The  initial  state  vector,  if  not  supplied 

by  the  user,  is  initialized  to  zero.  The  elements  of  the  input  vector, 

A 

u,  can  be  defined  by  selecting  standard  options  as  step,  ramp,  sinusoid, 
or  pulse  functions,  or  defined  by  a  discrete- time  series  of  point  values. 
The  control  variables  can  also  be  specified  as  originating  from  a  control 
law  as  discussed  in  the  section  on  optimal  multivariable  control  (Section 
2.2.8).  For  each  time  domain  response,  this  program  also  calculates  the 
numerical  value  of  the  quadratic  performance  index,  J,  which  is  optimized 
in  the  calculation  of  the  optimal  multivariable  control  matrices.  This 
performance  index  is  also  discussed  in  Section  2.2.8. 


. 
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The  basic  form  of  this  program  was  part  of  GSAP  [1,2,3]  and 
has  been  revised  slightly  for  use  in  GEMSCOPE. 

2-2.5  Transfer  Function  Model 

Most  of  the  current  research  work  in  areas  such  as  optimal 
control  is  performed  using  state  space  techniques.  However,  a  great 
deal  of  practical  design  and  analysis  of  control  systems  is  still  done 
using  classical  techniques  based  on  transfer  functions  and  block  dia¬ 
grams.  This  approach  has  been  the  subject  of  many  textbooks  such  as 
that  by  Coughanowr  and  Koppel  [5].  Therefore,  to  serve  as  a  convenient 
link  between  the  two  approaches,  a  program  was  written  which  would 
accept  a  problem  defined  in  terms  of  transfer  functions  and  a  coded 
equivalent  of  the  user's  block  diagram,  and  produce  an  equivalent  state 
space  representation  in  the  form  of  Equations  (2.1)  and  (2.2).  There 
are  no  arbitrary  restrictions  on  the  form  of  the  block  diagram  except 
the  dimensions  specified  for  vectors  and  matrices  in  the  digital  pro¬ 
gram.  Pure  time  delays,  isolated  differentiation  terms  and  isolated 
integral  terms  can  be  handled.  This  program  also  can  specify  control 
constants  for  standard  single  loop  control  laws  using  the  transfer 
functions,  the  block  diagram  and  several  different  criteria  such  as 
the  Ziegler-Nichols  frequency  response  technique,  the  Cohen  and  Coon 
response  curve  technique  and  others.  The  program  then  rearranges 
these  control  constants  into  a  matrix  form  so  that  it  can  be 

used  in  the  same  form  as  the  optimal  multivariable  control  law  dis¬ 
cussed  in  Section  2.2.8. 

This  program  was  originally  part  of  CSAP  [l,2,3]  and 
added  to  GEMSCOPE  by  Park  [27]. 


was 
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2.2.6  Model  Reduction 

Model  reduction  is  frequently  desirable  to  reduce  the  compu¬ 
tational  requirements  when  a  large  series  of  simulations  is  to  be  done, 
or  to  eliminate  states  which  are  not  of  interest  and/or  are  not  phys¬ 
ically  measurable.  Also,  when  the  model,  or  control  strategies  based 
on  it,  are  to  be  implemented  on  a  computer  controlled  process,  then 
model  reduction  will  reduce  the  amount  of  real-time  calculations. 
Therefore,  two  programs  have  been  included  in  GEMSCOPE  to  reduce  the 
order  of  a  state  space  model  by  calculating  a  reduced  order  model  whose 
state  vector  is  a  subset  of  the  vector  x.  Each  of  these  programs  permit 
the  user  to  reduce  the  order  of  either  the  continuous- time  model  of 
Equation  (2.1)  or  the  discrete-time  model  of  Equation  (2.3). 

In  the  first  program,  model  reduction  is  accomplished  using  a 
modal  analysis  of  the  high  order  model  by  eliminating  the  least  signif¬ 
icant  eigenvalues.  The  user  has  the  option  of  specifying  the  continuous¬ 
time  reduction  approaches  of  Marshall  [18],  Davison  [6],  Fossard  [8]  and 
Graham  [12]  as  well  as  the  equivalent  discrete- time  extensions  of  these 
approaches  which  are  discussed  in  Chapter  Three.  It  is  also  shown  in 
Chapter  Three,  that  using  these  modal  approaches,  the  continuous- time 
model  reduction  followed  by  conversion  to  the  discrete- time  form  is 
equivalent  to  discretization  of  the  high  order  model  followed  by  dis¬ 
crete-time  model  reduction..  Model  reduction  using  the  modal  analysis  is 
discussed  in  detail  in  Chapters  Three  and  Six. 

In  the  second  program,  a  least  squares  analysis,  and  a  sequence 
of  uniformly  distributed  random  numbers  are  used  to  calculate  the  reduced 
order  model.  Least  squares  is  used  to  "fit"  a  reduced  order  model  to 
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the  high  order  model  using  data  which  are  calculated  from  the  high  order 
model  and  the  random  numbers.  This  approach  is  discussed  in  detail  in 
Chapter  Four. 

2.2.7  Non-interaction 

For  the  control  of  multivariable  processes  it  is  frequently 
not  the  dynamic  relationship  between  a  given  pair  of  input  and  output 
variables  that  causes  difficulty  but  rather  the  interaction  between 
these  variables  and  other  state  variables.  One  design  approach  is  to 
develop  an  augmented  plant  model  such  that  there  is  no  interaction 
between  the  new  reference  variables  and  the  output  variables. 

The  decoupling  program  included  in  GEMSCOPE  is  a  modified 
version  of  one  developed  by  Gilbert  and  Pivnichny  [11].  The  program 
first  tests  to  see  whether  the  model  meets  the  necessary  and  sufficient 
conditions  for  non-interaction  by  state  feedback  as  described  by  Falb 
and  Wolovich  [7].  If  the  model,  given  in  the  form  of  Equations  (2.1) 
and  (2.2)  can  be  made  non-interacting  the  program  synthesizes  a  control 
law  in  the  following  form: 


u  =  Fx  +  Nv  (2.6) 

The  reference  input  vector  v  is  related  to  the  original  process  outputs 
by  the  following  equations  written  in  terms  of  Laplace  transforms 

;y(s)  =  W(s)  v(s)  (2.7) 

The  transfer  matrix  W(s)  is  diagonal  and  the  elements  are 
specified  by  the  user,  as  functions  of  two  vectors,  so  that  the  sys¬ 
tem  defined  by  Equation  (2.7)  has  the  desired  dynamic  response.  This 


fixes  the  elements  of  F  and  N.  The  basic  theory  behind  the  algorithm 
has  been  described  by  Gilbert  [10].  The  user  then  specifies  the  aug¬ 
mented  model  in  the  state  space  form  of  Equation  (2.1)  and  (2.2)  and 
re-enters  GEMSCOPE  for  further  simulation  and/or  control. 

This  approach  to  decoupling  control  has  been  applied  to  a 
fifth  order  distillation  column  model  by  Newell  [20,21]. 

2.2.8  Optimal  Multivariable  Controllers 

Optimal  multivariable  controllers  with  proportional-plus-in¬ 
tegral  feedback,  feedforward  compensation  for  constant  disturbances, 

d,  and  provision  for  driving  the  process  from  its  current  position  to 

SP 

setpoint  values,  ,  can  be  implemented  using  a  control  law  of  the  form 


j-1 

u(j)  =  Krax(j)  +  K1  23  y(i)  +  K^dU)  +  KSPZSP 

i=0 

where  u  and  d  are  defined  as  partitions  of  vector  ^  as 


(2.8) 


A 

u 


(2.9) 


and  K 


are  the  controller  matrices. 


For  the  special 


case  of  diagonal  control  matrices.  Equation  (2.8)  represents  a  set  of 


conventional  single  variable  control  laws. 


A  program  is  included  in  GEMSCOPE  to  allow  the  user  to  supply 
these  controller  matrices  or  to  calculate  them  for  the  partitioned  form 
of  Equation  (2.3)  as 
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x(j  +  1)  -  4>x(j)  +  A  u(j)  +0  d(j)  (2.10) 

A 

where  A  =  (A,0).  This  calculation  uses  recursive  relations  developed 
by  applying  the  concepts  of  discrete  dynamic  programming  to  the  model 

of  Equation  (2.10)  and  a  quadratic  performance  criterion  of  the  follow¬ 
ing  form: 

J  =  0N(x( N)  -  xSV  S  (x(N)  -  xSP) 

N  (2.11) 

+  0  |  (x(i)  -  xSP)T  Q  (x(i)  “  xSP)  +  uT(i-l)Ru(i-l) 

i=l 

where  /3,  a  time  weighting  factor,  and  Q,  R,  and  S,  state  or  control 
weighting  matrices,  are  specified  by  the  user.  The  dynamic  program¬ 
ming  solution  calculates  the  proportional  feedback,  feedforward,  and 
setpoint  controller  matrices  directly  from  the  model  in  Equation  (2.10), 
and  evaluates  the  integral  matrix,  from  an  augmented  state  space  model 
where  the  original  state  vector,  x  is  augmented  by  xP  where  x^  is 
defined  by 

j 

2SI(j)  =  T  ^  y(i)  (2.12) 

i=0 

and  where  y(i)  is  defined  by  Equation  (2.2)  with  E*  =  0.  Newell  et  al, 
[20,22,23]  showed  that  this  method  of  generating  control  matrices  is 
relatively  straightforward,  practical,  easy  to  implement  and  has  given 
outstanding  results  in  both  simulated  and  experimental  control  studies. 
Optimal  control  laws,  calculated  by  this  method,  are  used  in  Chapters 
Five  and  Six. 

2.2.9  Control  Law  Reduction 


When  optimal  control  theory  is  used  to  obtain  a  multivariable 


G-j 
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control  law  for  a  linear  state  space  model  in  the  form  of  Equation  (2.1) 
or  (2.3),  the  resulting  control  law  is  a  function  of  the  entire  state 
vector.  This  is  shown  by  the  control  law  in  Equation  (2.8).  In  many 
cases,  the  state  vector  contains  so  many  elements  that  the  resulting 
control  law  is  too  expensive  to  apply  to  the  control  of  the  actual 
process  using  a  control  computer.  This  expense  is  due  to  the  large 
on-line  computer  time  required  to  measure  all  the  elements  of  this 
large  state  vector,  or  to  calculate  any  which  cannot  be  measured. 

In  order  to  reduce  the  complexity  of  the  control  law,  GEMSCOPE 
contains  two  programs  to  calculate  low  order  control  laws  which  are  a 
function  of  a  subset  of  the  original  state  vector.  Each  of  these  pro¬ 
grams  make  use  of  the  high  order  control  law.  The  resulting  reduced 
order  control  law,  being  a  function  of  fewer  state  variables,  requires 
less  on-line  computer  time  for  measurement,  or  estimation,  of  the 
required  process  variables.  They  also  include  the  proportional  feed¬ 
back,  integral  feedback,  feedforward,  and/or  setpoint  control  modes 
which  were  present  in  the  high  order  control  law. 

The  first  program  uses  a  modal  analysis  of  the  high  order 
system  to  approximate  the  states  which  are  not  wanted  in  the  reduced 
order  control  law,  as  a  function  of  the  process  inputs  and  the  states 
to  be  retained  in  the  reduced  order  control  law.  This  approach  is 
developed  in  detail  in  Chapter  Five.  In  the  second  program,  a  least 
squares  analysis  is  used  to  fit  a  reduced  order  control  law  to  the 
high  order  control  law  using  data  which  are  calculated  from  the  high 
order  control  law  and  a  sequence  of  uniformly  distributed  random 
bers.  This  approach  is  discussed  in  Chapter  Six. 
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A  third  approach  for  calculating  a  control  law  of  reduced 
order  makes  use  of  two  of  the  GEMSCOPE  programs  already  described. 

A  reduced  order  model  can  be  obtained  using  the  programs  outlined 
in  Section  2.2.6  and  the  optimal  control  law  for  this  model  can  be 
calculated,  as  discussed  in  Section  2.2.8. 

Each  of  these  three  approaches  to  reduced  order  control  law 
design  are  discussed  in  Chapter  Six,  where  it  is  shown  that  the  ap-, 
proach  using  the  modal  analysis  to  reduce  the  high  order  control  law 

gives  the  best  results  when  applied  to  the  control  of  the  pilot  plant 
evaporator  outlined  in  Appendix  A. 

2‘2-10  Utilization  of  a  Discrete  Control  Vector 

Provision  has  been  made  within  GEMSCOPE  to  carry  out  time 
domain  simulation  using  control  variables,  u,  specified  by  a  dis¬ 
crete  time  series  stored  in  a  computer  file.  The  values  for  u 
can  be  supplied  by  the  user  or  calculated  by  another  program.  Nieman 
[24]  developed  programs  which  use  quasilinearization  and/or  linear 
programming  techniques  to  generate  the  optimal  control  policy  to  drive 
a  process  from  one  state  to  another  in  such  a  manner  as  to  minimize 
a  criterion  such  as  minimum  time,  or  minimum  sum  of  the  absolute 
errors.  These  programs  can  be  used  to  generate  a  discrete  series  of 
control  values,  which  can  be  read  into  GEMSCOPE  and  used  to  generate 
a  time  domain  response. 

2.3  DISCUSSION 

The  programs  discussed  above  form  a  useful  package  of  simu¬ 
lation  and  design  programs  for  linear  systems.  However,  there  is 
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still  a  need  for  further  development  of  control  and  numerical  tech¬ 
niques  that  will  make  the  package  more  flexible  and  will  extend  the 
applications  to  larger  systems,  non-linear  systems,  stochastic  sys¬ 
tems,  etc.  Some  possible  future  additions  to  GEMSCOPE  include: 

a)  Programs  to  design  state  vector  observers  and  filters.  This 
was  studied  by  Hamilton  [13]  and  could  include  the  Luenberger 
observer  [16]  and  the  Kalman  Filter  [14]. 

b)  Altering  the  time  domain  response  program  in  conjunction  with 
the  design  of  state  estimators,  to  include  the  provisions  for 
filtering  measurements  and  for  adding  random  noise  to  inputs 
and/or  states  so  noisy  processes  could  be  simulated. 

c)  Programs  to  include  time  delays  in  the  multivariable  feedback 
control  law,  as  was  studied  by  Alevisakis  [4]. 

d)  Programs  to  calculate  an  adaptive  control  algorithm  along  with 
the  calculation  of  the  time  domain  response  of  the  system.  One 
such  algorithm  was  considered  by  Oliver  [26]. 

e)  Programs  to  design  multivariable  control  laws  by  techniques 
other  than  the  dynamic  programming  approach  described  above. 
Some  possible  schemes  include  closed- loop  pole  placement  and 
the  newer  frequency  domain  techniques  described  by  MacFarlane 
[17]. 

f)  A  general  numerical  linearization  program  which  could  calculate 
the  linear  state  space  model  in  the  form  of  Equation  (2.1) 
from  a  non-linear  model  provided  by  the  user  in  the  form  of  a 


subroutine. 


* 


2.4  CONCLUSIONS 


GEMSCOPE  forms  the  base  of  an  easily  expandable  package  of  sirri' 
ulation  and  design  programs  for  linear  systems.  It  is  designed  to  fac¬ 
ilitate  the  addition,  modification  or  deletion  of  individual  programs 
and  to  make  full  use  of  the  support  programs  supplied  as  part  of  the 

computer  operating  system.  Several  possible  additions  to  GEMSCOPE  have 
been  discussed. 


CHAPTER  THREE 


MODEL  REDUCTION  FOR  DISCRETE-TIME  DYNAMIC  SYSTEMS 


ABSTRACT 

The  model  reduction  techniques  of  Marshall,  Davison  and  Fossard 
linear  continuous— time  models  are  extended  for  use  in  reducing  the 
order  of  linear  discrete-time  models.  Two  approaches  for  the  reduction 
of  high  order  continuous— time  models  to  low  order  discrete-time  models 
are  presented  and  evaluated.  The  resulting  reduced  order,  discrete- 
time  models  obtained  by  the  two  approaches  are  shown  to  be  equivalent. 

A  simple  numerical  example  is  solved  to  demonstrate  the  techniques 
presented. 
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3 . 1  INTRODUCTION 


Linear  state  space  models  have  been  extensively  used  to  math¬ 
ematically  describe  a  wide  variety  of  physical  systems,  such  as  air¬ 
craft,  chemical  plants,  and  refineries.  In  many  situations,  dynamic 
models  consisting  of  a  large  number  of  differential  or  difference 
equations  can  be  derived  from  theoretical  considerations.  Often, 
such  models  are  so  large  as  to  be  inconvenient  or  impractical  for  many 
purposes,  including  simulation  and  control  system  design  and  implemen¬ 
tation.  Consequently,  a  need  exists  for  systematic  procedures  for 
deriving  reduced  order  dynamic  models  from  high  order  models. 

Model  reduction  techniques  for  continuous- time,  state  space 
models  have  been  the  subject  of  many  recent  investigations  with  modal 
approaches  by  Marshall  [14],  Davison  [9],  Chidambara  and  Davison  [5, , 
6,7],  Graham  [12]  and  Fossard  [10]  receiving  the  most  attention.  By 
contrast,  only  a  few  investigations  have  considered  the  model  reduc¬ 
tion  of  discrete-time  systems.  Anderson  [4]  presented  a  method  for 
deriving  a  discrete- time  reduced  order  system  based  on  a  least 
squares  fit  of  the  response  of  the  reduced  order  model  to  that  of 
the  original  high  order  model.  He  has  also  presented  revisions  to  his 
method  to  reduce  computational  difficulties  (Anderson  [3]),  to  account 
for  pre-specified  elements  of  the  reduced  order  system  (Anderson  [2]) 
and  to  produce  better  agreement  at  the  final  steady  state  (Anderson 
[1]).  In  this  investigation  the  modal  methods  of  Marshall  [ 14] 

Davison  [9]  and  Fossard  [10]  are  extended  to  discrete-time  models. 

The  discrete-time  model  reduction  techniques  will  be  derived 
m  Section  3.2.  In  Section  3.3,  the  problem  of  reducing  a  high  order 


. 

am 
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continuous- time  model  to  a  low  order  discrete- time  model  is  considered. 
Finally ,  a  numerical  example  illustrating  the  model  reduction  techniques 
is  presented  in  Section  3.4. 
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3.2  MODEL  REDUCTION  FOR  DISCRETE-TIME  SYSTEMS 


3.2.1  Problem  Formulation 

Consider  the  stable,  time- invariant,  discrete- time  model; 

A  A 

x(j+l)  =  <£  x(j)  +  A  u(j)  (3.1) 


which  can  be  partitioned  as: 


2S1(j+1) 

~*i 

ti 

2Ex(j) 

rA 

Si 

= 

+ 

A 

x2(j+1)_ 

_=3 

i2(j)- 

_  =  2- 

(3.2) 


where  x  =  n-dimensional  state  vector. 

=  ^-dimensional  vector  of  states  to  be  retained  in  the  reduced 
order  model. 

*2  =  (n-/) -dimensional  vector  of  states  to  be  eliminated. 

A 

u  =  (nH-q) -dimensional  vector  of  inputs  which  will  include  the 

m-dimensional  control  vector  and  the  q-dimensional  disturbance 

vector . 

.  J  A 

<P  and  A  =  constant  coefficient  matrices  of  appropriate  dimensions, 
j  =  a  counter  for  time  intervals  such  that  x(j)  denotes  x(t)  at 
t  =  jT,  where  T  is  the  discrete  time  interval. 


The  objective  is  to  determine  a  reduced  order  model  of 
specified  order,  /,  which  is  of  the  form 

2^0+1)  =  =R  ~i  (d  >  + 


(3.3) 


1~' 
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where  0R  and  are  constant  matrices  of  appropriate  dimensions  and 
must  be  determined. 

Since  the  model  reduction  techniques  to  be  considered  are  based 
on  a  modal  analysis,  the  required  preliminary  material  will  be  presented. 

As  n  x  n  matrix,  M,  exists  which  transforms  Equation  (3.1) 
into  the  Jordan  canonical  form  (Gantmacher  [ 11 ]  ).  Define, 

x(j)  =  M  z(j)  ('3.4') 

or  in  partitioned  form, 

25L(j) 

2S20) 

where  z  is  the  canonical  state  vector,  z±  is  an  i- dimensional  vector, 
and  z ^  is  an  (n-/) -dimensional  vector.  It  is  assumed  that  matrix  M 
is  arranged  so  that  its  columns  are  ordered,  from  left  to  right,  in 
order  of  decreasing  significance  of  the  corresponding  eigenvalues  of  <fi. 
(This  can  always  be  achieved  by  appropriate  column  interchange.)  Since 
the  system  in  Equation  (3.1)  is  assumed  to  be  stable,  all  the  eigenvalues 
of  0  lie  in  the  unit  circle  (Kuo,  [13]).  Thus,  the  columns  of  M  on  the 
left  correspond  to  the  eigenvalues  of  0  which  affect  x  the  most,  and 
those  on  the  right  correspond  to  eigenvalues  of  0  which  affect  x  the 
least.  If  matrix  <£  has  n  distinct  eigenvalues,  the  columns  of  M  are 
then  the  eigenvectors  of  <P  and  M  is  referred  to  as  the  Modal  Matrix. 

When  Equation  (3.5)  is  substituted  into  Equation  (3.2),  the  Jordan  can¬ 


?L  =2 

=3 


2X(j) 

z2(j) 


(3.5) 


onical  form  results: 


■ 


. 


I 


29 


z.1(j+l) 
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HR 
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L_j. 
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r  A  -i 

k 

1 
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HR 
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l2(j)_ 

+ 

x. 

u(j) 


(3.6) 


A 


where  a,  a  block  diagonal  matrix,  and  6,  are  defined  as 


a  = 


a.  0 
=  1  = 


£'  22. 


=  M*1  0  M 


(3.7) 


t- 


rA 

ii 

f2 


(3.8) 


and  V  is  defined  as 


V  = 


V  V 
=  1  =2 


V  V 

L=3  =4. 


=  M 


-1 


(3.9) 


Because  of  the  way  the  columns  of  M  have  been  arranged,  the  Jor¬ 
dan  canonical  matrix,  a,  has  the  eigenvalues  of  <f>  which  most  affect  x^ 
the  matrix,  o^,  and  those  which  least  affect  x^  in  the  n -ji  x  n-j^ 
matrix,  a Hence,  it  follows  from  Equation  (3.6),  that  z^  represents 
the  modes  which  most  affect  x  ,  and  z^  represents  the  modes  which  least 
affect  x-^*  It  is  also  apparent  from  Equation  (3.6)  that  the  modes  rep¬ 
resented  by  _z ^  and  z2  are  noninteracting.  For  systems  with  n  linearly 
independent  eigenvectors,  a  ±s  a  diagonal  matrix  and  z±  and  z£  will 
always  be  noninteracting.  However,  for  systems  with  less  than  n  linearly 
independent  eigenvectors,  the  partitioning  of  the  original  system  in 
Equation  (3.2)  must  be  such  that  it  does  not  split  a  Jordan  block  in  a, 
so  that  z;^  and  z^  will  remain  non-interacting. 


. 


3.2.2  Marshall's  Analysis 

In  Marshall's  analysis  for  continuous- time  systems  (Marshall 
[14]),  the  reduced  model  is  derived  by  neglecting  the  dynamics  of  the 
modes  of  the  original  high  order  model  represented  by  _z  2 .  This  approach 
retains  the  significant  system  modes  and  offers  the  advantages  of  sim¬ 
plicity  and  zero  steady  state  error  for  constant  inputs.  Graham  [ 12 ]  ha 
shown  that  the  model  reduction  method  developed  independently  by 

Chidambara  (Chidambara  and  Davis  on  [6])  is  equivalent  to  Marshall's 
method  (Marshall  [14]).  In  this  section,  Marshall's  analysis  will  be 
extended  to  the  reduction  of  discrete-time  models. 

Equation  (3.2)  implies  the  relation, 

%(j+l)  =  ti  +  £2  -2(j)  +  !l  (3.10) 

An  expression  for  x2(j)  can  be  obtained  by  combining  Equations  (3.4) 
and  (3.9)  to  give 

*2(j)  =  I^~l  (z2(j)  "  V3  xx(j))  (3.11) 

As  an  approximation,  assume  that  z^Cj)  immediately  attains  the  new 
steady  state  value  (i.e.,  instantaneous  response).  Since  z ^  does  not 
affect  x^  significantly ,  the  effect  of  this  approximation  on  x-^  is 
small.  Thus,  from  Equation  (3.6), 

z2(j+l)  =  z2(j)  =  (I-a2)_1  |2  u(j)  (3.12) 

Now,  combining  Equations  (3.10),  (3.11)  and  (3.12)  gives  the  reduced 
order  system  of  Equation  (3.3)  with, 
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*k  ‘  *i  -  fc  S*'1  S3  (3'13) 

(3a4) 

Thus,  Marshall's  analysis  for  continuous- time  models  can  be 
extended  to  provide  a  straightforward  method  for  reducing  the  order  of 
discrete-time  models.  For  constant  inputs,  the  reduced  and  original 
systems  reach  the  same  steady  state  value  of  x  .  For  other  types  of 
sustained  inputs  (e.g.,  a  ramp  or  a  sinusoid)  there  is  a  continuing 
error  between  the  responses  of  the  reduced  and  original  systems 
(Graham  [12]).  However,  this  error  is  often  small. 

3.2.3  Davison's  Analysis 

Davison's  analysis  for  continuous- time  systems  (Davison  [9]) 
is  based  on  the  assumption  that  the  contribution  of  the  insignificant 
eigenvalues  of  the  original  system  to  the  system  response  is  small  and 
may  be  neglected.  In  Davison's  analysis,  a  reduced  system  results  in 
which  the  dominant  modes  are  excited  in  the  same  proportion  as  in  the 
original  system. 


The  solution  to  Equation  (3.1)  can  be  written  as, 

Ij 

£^~k  §  u(k) 

k=0 


(3.15) 


where  x(°)  is  the  initial  condition  for  vector  x- 
consider  the  case  with  x(0)  =  0.  Then, 


E,j-k  A  A  . 
£  A  u(k) 

k=0 


Following  Davison  ['9 ] , 


x(j+1)  = 


(3.16) 


. 


32 


Using  the  theorem  in  Appendix  D.l  and  Equation  (3.7),  Equation  (3.16) 
can  be  written  as. 


x(j+l)  =  M  ^  M_1  A  u(k) 

k=0 

Equation  (3.17)  becomes,  after  partitioning, 

0 


x^j+l) 


x2(j+1) 


=  1  =2 

M  M 
=3  =4 


e 

k=0 


2i"k 


o  ?J-k 


Si  h 
h  ^4 


AA, 

A  u(k) 


(3.17) 


(3.18) 


As  an  approximation,  assume  that  the  eigenvalues  in  a 2  have  no  effect 
on  the  system  response,  and  therefore  can  be  neglected  (i.e.,  a 2  =  0) . 
Equation  (3.18)  becomes. 


x1(j+l)  M1 
x2(j+l)  M3 

Define  an  /x  1  vector  _£  as 

j 

i-E 

k=0 


j 

E 

k=0 


«fk(Vx,  V2)  A  u(k) 


^(Vi,  V2)  A  u(k) 


(3.19) 


and  Equation  (3.19)  becomes, 


x-j^j+l) 

II 

Si 

x2^+1^ 

l_ _ 

Rearranging  Equation  (3.20)  gives 


(3.20) 


i  -  M^1  x^j+l) 


and 


-1 


x2(j+D  =  ^  Mx  x^j+1) 


(3.21) 
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As  in  Davison  [9],  and  are  calculated  by  the  following  analysis. 
Equation  (3.21)  is  substituted  into  the  unforced  form  of  Equation  (3.10) 
(i»e. ,  u(j)  =  0)  to  give, 

x^j+l)  =  (0L  +  d>2  M3M1"1)  x1(jl 

so  that,  =  +  £2  =3=i_1  (3.22) 


A 


To  calculate  the  solutions  of  the  original  model  and  the 

reduced  order  model  are  compared  with  x(0)  =  0. 


The  solution  of  the  original  model  in  Equation  (3.1)  for  x 
ignoring  a, ^  is, 

x^j+l)  =  mJ  £ 


k=0 


+  ^(k) 


(3.23) 


The  solution  of  the  reduced  order  model  is, 

j 


*i(j+D  -  £  J'k*  -■ 


k=0 


However,  using  the  equations,  <P  =  MaV  and  VM  =  I,  it  can  be 


shown  that 


h  =  Si  2i  Si 


-1 


so  this  reduced  order  solution  becomes, 


,  J 

J 

x1(j+l)  = 

^Mf-1  ^  u(k)  | 

k=0 

1) 

(3.24) 


Equating  Equations 


(3.23)  and  (3.24)  gives  the  following  expression  for 
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Ar  =  Mi  (V^!  +  V^)  =  Mx  |x  (3.25) 

Thus^  a  simple,  reduce:d  order,  discrete- time  system  has  been  derived 
using  Davison  s  modal  approach.  This  analysis  could  also  be  extended 
to  include  the  various  modifications  of  Davison's  Method  (Davison  [g]  , 
Chidambara  and  Davison  [6]). 


It  can  be  shown,  using  the  partitioned  form  of  the  expression 
~  that  the  transition  matrix  in  Equation  (3.22)  is  equivalent  to 
the  transition  matrix  in  Equation  (3.13),  obtained  by  Marshall's  analysis. 


3.2.4  Fossard's  Analysis 

Fossard  [ 10 ]  presented  a  revision  to  Davison's  continuous- time 
result  to  correct  the  steady  state  error  which  results  between  the 
reduced  order  model  and  the  original  high  order  model  after  a  step  change 
in  the  input.  This  correction  is  developed  here  for  the  discrete-time 
form  of  the  model. 

Consider  the  value  of  x^(j)  from  Equation  (3.5)  as 

x^j)  =  Mx  zx(j)  +  M2  z£(j)  (3.26) 

The  steady  state  values  of  z  and  z2  are  calculated  from  Equation  (3.6) 
as 


%ss  = 

ii 

U 

~ss 

(3.27) 

— 2ss  = 

A 

k 

A 

U 

“  ss 

(3.28) 

Thus,  the  steady  state  X-^  vector  can  be  expressed  as 

%ss  =  Si  w1  ti  is  +  m-^'1  &z  k. 


(3.29) 
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The  steady  state  vector,  xlss,  evaluated  using  Davison's 
resulting  reduced  order  model  can  be  obtained  in  a  similar  manner 


Thus, 


%(j)  =  z1(j) 


-1  £ 

A 

5. 

u 

=  1 

— ss 

. -1 

A 

A 

oc ) 

6n 

u 

=  1 

=  1 

— ss 

(3.30) 

(3.31) 

(3.32) 


To  have  the  steady  state  in  Equation  (3.32)  agree  with  that  in 
Equation  (3.29)  the  term 

M2  (J-ft)'1  i2  U(j) 

must  be  added  to  Equation  (3.32).  This  is  accomplished  in  the  following 
reduced  order  discrete- time  model  which  corresponds  to  the  continuous¬ 
time  result  presented  by  Fossard. 

2£^(j+l)  =0Rx^(j)  +|R^(j)  (3.33) 

%(j)  =  2^  (j)  +  1R  u(j)  (3.34) 

where 

Sr  =  =2  0-35) 

,  A 

=R  anc*  ~  t^ie  results  °f  Davison's  analysis  and 

are  given  in  Equations  (3.22)  and  (3.25) 
respectively. 

3-3  REDUCTION  OF  HIGH  ORDER  CONTINUOUS  MODELS  TO  LOW  ORDER 

DISCRETE  MODELS 


Models  of  physical  systems  are  often  represented  as  large 
of  linear  differential  equations,  in  the  general  form. 


sets 
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AA 

x  =  Ax  +  Bu  (3.36 

where  x  and  u  are  defined  in  Section  3.2.1,  and  A  and  B  are  matrices 
of  appropriate  dimensions.  For  both  simulation  studies  and  the  design 
and  implementation  of  computer  control  systems,  it  may  be  desirable  to 
have  available  a  low  order  discrete- time  model. 

Two  general  approaches  for  deriving  a  low  order  discrete- time 
model  from  a  high  order  continuous- time  model  are  shown  in  Figure  3.1. 
The  high  order  model  can  first  be  reduced  and  then  the  resulting  low 
order  model  discretized  (Approach  I),  or  alternatively,  the  high  order 
model  can  be  discretized  and  this  discrete- time  model  can  then  be 
reduced  (Approach  II).  The  continuous- time,  reduced  order  methods  of 
Marshall  [14],  Davison  [9]  and  Fossard  [10]  will  be  needed  for  Approach 
I  and  will  now  be  summarized. 


The  transformation  matrix  M  exists  (see  the  theorem  in  Appendix 
D.l)  which  transforms  Equation  (3.36)  into  the  Jordan  canonical  form  as 


where 


x  = 


Mz 


AA 

z  =  Jz  +  Gu 


A 

u 


(3.37) 


(3.38) 


(3.39) 


A 

G 


A 

VB 


(3.40) 


To  calculate  the  continuous- time  reduced  order  model,  the  high 
order  model  in  Equation  (3.26)  and  the  transformation  matrix  M  must  be 
arranged  and  partitioned  as  was  discussed  in  Section  3.2.1. 


37 


w 

g 

M 

H 


CO 

g 


M 

EH 

O  P 
U  W 


Pi 

w 


Q 

8 


5? 

W  W 
X  Eh 
W 
<3  Pi 
o 

O  CO 
g  M 
M  Q 
U 

P  Pi 
Q  W 

o 

Pi 

o 

pm 

co 

w 

e 


rs 

o 

p 


_  o 

<3  Eh 
O 

Pi  p 
PM  W 


p-l 

<3 


Q 

s 


CO 

§ 

6 

H 

Pm 


38 


The  reduced-order  model  is  of  the  form, 


A  A 

% '  tei +  te 


where  Marshall's  method  [14]  gives 


=  M ,  J  M  - 

=  R 

=1  =1  =1 

A 

A 

Bn 

=  B,  -  A„V, 

=R 

=1  =2=4 

-1 


-1  -1  A 

J„  G„ 


and  Davison's  method  [9]  yields, 


(3.41) 


(3.42) 

(3.43) 


h  ■  *1  +  =2  Ml 

A  A 

Sr  “  Si  gx 


-l 


Fossard's  method  produces  a  reduced  order  model  in  the  form 


A  A 

*£  -  Ar  +  Br  u 

%  =  H  +  Ir  s 


A 


where  AR  and  BR  are  given  in  Equations  (3.44)  and  (3.45)  and 


:R 


-1  A 

Ir  =2  =2  =2 


(3.44) 

(3.45) 

(3.46) 

(3.47) 

(3.48) 


It  is  shown  in  Appendix  B.3  that  the  discrete- time  form  of  E 

— R 

presented  in  Equation  (3.35)  is  identical  to  the  continuous- time  form 
presented  in  Equation  (3.48).  This  is  expected  since  Equations  (3.34) 
and  (3.47)  are  identical  in  form  and  must  hold  for  all  time. 

Approaches  I  and  II  differ  in  the  order  in  which  the  individ¬ 
ual  steps  of  model  reduction  and  discretization  are  carried  out.  From 
the  nature  of  these  operations,  intuitively  one  expects  that  both 
paths  will  result  in  an  identical  reduced  order  discrete- time  model. 


• 
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This  equivalence  is  shown  analytically  in  Appendices  B.l,  B.2  and  B.3 
by  making  use  of  the  definitions  of  <£  and  A  in  terms  of  A  and  B  as, 


^  =  exp  (AT) 


r T 

1  exp (A(T- 7) )Bdr 

Jo 

A 

and  BR  as, 

<fR  =  exp(ART) 

T  A 

exP  (AR(T-  T))BRdr 


(3.49) 


(3.50) 


The  equivalence  of  Approaches  I  and  II  will  also  be  demonstrated  by  the 
numerical  example  in  Section  3.4. 

Since  the  order  of  model  reduction  and  discretization  does  not 
affect  the  resulting  reduced  order  discrete-time  model,  one  should  con¬ 
sider  which  approach  is  the  most  convenient  for  a  particular  applica¬ 
tion.  If  the  only  requirement  is  the  derivation  of  a  low  order  discrete¬ 
time  model  from  the  original  high  order  continuous-time  model,  then 
Approach  I  should  be  followed,  since  its  computational  requirements 
are  less.  The  continuous- time  reduction  and  the  discrete- time  reduction 
involve  the  same  number  of  matrix  operations.  However,  the  derivation 
of  the  high  order  discrete- time  system  using  Equation  (3.49)  requires 
many  more  calculations  than  the  derivation  of  the  low  order  discrete¬ 
time  system  in  Equation  (3.50).  Thus,  Approach  I  is  computationally 
more  efficient  than  Approach  II.  However,  if  the  availability  of  a 
high  order,  discrete-time  model  is  desirable,  perhaps  for  the  evalua- 


■ 
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tion  of  proposed  control  systems,  Approach  II  should  be  used  since  the 
high  order,  discrete- time  model  must  be  derived  in  any  event. 

It  is  interesting  to  compare  the  results  of  the  discrete- time 
model  reduction  with  the  results  of  the  continuous- time  reduction.  It 
can  be  shown  by  elementary  matrix  algebra,  using  the  partitioned  forms 
of  the  expressions  <p  =  M  a  V  and  MV  =  I,  that  the  discrete  transition 
matrix  in  Equation  (3.13)  resulting  from  Marshall's  analysis  is 
equivalent  to  the  following, 

gR  -  ti  “  t2  X4"1  ^3  =  Mi  Si  Mi"1  (3.51) 

Thus,  by  comparing  Equation  (3.51)  with  (3.42),  Equation  (3.14)  with 
(3.43),  Equation  (3.22)  with  (3.44)  and  Equation  (3.25)  with  (3.45), 
it  can  be  seen  that  there  is  a  direct,  term-by-term  correspondence 
between  the  discrete-time  and  the  continuous-time  results. 

3 . 4  EXAMPLES 


3.4.1  Second  Order  Example 

A  simple  problem  considered  by  Chidambara  and  Davison  [6]  and 
by  Graham  [12],  will  be  used  to  demonstrate  both  the  discrete-time  model 
reduction  techniques  and  the  equivalence  of  the  two  approaches  shown  in 
Figure  3.1.  Consider  the  second  order  continuous- time  system  in 
Equation  (3.52),  which  is  to  be  reduced  to  a  first  order  discrete¬ 
time  system  in  x-i  • 


V 

’0  l" 

"0" 

X  = 

X2 

= 

-4  -5 

X 

+ 

1 

—  mm 

_  _ 

x 


u 


(3.52) 


. 


. 


, 
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The  eigenvalues  of  Equation  (3.52)  are  -1  and  -4  so  that. 


J  = 


-1  0 
0  -4 


and 


M  = 


1/3  -1/3 

n 

>i 

•N 

4  1 

-1/3  4/3 

1  1 

The  equivalent  discrete- time  system,  with  sampling  period  T  is 


x(j+l)  =  1/3 


4e 


-T  -4T  -T  -4T 


-e 


e  -e 


/  "T  "4T  “T  ,  -4T 
-4e  +4e  ,  -e  +4e 


x(j)  +  1/3 


3/4-e"T  +l/4e'4T 


-T  -4T 
e  -e 


u(j) 


(3.53) 


Equation  (3.53)  has  eigenvalues  of  e~T  and  e"4T. 


Marshall's  Analysis 

Marshall's  analysis  for  Approach  I,  using  Equations  (3.42)  and 
(3.43)  gives, 

£r  -  “i  jj.  sr1  ■  i 

A  A  -1  -1  A 

Sr  “  Si  -  A2V4  J2  g2  -  1/4 

Thus,  the  low  order  continuous -time  model  is 

=  -xL  +  u/4  (3.54) 

with  the  corresponding  discrete- time  model  as 

Xl(j+1)  =  e_T  xL(j)  +  (l-e'T)  u(j)/4  (3.55) 

For  Approach  II,  using  Equations  (3.13),  (3.14),  and  (3.53) 
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<t> 

=R 

=  tl 

-  <6  V 
=2=4 

-T 

=3  =  e  ^ 

A 

A 

-1 

(i-S 2^ 

=R 

"Si 

+  =2=4 

(l-e~T)/4 


and  the  reduced  order,  discrete- time  model  from  Approach  II  is, 

^(j+1)  =  e  T  X;L(j)  +  (l-e"T)u(j)/4  (3.56) 

which  is  identical  to  Equation  (3.55). 


Davison's  Analysis 

Davison's  analysis,  for  Approach  I,  using  Equations  (3.44)  and 
(3.45)  gives, 

h  ’  h  +  A 2  S3  M^1  =  -1 

A  A 

=R  =  ^1  Si  =  ^ 

The  low  order,  continuous- time  model  is, 

=  -x^  +  u/3  (3.57) 

with  the  corresponding  discrete-time  model  as 

x-L(j+l)  =  e  Tx1(j)  +  (l-e_i)  u(j)/3  (3.58) 

For  Approach  II,  using  Equations  (3.22),  (3.25),  and  (3.53) 

- 1  -T 

4>  —  <t>  +  d>  MM  —  e 

=R  =1  =2  =3=1 

Ar  *=  Mi  (V^  +  V^)  =  (l-e-T)/3 
and  the  reduced  order,  discrete- time  model  from  Approach  II  is, 

x-L(j+l)  =  e'T  xx(j)  +  (l-e“T)  u(j)/3  (3.59) 


. 
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which  is  identical  to  Equation  (3.58). 


Fossard's  Analysis 


Fossard ' 


s  analysis  produces  the  same  AR,  BR,  <£R  and  Ar  as  does 


Davison's  analysis  above. 


above.  ER  can  be  calculated  using  two  results.  The 


discrete- time  result  can  be  obtained  using  Equation  (3.35)  as 


1r  =  M2  (I-a2)_1  |2  =  -  1/12 


and  the  continuous- time  result,  Equation  (3.48),  produces 


These  are  identical  results  as  is  expected. 

3.4.2  Evaporator  Model 

As  a  check  of  the  computer  program,  the  tenth  order  continuous¬ 
time  evaporator  model  presented  in  Appendix  A  was  reduced  to  a  third 
order  discrete- time  model  by  the  two  approaches  of  Figure  3.1  using 
Marshall's  analysis.  In  both  cases  Wl,  W2  and  C2  were  the  states  retained 
in  the  reduced  order  model  and  x^  =  [Wl,  W2,  C 2 ]  .  The  eigenvalues 
retained  in  the  reduced  order  model  were  the  three  closest  to  zero  in 
the  continuous-time  case  and  the  three  closest  to  one  in  the  discrete¬ 
time  case.  Figure  3.2  shows  that  the  response  of  the  two  resulting  third 
order  models  are  identical  when  forced  by  a  +20%  change  in  feed  flow 
rate.  The  third  order  evaporator  models  calculated  by  Marshall's  analysis 
are  shown  in  Table  3.1,  where  the  similarity  of  the  two  discrete-time 
results  can  be  observed. 
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3.5  CONCLUSIONS 


Existing  model  reduction  techniques  of  Marshall,  Davison,  and 
Fossard  for  continuous- time  systems  have  been  extended  for  use  with 
discrete-time  systems.  Two  approaches  for  reducing  a  high  order 
continuous- time  model  to  a  low  order  discrete-time  model  have  been 
evaluated.  The  reduction  of  a  continuous- time  model  followed  by 
discretization  is  computationally  more  efficient,  but  analytically 
equivalent  to,  discretizing  the  high  order  system  and  then  reducing 
the  high  order  discrete- time  model  using  the  methods  outlined  in  this 
chapter . 


. 


CHAPTER  FOUR 


MODEL  REDUCTION  BY  LEAST  SQUARES  USING  RANDOM  DATA 

ABSTRACT 

A  least  squares  analysis  is  applied  to  model  reduction  using 
data  which  is  specified  by  the  high  order  model  and  a  sequence  of 
uniformly  distributed  random  numbers.  This  approach  produces  a 
reduced  order  model  that  is  better  than  the  model  calculated  by 
the  traditional  trajectory  fitting  approach  to  least  squares  when 
applied  to  the  reduction  of  a  tenth  order  evaporator  model.  This 
approach  also  eliminates  some  of  the  difficulties  found  in  the 
trajectory  fitting  application  of  least  squares  to  model  reduction. 

This  least  squares  approach  is  applied  to  the  reduction  of 
both  the  continuous-time  and  the  discrete- time  evaporator  models. 

The  problem  of  calculating  a  reduced  order  discrete^time 
model  from  a  high  order  continuous- time  model  is  considered.  It 
is  shown  that,  using  the  same  sequence  of  random  numbers,  different 
reduced  order  discrete-time  models  will  be  calculated  depending  upon 
whether  the  model  reduction  or  the  model  discretization  step  is 
calculated  first. 
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4.1  INTRODUCTION 

Many  physical  systems,  such  as  chemical  plants,  have  been 
mathematically  described  by  linear,  state  space  models,  which  con¬ 
sist  of  large  numbers  of  differential  or  difference  equations.  These 
models  often  contain  so  many  equations  that  they  cannot  economically 
be  used  for  simulation  studies  or  for  control  system  design.  Thus, 
a  need  exists  for  procedures  to  eliminate  some  of  the  less  important 
variables  from  the  high  order  model  to  obtain  a  dynamic  model  of 
reduced  order. 

4.1.1  Previous  Work 

Many  different  approaches  to  model  reduction  have  been  con¬ 
sidered.  In  the  derivation  of  most  complex  chemical  process  models, 
an  intuitive  reduction  procedure  is  applied  when  the  modelled  assumes 
that  a  particular  variable  responds  fast,  and  so  uses  a  steady  state 
relation.  This  approach  relies  upon  experience  with  the  particular 
process  being  modelled. 

Mathematical  approaches  to  model  reduction  using  a  modal 
analysis  and  a  least  squares  analysis  have  been  considered.  Marshall 
[13],  Davison  [6],  Fossard  [7]  and  others  used  a  modal  analysis  of 
the  high  order  continuous -time  model  in  order  to  calculate  reduced 
order  models.  These  methods  were  extended  for  use  with  discrete¬ 
time  models  in  Chapter  Three. 

As  an  alternative  mathematical  approach ,  Anderson  [4]  cal¬ 
culated  a  reduced  order,  discrete- time  model  using  a  least  squares 
analysis,  to  fit  the  response  of  the  reduced  order  model  to  that  of 


' 


. 

. 


the  original  high  order  model.  This  approach  was  extended  for  the 
reduction  of  continuous-time  models  by  Anderson  and  Kwan  [5].  Anderson 
[3]  reformulated  his  original  procedure  in  order  to  reduce  the  com¬ 
puter  storage  requirements.  Anderson  [2]  also  discussed  the  applica¬ 
tions  of  least  squares  to  model  reduction  when  some  of  the  elements 
in  the  reduced  order  model  matrices  can  be  prespecified.  In  addition, 
Anderson  [1]  considers  improving  the  agreement  between  the  high  order 
and  reduced  order  models  in  specific  regions  of  the  response,  by  addr¬ 
ing  extra  data  points  in  the  region  of  interest.  This,  however, 
would  be  at  the  expense  of  poorer  agreement  in  other  regions  of  the 
response.  Graham  et  al.  [8,9]  use  least  squares  along  with  quasi¬ 
linearization  for  the  reduction  of  high  order,  nonlinear  models. 

They  apply  a  parameter  significance  test  described  by  Rauch  [17] 
to  determine  any  of  the  unknown  parameters  which  can  be  neglected 
and  set  equal  to  zero.  Their  method  is  also  applicable  to  linear 
models . 

4.1.2  Difficulties  in  Previous  Least  Squares  Approach 

An  important  aspect  of  the  use  of  the  least  squares  approach 
to  model  reduction  is  the  type  of  data  which  is  used.  Anderson  [4] 
uses  data  obtained  from  a  time  domain  response  of  the  high  order 
model,  but  he  does  not  indicate  for  what  conditions  this  data  should 
be  generated.  As  was  discussed  by  Nicholson  and  Anderson  [15],  if  a 
time  domain  response  is  obtained  by  input  forcing,  the  inputs  must  be 
linearly  independent  to  avoid  a  linearly  dependent  set  of  equations 
in  the  least  squares  solution.  Thus,  a  step  in  each  input  cannot 
be  used.  Anderson  [15]  points  out  that  a  separate  reduced  order  model 


. 

■ 


. . .  while  such 


could  be  calculated  for  each  input,  but  adds  that  " 
results  are  not  meaningless,  they  are  neither  desirable  nor  use¬ 
ful,  McGinnis  [14]  took  this  approach  in  deriving  reduced 

order  models  of  a  distillation  column.  He  obtained  a  different 
reduced  order  (2  x  2)  transition  matrix  for  each  of  his  five  inputs. 
Then,  by  trial  and  error,  he  was  able  to  find  a  transition  matrix 
that  gave  satisfactory  results  for  all  his  inputs. 

If  a  reduced  order  model  is  calculated  for  good  agreement 
with  a  particular  response  caused  by  one  set  of  input  conditions,  it 
would  be  biased  towards  that  particular  set  of  conditions.  The 
response  of  the  resulting  reduced  order  model  would  then  have  poorer 
agreement  with  the  response  of  the  high  order  model  generated  by  any 
other  set  of  inputs. 

The  following  sections  of  this  chapter  summarize  the  least 
squares  approach  to  model  reduction  and  discuss  the  use  of  random 
numbers  as  data  for  the  least  squares  solution,  to  eliminate  the 
problems  discussed  in  this  section. 

4.2  LEAST  SQUARES  WITH  RANDOM  DATA 

The  fitting  technique  used  in  this  work  is  the  classical 
least  squares  technique  as  described  in  numerical  analysis  textbooks, 
for  example  by  Lapidus  [12]  and  as  used  for  model  reduction  by 
Anderson  [2  -  4].  The  only  difference  between  Anderson's  work  and 
that  presented  here,  is  the  method  of  specifying  the  data  which  is 
used  in  the  least  squares  fit.  The  least  squares  approach  will  be 
outlined  for  the  reduction  of  discrete-time  models  and  will  be  sum- 
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marized  in  Section  4.2.5  for  the  reduction  of  continuous -time  models. 

The  linear,  time- invariant  discrete-time  model  is  in  the 

form 

x(j  +  1)  =  0  x(j)  +  A  u(j)  (4.1) 

where  x  =  n-dimensional  state  vector 

A 

u  =  (m  +  q) -dimensional  input  vector  which  includes  the  m- 
dimensional  control  vector  and  the  q-dimensional  dis¬ 
turbance  vector. 

A 

0,  A  =  constant  coefficient  matrices  of  appropriate  dimen¬ 
sions  . 


The  state  vector  can  be  partitioned  as 


(4.2) 


where 


—2 


/-dimensional  vector  of  states  to  be  retained  in  the 
reduced  order  model 

(n-/) -dimensional  vector  of  the  remaining  states. 


The  desired  reduced  order  model  is  in  the  form  of 

^(j  +  1)  =  X-l  ( j  )  (4.3) 

where  0^  and  A^  are  constant  coefficient  matrices  of  appropriate 

dimensions.  This  section  will  discuss  the  evaluation  of  0  and  ^ 

— R  — R 

under  various  conditions.  Sections  4.2.1,  4.2.2  and  4.2.3  will  deal 
with  the  case  when  the  high  order  system  is  stable.  This  is  the  case 


when  the  eigenvalues  of  <£  lie  between  aero  and  one.  Section  4.2.4 
will  discuss  the  case  where  the  model  contains  integrating  states  cor¬ 
responding  to  eigenvalues  of  one.  Section  4.2.5  outlines  the  results 
for  continuous-time  systems  while  Section  4.2.6  discusses  other  appli¬ 
cations  of  the  theory. 


4.2.1  Calculation  of  for  a  Stable  System 

To  calculate  consider  the  following  unforced,  partitioned 

forms  of  Equations  (4.1)  and  (4.3) 


Xx(j  +  1) 

ii  tz 

x2(j  +  1) 

_=3  =4_ 

(4.4) 


and 


Xx(j  +  1)  =  0R 


(4.5) 


An  expression  for  can  be  written  as  the  least  squares  solution  of 


where 


W  =  Z  <£r 


(4.6) 


That  is 


<at 

IR 


T  -IT 
(Z  Z)  Z  W 


(4.7) 


T 

w  (1) 

Ft  1 

x}a) 

w  = 

T 

w  (2) 

• 

and  Z  = 

T 

x{  (2) 

• 

• 

•  T 

w  (k)  _ 

(4.8) 


ttl 

and  where  k  is  the  number  of  data  vectors.  The  j  vector  w(j)  is 


defined  as 


. 
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w(j)  =  x^j  +  1)  =  x-j^Cj)  +  <£2  x2(j)  (4.9) 


The  solution  in  Equation  (4.7)  is  identical  to  the  result 
presented  by  Anderson  [4]  when  the  latter  is  applied  to  an  unforced 
model.  However,  the  method  of  specifying  the  data  which  is  repre¬ 
sented  by  Z  and  W  is  different.  Anderson  fitted  the  reduced  order 
model  to  an  open  loop  response  of  the  high  order  model.  He  did  not, 
however,  specify  the  disturbance  used  to  generate  this  response. 

Thus,  in  his  presentation,  the  vectors  which  make  up  the  rows  of 
T  T  T 

Z  (x^(l),  x-^(2),  x-^(3),  etc.)  are  consecutive  points  along  a  partic¬ 
ular  open-loop  response.  In  the  work  outlined  here,  the  elements  of 
the  vectors  x^  ( j  )  and  x2  (j  )  ,  but  not  x-^j  +  1),  used  in  Equations  (4.8) 
and  (4.9),  are  uniformly  distributed  random  numbers.  This  was  done 
in  order  to  overcome  some  of  the  difficulties  with  the  trajectory 
fitting  approach.  With  the  use  of  random  numbers,  it  is  not  nec¬ 
essary  to  obtain  a  linearly  independent  set  of  inputs  to  use  in  the 
calculation  of  the  time  domain  response.  Also,  using  the  random 
numbers,  data  from  the  entire  range  of  possible  states  is  used.  This 
will  eliminate  the  bias  towards  a  particular  set  of  inputs  which  is 
found  in  the  trajectory  fitting  approach. 


4.2.2  Calculation  of  A  for  a  Stable  System 

Once  0  has  been  obtained,  ^  is  calculated  to  provide 
steady  state  agreement  between  the  high  order  model  and  the  reduced 
order  model  for  step  changes  in  u.  This  results  in  the  following 


expression 


, 
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=  "  £r>  [(I  -  0)  (4.10) 

where  the  notation  [E]^  represents  a  matrix  formed  by  the  first  £  rows 
of  E. 


4.2.3  Calculation  of  Both  0^  and  A^  Using  Least  Squares 

An  alternate  approach  for  calculating  a  reduced  order  model 
using  least  squares,  is  to  evaluate  both  0_  and  A„  at  the  same  time. 

To  do  this,  the  general  approach  outlined  in  Section  4.2.1  is  followed, 
except  that  both  0^  and  A^  are  calculated  by  Equation  (4. 7)  as 


where 


Z  = 


xJCD,  uT(l) 

x*(2),  4T(2) 


‘T...  AT.. 
X-L  GO  ,  u  (k) 


and  the  rows  of  W  are  calculated  using 


(4.11) 


(4.12) 


w(j)  =  xx(j  +1)  =  xx(j)  +  ^  x2(j)  +  ^  ft(j)  (4.13) 

Here,  all  the  elements  of  an<^  —  are  sPecified  as  random 

numbers . 


The  resulting  reduced  order  model,  as  calculated  by  Equa¬ 
tion  (4.11),  is  not  guaranteed  to  provide  steady  state  agreement 
between  the  high  order  and  the  reduced  order  models.  A  reduced 
order  model  with  better  steady  state  agreement  can  be  calculated  by 


. 
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adding  steady  state  points  to  matrices  Z  and  W,  as  is  discussed  by 
Anderson  [1],  but  this  will  be  at  the  expense  of  poorer  transient 
agreement  between  the  two  models.  A  second  approach  to  provide  steady 
state  agreement  is  to  recalculate^  using  Equation  (4.10). 

4.2.4  Reduction  of  Systems  with  Integrating  States 

For  systems  with  at  least  one  integrating  state  (with  at 
least  one  eigenvalue  of  <£  equal  to  one),  the  model  reduction  approach 
outlined  in  Sections  4.2.1  and  4.2.2  must  be  revised  slightly.  The 
reduced  order  state  vector  x.  can  be  partitioned  as 

*-£) 

corresponding  to  the  integrating,  X.JJ  and  non-integrating,  x^j  state, 

variables  of  x  .  Matrix  <£  can  be  similarly  partitioned  as 
1  — -K 


(4.15) 


The  partition  corresponding  to  the  integrating  states  can  be  pre¬ 
specified  as 


(4.16) 


The  dimension  of  the  identity  matrix  in  Equation  (4.16)  is  the  same 
as  the  number  of  integrating  states.  The  partition  of  0R  corres¬ 
ponding  to  the  non-integrating  states  is  calculated  using  least 
squares  as 


W 


(4.17) 
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In  this  case 


Z  = 


4(1) 

4(2> 

4(k) 


til 

and  the  j  row  of  W  is  defined  as 


(4.18) 


^(3)  =  Sin  Vj)  +  S2  *2(j) 


(4.19) 


where  0^  has  been  partitioned  as 


Si  =  (Su>  SiN) 


(4.20) 


corresponding  to  the  partitioning  of  x  in  Equation  (4.14).  The 
elements  of  vectors  x^  and  x^  are  specified  as  random  numbers. 

Once  0R  has  been  calculated,  AR  is  obtained,  by  columns 
for  agreement  at  large  times,  as 


where 


4i 

-R 


x^(j  +  1)  -  0R  x^(j) 


A 

u. 

l 


(4.21) 


£ 

A 

u. 

1 


X 


.th  -  -  A 

l  column  of  A^ 

.th  _  A 

i  element  of  u 


=  x^  for  large  values  of  the  counter  j  for  a  step  change 
in  u.  (u^  =  0,  k  7*  i) 


The  values  of  the  elements  of  x^ ( j  +  1)  and  x^(j)  corresponding  to 
the  non-integrating  states  are  equivalent  and  are  the  steady  state 
values  of  these  states  for  a  non-zero  value  of  u.  (u^  =  0,  k  ^  i). 


, 
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The  values  of  the  elements  of  (j  +  1)  corresponding  to  the  inte¬ 
grating  states  are  calculated  as 


x^(j  +1)  =  ti  2L*(j)  +  t2  -2^^  +  -l  “i  (4.22) 

where  A ^  =  i*^  column  of  A^  and  the  elements  of  jx^(j)  corresponding 
to  the  integrating  states  are  specified  by  random  numbers  and  the 
elements  of  x^O)  corresponding  to  the  non-integrating  states  are  the 

above  mentioned  steady  state  values  for  the  same  non-tzero  value  of 

A 
u. . 

l 

Equation  (4.21)  provides  for  steady  state  agreement  of  the 
non- integrating  states,  and  agreement  of  the  rate  of  change  of  the 
integrating  states  after  the  transients  die  out. 


4.2.5  Least  Squares  Reduction  of  Continuous -Time  Models 

The  above  procedures  can  be  applied  to  the  reduction  of 
continuous-time  models  in  the  form: 


to  give  a  reduced-order  model  in  the  form: 

A  A 

2£l  =  +  iR  “ 


(4.23) 


(4.24) 


For  the  case  of  a  stable  system,  with  all  the  eigenvalues  of  A 
negative 

=  (ZT  Z)"1  ZT  W  (4.25) 

where  Z  and  W  are  defined  as  in  Equation  (4.8)  except  that 
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T 

w  (j)  =  xx  =  A±  X-^j)  +  A2  x2(j)  (4.26) 


Also,  for  agreement  between  the  responses  of  the  high  order  and  the 
low  order  models  at  steady  state 


(4.27) 


For  the  reduction  of  continuous -time  models  which  contain 
integrating  states  (at  least  one  eigenvalue  of  A  equals  zero) ,  this 
procedure  can  be  modified  in  a  similar  manner  as  the  discrete-time 
procedure  was  modified  in  Section  4.2*4. 

4.2.6  Other  Applications  of  Least  Squares  with  Random  Data 

The  use  of  random  data  with  the  least  squares  analysis 
could  be  used  whenever  the  data  to  be  fitted  can  be  generated  from  a 
mathematical  relationship.  For  example,  this  would  be  the  case  if 
a  polynomial  approximation  to  a  part  of  a  trigonometric  function 
was  needed.  This  approach  is  described  in  Chapter  Six,  where  it  is 
used  to  generate  incomplete  state  feedback  control  laws  when  the 
more  complex,  high  order  control  law  is  known. 

4.3  APPLICATION  TO  THE  TENTH  ORDER  EVAPORATOR  MODEL 

The  least  squares  reduction  method  using  random  data,  as 
described  above,  has  been  applied  to  the  reduction  of  the  tenth 
order  evaporator  model  in  its  open-loop  form,  as  presented  in 
Appendix  A,  and  in  a  closed-loop  form,  with  Bl  controlling  Wl  and 
B2  controlling  W2.  The  closed-loop  form  was  obtained  using 


■ 


B1  =  0.9203  W1 


B2  =  2.768  W2 

These  two  feedback  relations  change  the  unity  eigenvalues  of  and 

change  W1  and  W2  from  integrating  states  to  non-integrating  states. 
This  closed-loop  model  is  also  shown  in  Appendix  A.  All  the  reduced 

order  models  are  presented  in  Appendix  F. 

4.3.1  Computational  Aspects 

The  random  numbers  required  for  the  approach  of  this  chap¬ 
ter  were  calculated  using  the  IBM  Scientific  Subroutine  Package  (SSP) 
program  PANDU,  as  described  by  IBM  [10,11],  The  numbers  calculated 
by  RANDU  are  referred  to  as  random  numbers  even  though  they  are 
really  only  pseudo-random,  having  been  calculated  by  a  mathematical 
algorithm  and  not  by  a  truly  random  process.  However,  since  these 
numbers  are  generated  systematically  using  a  computer  routine,  it 
is  possible  to  duplicate  any  sequence  for  purposes  of  comparison. 

The  random  numbers  were  uniformly  distributed  between  -1  and  +1, 
which  corresponds  to  process  values  of  from  zero  to  twice  the 
process  steady  state,  since  the  variables  in  the  model  are  in 
normalized  perturbation  form. 

The  computer  storage  requirements  for  the  least  squares 
program  was  reduced  by  losing  the  partitioned  approach  to  the  evalu¬ 
ation  of  Equation  (4.7)  suggested  by  Anderson  [3].  In  this  approach 
matrices  Z  and  W  are  partitioned  as 
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l  Z  " 

~«r 

z  = 

S2 

and 

W  = 

s2 

S3 

• 

s, 

• 

• 

• 

• 

• 

(4.28) 


where  the  are  square  matrices  of  order  H  and  the  W.  matrices  have 
H  rows  each.  Equation  (4.7)  can  now  be  written  as 


=i 


0  = 
■=R 


T  T  T 

h  h  +  hh  +  =3  h  + 


y  [i  i 


=1  +  h  +  S3  S3  + 


(4.29) 


This  reduces  the  computer  storage  requirement  by  a  great  deal  since 

each  of  the  terms  on  the  right  hand  side  of  Equation  (4.29),  ZT  Z. 

T 

and  W  can  be  evaluated  as  the  data  is  collected.  Furthermore, 

T 

a  0R  matrix  can  be  calculated  for  each  Z^  Z^  term  which  is  obtained. 
0R(j)  is  now  defined  to  include  j  terms  in  each  of  the  sums  in  Equa- 


tion  (4.29)  as 

"  j 

-1 

j 

E  si  s± 

E  si  si 

i=l 

H* 

II 

h- 1 

1 _ 

(4.30) 


Now,  new  values  of  0  are  calculated  until  0  (j)  and  0n (j  -  1) 

— K  — R  =R 

become  sufficiently  close  as 


-  1) 


<  r 


(4.31) 


where  f  is  a  convergence  criterion  parameter  specified  by  the  user 
and 
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(4.32) 
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is  an  element  of  E.  For  most 


where  E  is  any  (n  x  m)  matrix  and  e, 

ij 
_6 

models  calculated,  a  value  of  z,  =  10  was  used.  This  resulted  in 
values  of  k,  in  Equation  (4.8),  of  between  151  and  516. 

4.3.2  Evaporator  Model  Results 

The  method  of  least  squares  reduction  using  random  data 
has  been  applied  to  the  tenth  order  evaporator  model.  Each  of  the 
reduced  order  models  calculated  was  third  order  with  the  three  most 
important  states,  Wl,  W2  and  C2  in  their  state  vector.  The  results 
will  be  discussed  in  this  section. 

Unforced  Response  of  <£ 

—  IS. 

In  the  proposed  reduction  method,  only  d>  is  calculated 

—K 

using  least  squares.  Thus,  the  unforced,  open-loop  response  of  <£ 

— K 

to  a  non-zero  initial  state  should  be  compared  to  that  of  <£  for  the 
same  non-zero  initial  state.  This  comparison  is  shown  in  Figure  4.1 
for  the  initial  conditions  of  W1(0)  =  +20%  and  C2(0)  =  +15%.  Both 
holdups,  Wl  and  W2  remain  at  their  initial  values  (except  for  a  2% 
error  in  W2).  This  is  expected  since  these  are  integrating  states 
and  are  only  slightly  affected  by  these  initial  state  values.  The 
C2  responses  of  the  two  models  agree  very  well.  However,  the  respons 
of  C2  for  the  reduced  order  model  is  delayed  behind  that  of  the  tenth 
order  model.  This  is  not  observed  in  any  of  the  input  forcing  simu¬ 
lations  . 

Different  Random  Number  Sequences 

The  IBM,  SSP  routine,  RANDU,  that  was  used  to  calculate  the 
random  numbers,  requires  an  odd  integer  as  input  data  to  initialize 
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the  random  number  sequence.  Each  odd  integer  starts  a  different 
sequence.  Thus,  the  response  of  models  calculated  using  different 
number  sequences  must  be  compared.  Figure  4.2  shows  that  the  response 
of  three  different  reduced  order  models,  calculated  with  different 
number  sequences,  agree  very  well.  In  fact,  two  of  these  models 
produce  identical  responses  to  the  +20%  feed  flow  change.  Further¬ 
more,  all  three  responses  agree  well  with  the  response  of  the  tenth 
order  model  and  lead  the  tenth  order  response  as  would  be  expected. 

The  reduced  order  models  whose  responses  are  shown  in  Figure  4.2, 
are  open-loop  models  with  two  integrating  states  and  are  calculated 
using  the  approach  outlined  in  Section  4.2.4. 

Continuous  vs.  Discrete  Reduction 

A  reduced  order  discrete-time  model  can  be  calculated  from 
a  high  order  continuous-time  model  along  two  paths,  as  shown  in 
Figure  4.3.  Path  I  obtains  the  discrete-time  form  of  the  high  order 
model  and  reduces  this  using  the  discrete  model  reduction.  Path  II, 
reduces  the  order  of  the  continuous-time  model  and  then  discretizes 
the  resulting  reduced  order  continuous-time  model  to  get  the  desired 
reduced  order  discrete-time  model.  When  a  modal  approach  to  reduction 
is  used,  Chapter  Three  showed  that  the  results  by  Paths  I  and  II 
produce  an  identical  reduced  order  discrete-time  model.  This  is  not 
the  case  for  least  squares  reduction,  as  is  shown  in  Appendix  B.4. 
Figure  4.4  demonstrates  the  different  responses  for  the  two  reduced 
order  discrete- time  models  calculated  by  the  two  paths  of  Figure  4.3, 
using  the  same  random  number  sequence.  However,  it  should  be  pointed 
out  that  the  difference  between  the  two  reduced  order  responses  in 
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FIGURE  4.3  PATHS  FOR  CALCULATING  A  LOW  ORDER  DISCRETE-TIME 
MODEL  FROM  A  HIGH  ORDER  CONTINUOUS-TIME  MODEL 
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Figure  4.4  is  about  the  same  as  the  difference  due  to  different 
number  sequences  shown  in  Figure  4.2. 

Reduction  of  the  Closed-Loop  Model 

In  order  to  test  the  proposed  method  of  model  reduction  on 
a  model  of  a  different  form,  and  also  to  compare  it  with  other  least 
squares  approaches,  a  closed-loop  evaporator  model  was  formed  as 

noted  above.  This  model  did  not  have  any  integrating  states  since  the 
feedback  control  loops  added  changed  the  two  unity  eigenvalues  of  0 

to  lie  between  zero  and  one.  Figure  4.5  shows  that  the  proposed 
method  produces  a  model  whose  response  is  very  close  to  that  of  the 
tenth  order  model.  The  0^  matrix  for  this  third  order  model  was  cal¬ 
culated  by  Equation  (4.7)  and  the  ^  matrix  was  calculated  by  Equation 
(4.10).  This  model  also  agrees  with  the  tenth  order  model  at  steady 
state.  The  response  of  the  model  whose  0^  and  were  calculated 
together  using  least  squares  and  random  data  by  Equation  (4.11)  to 
(4.13)  and  whose  was  not  recalculated  for  steady  state  agreement, 
shows  an  expected  large  steady  state  error  in  Figure  4.5.  The  response 
of  the  model  for  which  0^  was  calculated  to  fit  an  open-loop,  unforced 
response  of  the  tenth  order  model  due  to  an  initial  non-zero  state  and 
^  was  calculated  for  steady  state  agreement,  shows  much  faster 
dynamics  than  do  any  of  the  other  models  shown  in  Figure  4.5.  The 

initial  state  of  the  trajectory,  to  which  0  of  the  last  model  was 

— R 

fitted,  was  chosen  by  random  selection.  This  may  have  caused  the 
difficulty  by  not  having  all  the  modes  of  0  sufficiently  excited. 
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4.4  DISCUSSION 


There  are  several  advantages  to  the  method  of  reduction 
using  least  squares  with  random  data. 

The  first  advantage  is  that  the  resulting  model  is  not 
biased  towards  any  particular  set  of  conditions.  Often,  when  a 
reduced  order  model  is  calculated,  using  least  squares,  to  fit  a  part¬ 
icular  response,  it  works  well  for  that  response,  but  not  particularly 

well  f°r  any  other  response.  An  example  is  the  model  whose  response 
is  shown  as  the  plus  signs  in  Figure  4.5.  It  gave  poor  agreement  with 
the  high  order  model  for  this  case  of  +20%  change  in  feed  flow.  How¬ 
ever,  as  Figure  4.6  shows,  it  agreed  well  with  the  high  order  model 
under  the  conditions  for  which  it  was  calculated. 

The  second  advantage  of  the  random  data  approach  is  its 
ease  of  application.  The  results  shown  in  this  chapter  that  were  cal¬ 
culated  using  this  method  were  obtained  with  one  computer  run  each. 
However,  for  the  one  model  calculated  by  fitting  the  unforced  response, 
several  computer  runs  were  required  before  it  was  determined  how  many 
points  along  the  response  should  be  considered.  The  convergence  test 
of  Equation  (4.31)  could  not  be  used  because  this  caused  the  use  of 
too  many  steady  state  points  and  the  corresponding  response  had  even 
faster  dynamics  than  the  response  shown  in  Figure  4.5  by  the  plus  signs. 

A  third  advantage  of  the  proposed  approach  is  in  the  con- 

T 

ditioning  of  the  matrix  (Z  Z)  which  must  be  inverted.  As  indicated 
by  Ralston  [16,  page  233],  a  matrix  is  ill-conditioned  if  the  inverse 
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of  the  matrix  (normalized  so  that  the  largest  element  is  one)  has 

very  large  elements  with  respect  to  one.  For  each  case  run  in  this 

T 

work,  the  inverse  of  this  normalized  form  of  (Z  Z)  was  calculated. 

For  all  cases  when  the  random  data  approach  was  used,  the  largest 

element  of  this  inverted  matrix  was  about  1.3,  the  same  order  of 

magnitude  as  the  largest  element  before  inversion.  On  the  other  hand, 

in  the  open-loop  response  fitting  cases  tried,  the  largest  element 

4 

of  this  inverted  matrix  was  at  least  of  order  10  ,  and  often  much 
higher.  Thus,  for  the  examples  used  in  this  work,  the  approach  using 
random  numbers  produced  a  much  better  conditioned  set  of  equations  to 
solve . 

4.5  CONCLUSIONS 


The  approach  to  least  squares  model  reduction  calculates 
better  reduced  order  models  of  the  tenth  order  evaporator  using  random 
data  than  using  data  from  an  open-loop  model  trajectory.  At  the 
same  time,  this  random  data  approach  is  easier  to  apply. 

When  a  reduced  order  discrete-time  model  is  calculated  from 
a  high  order  continuous-time  model,  different  results  are  obtained 
using  least  squares  depending  upon  whether  the  reduction  or  the 
discretization  step  is  performed  first. 
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CHAPTER  FIVE 


MODAL  APPROACH  TO  CONTROL  LAW  REDUCTION 

ABSTRACT 

A  method  is  proposed  for  reducing  the  complexity  of  state 
feedback  control  laws  for  high  order  systems.  Using  a  modal  approach, 
selected  state  variables  are  eliminated  to  give  an  incomplete  state 
feedback  control  law.  The  proposed  method  for  control  law  reduction  is 
applied  to  an  optimal  control  system  for  a  tenth  order  model  of  a  pilot 
scale  evaporator.  Simulation  and  experimental  results  show  that  this 
modal  approach  is  practical  and  results  in  better  control  than  methods 

based  on  reduced  order  models  or  conventional  multiloop  design  tech* 
niques . 
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5.1  INTRODUCTION 


The  object  in  this  investigation  is  to  develop  a  systematic 
approach  for  reducing  the  complexity  of  multivariable  controllers  for 
high  order  systems.  It  is  assumed  that  the  system  to  be  controlled  can 
be  adequately  described  by  a  linear,  time -invariant  state  space  model 
and  that  satisfactory  control  can  be  achieved  using  a  multivariable 
state  feedback  controller  which  may  also  include  integral  feedback, 
feedforward,  and  setpoint  control  modes.  Simplification  of  the  control 
law  is  achieved  by  eliminating  selected  state  variables  to  give  an 
incomplete  state  feedback  control  law.  The  motivation  for  this  design 
approach  is  that  in  most  practical  applications,  any  degradation  in 
system  performance  will  be  offset  by  the  practical  advantages  of  fewer 
measurements,  or  state  estimates,  being  required  and  reduced  demands  on 
the  on-line  computer.  The  simplified  control  law  is  determined,  a 
priori,  by  off-line  calculations. 

The  first  step  in  the  proposed  approach  is  to  select  a  control 
law  which  will  give  satisfactory  control  performance  when  applied  to  the 
best  available  model  of  the  process.  In  this  work,  the  control  law  is 
designed  by  applying  optimal  control  techniques  to  a  linearized  process 
model.  The  resulting  controller  is  a  linear,  time-invariant  feedback 
control  law  that  utilizes  the  entire  state  vector.  Simplification  of 
this  control  law  is  accomplished  by  eliminating  selected  state  variables 
from  the  feedback  control  terms.  The  problem  of  deriving  an  incomplete 
state  feedback  control  law  from  a  state  feedback  control  law  will  be 
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referred  to  in  the  study  as  "control  law  reduction",  and  the  resulting 
controller  will  be  referred  to  as  a  reduced  order  control  law. 

5.1.1  Previous  Work 

The  design  of  incomplete  state  feedback  and  output  feedback 
control  systems  for  linear  systems  has  been  the  subject  of  a  large  num¬ 
ber  of  investigations  in  recent  years.  Optimal  control  formulations  in 
which  the  control  configuration  is  constrained  to  be  a  function  of  the 
output  vector  has  received  considerable  attention  by  Davison  [2],  Kosut 
[6],  Levine  and  Athans  [7]  and  Sims  and  Helsa  [15].  However,  the  cal¬ 
culation  of  these  optimal  control  policies  often  requires  a  significant 
computational  effort,  particularly  for  large  systems.  Furthermore,  the 
dependence  of  the  optimal  control  policy  on  the  initial  state  must  be 
taken  into  consideration,  as  discussed  by  Kosut  [6]  and  Levine  and 
Athans  [  7]  . 

An  alternative  approach  considered  by  Anderson  [1]  and  Nichol¬ 
son  [11-133  is  to  derive  a  low  order  state  space  model  by  model  reduc¬ 
tion  techniques  and  then  use  the  optimal  control  policy  for  the  low 
order  model  as  a  suboptimal  control  law  for  the  high  order  system. 

Rogers  and  Sworder  [ 14]  have  proposed  that  the  reduced  order  model 
and  its  optimal  control  policy  be  calculated  simultaneously  with  the 
additional  restriction  that  the  low  order  control  law  is  the  best 
suboptimal  controller  for  the  high  order  system.  An  iterative  solution 
of  nonlinear  matrix  equations  is  required  in  their  approach. 

None  of  the  above  approaches  makes  direct  use  of  an  additional 
item  of  information,  namely,  the  optimal  control  law  for  the  high  order 
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system.  The  strategy  adopted  in  this  investigation  is  to  derive  an 
incomplete  state  feedback  law  from  a  state  feedback  control  law  by 
using  a  modal  analysis  of  the  high  order  state  space  model  to  elimin¬ 
ate  selected  state  variables. 

5.2  CONTROL  LAW  REDUCTION 

Consider  the  linear,  time -invariant,  discrete- time  model  in 

the  form: 

x(j+l)  =  0x(j>  +  Au(j)  +  Qd(j)  (5.1) 

Z(j)  =  C  x(j)  (5.2) 

where  x  is  the  n-dimensional  state  vector, 

u  is  the  m-dimensional  control  vector, 
d  is  the  q-dimensional  disturbance  vector, 

^  is  the  p-dimensional  output  vector, 

A,  @,  C  are  constant  coefficient  matrices  of  appropriate 
dimensions, 

j  is  a  counter  such  that  x(j)  denotes  x(t)  at  t  =  jT, 
where  T  is  the  discrete  time  interval. 

It  is  assumed  that  the  high  order  control  law  has  the  form: 

u(j)  =  x(j)  +  K1  ^  Z(i)  +  KFF  d(j)  +  KSP  y_S?  (5.3) 

i=0 

SP  SP 

where  £  is  the  constant  setpoint  vector  such  that  y^  is  the  desired 

value  of  y  ,  i  =  1,  ...,  p 
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£  *  £  >  E  >  S  are  the  proportional  feedback,  integral  feed¬ 

back,  feedforward  and  setpoint  control  matrices  of  appropriate 
dimensions . 


The  control  law  of  Equation  (5.3)  is  chosen,  not  only  because 
it  can  stabilize  the  system,  but  also  because  it  is  a  practical  control 
law  which  can  be  easily  applied  using  an  on-line  digital  computer  as  was 
shown  by  Newell  [10].  The  control  matrices  in  Equation  (5.3)  can  be 
determined  by  any  applicable  design  technique.  In  this  work  they  were 
calculated  by  dynamic  programming,  minimizing  a  quadratic  performance 
criterion,  as  described  by  Newell  [ 10 ] - 

The  objective  of  this  investigation  is  to  simplify  the  control 
law  by  reducing  the  number  of  state  variables  such  that 

u(j)  =  K™  x^j)  +  x(i)  +  K^F  d(j)  +  k|P  ySP;  (5.4) 

i=0 

where  x-^  is  the  ^-dimensional  subset  of  x  which  is  to  be  retained  in 
the  reduced  order  control  law.  Practical  considerations  dictate  that 
the  elements  of  X-^  be  chosen  so  that  the  system  defined  by  Equation  (5.1) 
can  be  stabilized  by  a  control  law  of  the  form  of  Equation  (5.4). 

Criteria  for  selecting  the  state  variables  to  be  retained  in  x^  will 
be  discussed  later. 


5.2.1  Method  Development 

If  the  state  vector,  x,  is  partitioned  such  that 


(5.5) 


•- 


. 
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then  the  state  feedback  control  law  of  Equation  (5.3)  can  be  expanded 
in  terms  of  and  x^}  to  give: 


u(  j) 


=  xx(j)  +  Kf  x20)  +  K1  ^  Z(i) 

i=0 


+  KFF  d(j)  +  KSP 


SP 

1 


(5.6) 


The  use  of  an  exact  expression  for  jc  in  Equation  (5.6),  derived  from 
Equations  (5.1)  and  (5.2),  would  result  in  a  complicated  control  law. 
However,  the  use  of  an  approximation  for  x^,  in  terms  of  the  other 
variables,  would  result  in  the  desired  form  of  the  control  law. 

Several  possible  approaches  to  this  approximation  exist: 

1.  It  could  be  assumed  that  x^  =  j),  which  simplifies  the  analysis 
but  is  too  extreme  an  assumption  for  most  systems. 

2.  A  second  assumption  is  that  the  response  of  x  is  instantaneous 
and  hence  x2(j)  is  an  algebraic  function  of  x-^(j),  u(j),  and 
d(j).  However,  it  is  not  difficult  to  find  systems  where  the 
response  time  of  x  (j)  is  so  slow  that  it  cannot  be  neglected. 

3.  A  third  approximation  for  x.^  could  be  derived  using  a  modal 
analysis  of  the  high  order  system  represented  by  Equation  (5.1). 


The  modal  analysis,  which  was  selected,  is  analogous  to  that 
used  in  the  reduction  of  continuous- time  models  originally  proposed  by 
Marshall  [9]  and  as  extended  to  discrete- time  systems  in  Chapter  Three. 
The  basic  steps  are  as  follows: 


An  n  x  n  matrix,  M,  exists,  which  transforms  Equation  (5.1) 
into  its  Jordan  canonical  form,  as  shown  by  Gantmacher  [3].  This  sim¬ 
ilarity  transformation,  which  also  preserves  the  eigenvalues  of  d>, 
is  defined  by 


■ 
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x  =  M  z 


(5.7) 


or  in  partitioned  form 


X  ~ 

M,  M_ 

z _ 

—1 

=  1  =2 

-1 

x„ 

H,  M, 

z 

-2 

=3  =4 

2  J 

(5.8) 


where  _z  is  the  canonical  state  vector  and  z^  and  z have  the  same 
dimensions  as  and  respectively.  When  Equation  (5.8)  is  sub¬ 
stituted  into  Equation  (5.1),  the  following  Jordan  canonical  form 
results  (shown  here  in  partitioned  form) 


(j+1) 

__ 

2l 

0  ■ 

I 

/—S 

•r-) 

l"! _ 

+ 

V 

u(j)  + 

V 

(j+2) 

0 

=2_ 

£2(  j) 

M. 

J  2. 

d(j) 


(5.9) 


where  a,  b_}  and  v  are  defined  as: 


a  =  M  1  <t>  M  = 


=i  s 


L  =  =2 


(5.10) 


5  =  M'1  A  = 


*>  =  m"1  e  = 


=i 


=2. 


2i 


l22. 


(5.11) 


(5.12) 


®  i-s  a  block  diagonal  matrix  with  the  eigenvalues  of  d>  on  its 
principal  diagonal.  It  then  follows  directly  from  Equation  (5.9)  that 
each  element  of  z  is  associated  with  a  particular  eigenvalue  of  <t> . 


An  expression  for  x2  can  be  calculated  using  Equation  (5.8) 
and  matrix  V  which  is  defined  as, 


■ 
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-1 


V  = 


M 


(5.13) 


to  give 


(5.14) 


Equation  (5.14)  could  be  used  to  eliminate  from  the  high  order  con¬ 


trol  law  in  Equation  (5.6)  .  However,  the  resulting  control  law  would 
then  contain  z^.  As  an  approximation,  the  dynamics  associated  with 
can  be  neglected  and  z^(j)  can  be  replaced  by  its  steady  state 
value.  It  will  be  shown  later  that  this  approximation  is  a  reasonable 
one  if  the  columns  of  M  are  arranged  in  an  appropriate  fashion.  Thus, 
for  constant  input  vectors,  assume  that  z has  reached  its  final  steady 
state  value;  then  from  Equation  (5.9) 


(5.15) 


Equation  (5.15)  is  then  used  with  Equation  (5.14)  to  give 


=  *2  SO)  +  (I-gjf1  22  -  V3  xx 


(j)  (5.16) 


The  desired  reduced  order  control  law  can  now  be  obtained  by  combining 
Equations  (5.6)  and  (5.16)  such  that  the  coefficient  matrices  in  Equa¬ 
tion  (5.4)  are  defined  as: 


(5.17) 


■ 


-c 


* 
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SP 

=R 


u  SP 
K  K 


where 


K 


u 


Thus ,  a  reduced  order  control  law  has  been  derived  in  which 
the  feedback  contribution  is  a  function  of  only  a  selected  subset  of 
the  state  vector,  rather  than  the  entire  state  vector.  It  should  be 
noted  that  the  control  system  designer  already  has  the  freedom  to 
decide  what  variables  are  included  in  the  feedforward,  integral  and 
setpoint  terms. 


5.2.2  Continuous-Time  Systems 

A  similar  analysis  could  be  used  to  obtain  a  continuous- time 
reduced  order  control  law.  Consider  the  following  linear,  time-invar¬ 
iant,  continuous- time  model: 


x=Ax+Bu+Dd  (5.18) 

1  =  C  x  (5.19) 


where  x,  u,  d,  y,  and  C  are  defined  earlier  and  A,  B,  and  D  are  con¬ 
stant  coefficient  matrices  of  appropriate  dimensions.  This  system 
can  be  stabilized  by  the  high  order  control  law: 


u  =  K 


FB 


'*S‘  S« 


..  ,  „FF  ,  ,  SP  SP 

y  dt  +  K  d  +  K  y 


(5.20) 


The  resulting  reduced,  continuous-time,  control  system  is  as  follows 


2  ■  SR®  — l  +  Sr 


s: 


1  dt  +  Kf  d  +  K^P  ySP 


(5.21) 


' 


' 
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where 


FB 

£r 


K“  (K™ 


sT  C  X3> 


I 

£r 


KU  K1 


FF  u  FF  FB  - 1  - 1 

Sr  =  S  <£  '  S2  h  S2  S2> 


and 


K®P  =  KU  KSP 
—  K.  =  = 


Ku  =  (I  +  &  v.-1  j'1  g0)_1 


=2  =4  =2  =2' 

and  G ,  and  are  defined  from  the  continuous- time  canonical  sys¬ 

tem  as: 


z. 

0  " 

z. 

G, 

~H  ~ 

-1 

=  1 

-1 

+ 

=  1 

u  + 

=  1 

z_ 

0 

Jo 

Gn 

Ho 

L-2  J 

=  2J 

L— 2  J 

L  =2  - 

L=2  J 

(5.22) 


J  = 


Si  2 


g  g2 


=  V  A  M 


G  = 


£1 


l£2 


=  V  B 


H  = 


H 


LH2 


=  V  D 


Matrices  M  and  V  are  the  same  as  those  used  for  the  discrete- time  system, 
and  can  be  calculated  using  either  the  continuous- time  or  the  discrete¬ 
time  models  (See  the  theorem  in  Appendix  D.l). 


' 
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This  continuous- time  control  law  is  subject  to  all  the  re¬ 
strictions  which  were  imposed  on  the  discrete- time  control  law  derived 
earlier . 

5.2.3  Application  Considerations 

In  applying  this  method,  two  important  design  decisions  must 
be  made:  first,  which  subset  of  x  should  be  included  in  the  reduced- 
order  control  law  as  'x  and,  secondly,  which  subset  of  the  system's 
eigenvalues  should  be  retained  in  a^. 

The  designer  must  first  specify  which  states  are  to  be  re¬ 
tained  in  the  control  law.  This  choice  may  require  a  trial  and  error 
approach,  but  a  few  general  guidelines  are  available.  For  example, 
process  variables  which  are  to  be  tightly  controlled  or  constrained 
must  be  included  in  the  reduced  order  control  law.  Other  process  var¬ 
iables,  which  if  not  controlled  would  give  poor  system  performance, 
must  also  be  included  in  the  reduced  order  control  law.  Furthermore, 
the  state  variables  retained  in  x^  must  be  able  to  satisfactorily  stab¬ 
ilize  the  high  order  system  by  means  of  Equation  (5.4).  A  further  con¬ 
sideration  in  the  choice  of  which  states  variables  to  retain  in  x-^  is 
whether  a  state  variable  is  measurable  or  not.  It  is  obvious  that  other 
considerations  aside,  those  state  variables  which  can  be  most  easily 
or  economically  measured  are  the  ones  which  should  be  retained.  Un¬ 
measurable  states  should  be  included  in  x-^  only  if  they  are  essential 
for  satisfactory  control  since  a  state  estimation  procedure,  such  as 
the  Kalman  filter  [5]  or  the  Luenberger  observer  [8]  would  then  be 
required . 
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After  the  elements  of  ^  have  been  selected,  the  designer  must 
assign  the  eigenvalues  of  g  to  either  ^  or  «2  by  an  appropriate  re¬ 
arrangement  of  the  columns  of  M  in  Equation  (5.7).  In  general,  the 
eigenvalues  which  have  the  most  significant  effect  on  the  elements  of 
-1  should  be  Placed  in  The  two  most  important  factors  to  consider 

m  assigning  these  eigenvalues,  as  shown  by  the  following  analysis, 
are  the  magnitude  of  the  eigenvalues  of  <P  and  the  magnitude  of  the 

elements  of  the  modal  matrix,  M.  Equation  (5.7)  can  be  expanded  for 
the  itb  element  of  x  as 

Xi  =  mil  Z1  +  mi2  Z2  +  +  min  zn  ~  h  2,  .  .  . ,  -n  (5.23) 

where  m±.  is  an  element  of  M  and  x.  and  z.  are  the  ith  element  of  x 

and  z  respectively.  The  magnitude  of  m  determines  the  effect  of  the 
• th  .th  1J 

J  mode  on  the  i  state  variable.  Similarly,  if  g  has  distinct 
eigenvalues,  then  the  response  of  the  kth  canonical  state  , of  the 
unforced  system  can  be  written,  from  Equation  (5.9),  as 

zk<j+1)  =  ak  zk(J)  k=  1,  2,  ...,  n  (5.24) 

where  is  an  eigenvalue  of  <t>. 

For  an  eigenvalue  close  to  the  origin.  Equation  (5.24)  indi¬ 
cates  that  the  corresponding  mode  decays  rapidly  and  hence  its  effect 
on  a  particular  element,  x^  of  the  state  vector,  would  be  short¬ 
lived.  Thus,  Equations  (5.23)  and  (5.24)  imply  that  the  kth  mode 
will  most  significantly  affect  a  state  variable,  x.,  when  m.,  is  large 
and  «k  is  also  large.  In  this  case,  a  should  be  included  in  a  . 

K  =1 

In  other  situations,  mik  will  be  very  small  and  thus  zfc  will  have  a 
negligible  effect  on  x^  In  other  situations,  an  eigenvalue  of  g 


. 


will  be  very  small  in  magnitude  and  will  have  only  a  short-lived  effect 
on  x.,  reaching  its  steady  state  value  quickly.  In  these  two  situa¬ 
tions,  the  eigenvalues  should  be  incorporated  into  However,  in 

other  cases,  the  decision  may  be  less  obvious  since  a  large  value  of 
mik  may  correspond  to  a  small  value  of  <*k,  or  vice  versa.  After  the 
eigenvalues  have  been  partitioned  into  a±  and  the  columns  of  M 
must  be  rearranged  accordingly. 

If  the  columns  of  M  are  arranged  as  above  so  that  the  modes 
which  most  affect  x^  ar^  in  z^,  then  the  modes  which  least  affect  x^ 
will  be  in  z2.  This  means  either  that  the  m  values  relating  z2  to 
-q  are  small,  or  that  the  eigenvalues  in  are  small.  Thus,  the 
steady  state  assumption  for  z2  will  not  significantly  affect  x  since 
the  canonical  states  in  z2  either  respond  faster  than  those  in  z or 
do  not  affect  x^  due  to  the  small  values  of  their  m_  elements.  Thus, 
the  magnitude  of  the  eigenvalues  in  a  and  the  associated  m  values 

~Z  ij  ’ 

give  the  designer  a  qualitative  feel  for  the  accuracy  of  the  reduced 
order  system.  Iterative  trial  and  error  solutions  to  the  design  prob¬ 
lem,  with  or  without  closed-loop  simulations  or  experimental  runs  may 
be  required  to  obtain  a  satisfactory  reduced  order  control  law. 

One  additional  restriction  must  be  imposed  for  systems  with 
repeated  eigenvalues  which  do  not  have  linearly  independent  eigenvec¬ 
tors.  The  block  diagonal  matrices,  and  must  be  assigned  so  that 
-1  and  —2  are  non- interacting  as  was  assumed  in  partitioning  Equation 
(5.9).  This  will  always  be  the  case  for  systems  with  n  linearly  inde¬ 
pendent  eigenvectors,  since  then  a  is  a  diagonal  matrix.  However,  for 
systems  with  less  than  n  linearly  independent  eigenvectors,  the  part- 


; 
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itioning  of  the  original  system  must  be  such  that  it  does  not  split  a 
Jordan  block  in  a,  so  that  z^  and  z ^  will  remain  non-interacting.  This 
might  necessitate  increasing  or  decreasing  the  size  of  vector  X;L 
slightly.  However,  if  the  original  system  cannot  be  satisfactorily 
partitioned  to  keep  z ^  and  z^  non-interacting,  then  the  parameters  in 
the  original  model  could  be  changed  by  a  small  amount  so  that  the 
eigenvectors  would  be  linearly  independent. 

5-3  APPLICATION  TO  THE  CONTROL  OF  AN  EVAPORATOR 

The  proposed  method  was  evaluated  by  comparing  simulated  and 
experimental  response  data  for  a  pilot  plant  evaporator  which  is  des¬ 
cribed  in  Appendix  A.  The  simulated  responses  were  generated  from  a 
10  order  state  space  model  of  the  evaporator  which  was  derived  by 
Newell  [10].  The  following  control  laws  were  compared: 

(a)  a  third  order  control  law  designed  by  the  proposed  method 
using  a  tenth  order  optimal  control  law  as  the  starting  point. 

(b)  a  third  order  control  law  which  is  optimal  for  a  third  order 
state  space  model  and  a  quadratic  performance  index.  The  third  order 
model  was  derived  from  the  tenth  order  model  usings  a  model  reduction 
technique  based  on  the  modal  analysis  of  Marshall  as  presented  in 
Chapter  Three. 

(c)  a  control  law  which  is  optimal  for  the  tenth  order  model  and 
a  quadratic  performance  index. 

(d)  a  conventional  multiloop  feedback  control  scheme  using  con¬ 
troller  constants  which  were  determined  experimentally  by  Jacobson  [4]  . 

The  first  step  in  the  application  of  the  proposed  method  is 
to  determine  which  states  to  include  in  x^.  An  analysis  of  the  evap- 
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orator  showed  that  three  state  variables,  product  concentration,  C2, 
and  the  two  holdups,  W1  and  W2  must  be  retained  in  C2  must  be 

retained  since  regulation  of  C2  is  specified  to  be -the  primary  control 
objective.  W1  and  W2  must  be  retained  since  these  are  integrating 
states  (corresponding  to  eigenvalues  of  £  of  one)  which  would  tend  to 
exceed  physical  operating  limits  if  they  were  not  controlled.  The  set- 
point  and  integral  control  modes  were  also  chosen  to  be  a  function  of 
these  three  state  variables,  making  ^  =  x^. 

The  second  step  in  the  proposed  method  is  to  choose  which 
eigenvalues  to  retain  in  a by  an  appropriate  arrangement  of  the 
columns  of  M  (M  is  shown  for  the  evaporator  model  in  Appendix  A) .  If 
the  previous  choice  of  2^  =  [Wl,  W2,  C2]  is  made,  then  the  three  lar¬ 
gest  eigenvalues  should  be  retained.  This  follows  from  an  inspection 
of  matrix  M  since  the  three  slowest  modes  have  the  most  effect  on  the 
elements  of  x^. 

Once  the  state  vector  and  eigenvalues  have  been  partitioned 
and  the  columns  of  M  have  been  arranged,  the  high  order  control  law, 
matrix  V,  and  the  canonical  state  equation  are  partitioned  as  shown 
in  Equations  (5.6),  (5.13),  and  (5.9),  respectively.  Then,  the 
reduced  order  control  matrices  can  be  calculated  using  Equation  (5.17). 


The  third  order  control  laws  which  were  designed  by  the  pro¬ 
posed  method  are  presented  in  Tables  5.1  and  5.2.  Table  5.3  summar* 
izes  the  values  of  the  quadratic  performance  index,  J, 


N 


J 


j  =  l 


CD  T  Cp  p 

(x(j)  -  x  )  g(x(j)  -  x  )  +  u  (j-1)  R  u(j-l) 


(5.25) 


•  *  ^ 


TABLE  5.1 


PROPORTIONAL,  FEEDFORWARD  AND  SETPOINT  CONTROL 

MATRICES  DESIGNED  BY  CONTROL  LAW  REDUCTION 


2.467 

4.288 

4.128 


0.02163  -4.705 
-1.340  8.885 
9.760  9.528 


1.238 

0.9815 

0.9978 


-0.5639 

0.2177 

0.9877 


-0.4138 

-0.001227 

-0.001398 


-2.468 

-4.289 

-4.128 


-0.02119  5.276 

1.340  -9.100 

-9.763  -10.52 


TABLE  5.2 

PROPORTIONAL  PLUS  INTEGRAL  CONTROL 
MATRICES  DESIGNED  BY  CONTROL  LAW  REDUCTION 


2.977 

5.066 

0.09369 

-1.246 

-5.408 

8.100 

5.686 

12.29 

11.58 

0.4117 

0.8398 

0.02088 

-0.2887 

-0.4861 

0.8159 

0.8559 

1.945 

1.068 

' 
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obtained  by  controlling  the  tenth  order  model  with  three  of  the  con¬ 
trol  laws  mentioned  above  (controllers  (a),  (b)  and  (c)).  The  optimal 
tenth  order  control  law  is  shown  in  Appendix  A  and  was  calculated 
using  q  =  diag  (0,0,10,0,0,0,10,100,0,0)  and  R  =  diag  (0.05,0.05,0.05). 
The  optimal  control  matrices  for  the  third  order  model  used  the  same 
R,  with  Q  -  diag  (10,10,100)  and  are  shown  in  Tables  5.4  and  5.5.  The 
simulated  runs  presented  in  Table  5.3  are: 

(a)  Proportional  feedback  control  with  a  non-zero  initial  condi¬ 
tion  (i.e.,  x(0)  ^  0), 

(b)  Proportional-integral  feedback  control  with  a  non-zero  initial 
condition, 

(c)  Proportional  feedback  -  feedforward  control  with  a  +20%  feed 
flow  disturbance  (starting  from  x(0)  =  0), 

(d)  Proportional  feedback-setpoint  control  with  a  +7%  setpoint 
change  in  C2. 

The  results  summarized  in  Table  5.3  show  that  the  third  order 
control  law  calculated  by  the  proposed  method  was  better,  as  indicated 
by  the  lower  value  of  J,  than  the  third  order  controller  calculated 
from  the  reduced  order  model.  The  only  exception  to  this  is  for  the 
case  of  proportional  feedback- feedforward  control  where  the  two 
third  order  control  laws  give  almost  identical  values,  differing  by 
only  0.25%.  In  each  case,  the  high  order  optimal  control  law  gave  the 
lowest  value  of  J,  as  would  be  expected. 

The  transient  responses  for  proportional  feedback  control  in 
Figure  5.1  demonstrate  that  the  reduced  order  control  law  calculated 
by  the  proposed  method  gives  better  control  than  the  control  law 


. 
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TABLE  5. A 


PROPORTIONAL,  FEEDFORWARD  AND  SETPOINT  CONTROL 


MATRICES  DESIGNED  FROM  REDUCED  ORDER  MODEL 


4.904 

5.784 

4.093 


-0.4013  -11.92 

-1.600  4.425 

9.685  9.357 


1.238 

0.9832 

0.9983 


-0.5583 

0.2231 

0.9937 


-0.4128 

-0.00054 

-0.00122 


-4.904 

-5.784 

-4.094 


0.4017  12.47 

1.599  -4.648 

-9.686  -10.35 


TABLE  5.5 

PROPORTIONAL  PLUS  INTEGRAL  CONTROL 
MATRICES  DESIGNED  FROM  REDUCED  ORDER  MODEL 


5.49 

6.429 

5.519 


-0.1903  -12.00 

-1.386  4.487 

12.26  11.81 


0.9893  -0.05066 

1.156  -0.3255 

0.8254  1.935 


-1.175 

0.4373 

1.090 
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designed  from  the  third  order  model.  This  is  especially  evident  for 
the  primary  controlled  variable  C2,  where  a  slight  oscillation  develops 
when  the  control  law  designed  from  the  reduced  order  model  is  used. 
Figure  5.2  presents  a  similar  comparison  for  proportional  plus  integral 
control  law  with  the  same  initial  state  as  was  shown  in  Figure  5.1. 
Again,  the  proposed  method  is  the  better  low  order  controller  but  the 
responses  are  more  oscillatory  than  for  proportional  feedback  control, 
as  was  shown  in  Figure  5.1. 

Figure  5.3  and  the  performance  indices  in  Table  5.3  indicate 
that  the  two  third  order  feedback,  feedforward  control  laws  are  quite 
similar.  All  three  control  laws  in  Figure  5.3  corrected  for  the  +20% 
feed  flow  change  with  no  steady  state  error.  In  Figure  5.4  the  closed- 
loop  responses  of  three  feedback-setpoint  controllers  are  compared  for 
a  +7%  setpoint  change  in  C2.  Again,  the  control  law  designed  from  the 
third  order  model  results  in  a  more  oscillatory  response  than  the 
control  law  designed  by  the  proposed  method.  All  three  responses 
attained  the  desired  steady  state. 

In  general,  the  effectiveness  of  the  reduced  order  controller 

FF  SP  I  FB 

matrices,  KR  ,  KR  and  KR,  depends  upon  the  KR  with  which  they  are 

used.  This  point  was  investigated  for  a  proportional  plus  setpoint 

FB 

control  law  where  KR  was  calculated  from  the  reduced  order  model  and 

SP 

used  with  the  K  matrix  which  was  optimal  for  the  high  order  model. 

For  the  same  +7%  setpoint  in  C2,  this  approach  gave  J  =  37.89,  This 
large  value  was  due  mainly  to  the  resulting  steady  state  C2  value  of 
+11.5%  instead  of  the  desired  value  of  +7%. 


93 


^ — 


G  9 

NIW/Q3  NI  3iVdMD3J 


W  “ 

14  f 
1 

I 

\ 

/ 

\ 

!\ 

f 

-  /  > 

J 

I  . 

1 

■ 

P 

ZH - sz 

- 1 - ^ 

- h 


8 


m 

-  'f 


.. 


*  G  0 

Nlkvai  NI  M03J  WV31S 


51 


-4  o 


a 

j — 

u 

>H 

m 

2 

o 

oo 

h4 

2: 

•V 

PQ 

.H 

in 

o 

00 

jS 

KH 

h~ 

OS 

<3 

41 

m 

in 

o 

O 

00 

Pi 

■ — ^ 

H 

H 

(4 

O 

O 

U 

00 

8 


CD 


w 


n 


ps 

o 

•H 

4J 

O 

3 

T3 


<u 

o 

•H 

4J 

U 

o 

aj 

3 

rH 

a; 

iH 

tH 

4 

O 

a 

4 

rH 

•H 

4-J 

QJ 

4-1 

PS 

T3 

a 

O 

O 

o 

a 

0 

XX  OX  fi 

TJ3ATI  3d  NT  N3NCH 


“4“ 


liif 

* 


-  ft 


'T 

01 


U1 

ui 

i— 

3 


+- 


4 


9X  ET  rjX 

DdH  -NI  NI  13AT! 


c 


CN 

w 

Pi 

B 

M 

Pm 


94 


NIW/B'l  NI  3iVaMDli 


NIW/B~I  NI  MDli  1/W31S 


o 

Pi 

eh 

a 


o 

O  (T> 


O  *H 
Pn  O 
Q  CO 
W 

M  « 

fL,  v£> 


U1 

Lu 


cn 

m 

UO 

w 

Pi 

1 

w 


a 


O 

P 

a 

<U 

o 

U 

•H 

u 

s 

o 

cd 

p 

iH 

'O 

QJ 

r— 1 

rH 

U 

cd 

O 

e 

rH 

•H 

+J  Pi 

(U 

4J 

a  o 

•p 

p- 

O  'H 

o 

O 

O  +J 

e 

z 

i - 1-  - 2 

K - 1 - 

H 

35 

K 

/ 

, 

\  .  i 

1 

\ 

i  { 

/ 

) 

,y/ 

{ 

-  - 1 - -  2 

> 

f-= - 1 - i - 

s 


cn 


"W 


nj 


cu 

■H 


SIT 


&:t. 


ot 


'TO ATI  'Jd  N1  N3N01 


OEM  -NI  NI  13 A3' I 


TIME  IN  MINUTES  TIME  IN  MINUTES 

FIGURE  5.4  SIMULATED  COMPARISON  OF  PROPORTIONAL  FEEDBACK  PLUS  SETPOINT  CONTROL 
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In  process  control  applications,  feedback  control  laws  are 
primarily  used  to  correct  for  unmeasured  disturbances  which  enter  the 
system.  Thus,  it  was  of  interest  to  examine  the  effectiveness  of  a 
controller  designed  by  the  proposed  method  in  such  situations.  Figure 
5.5  shows  the  response  of  the  high  order  evaporator  model  when  upset 
by  a  +20%  feedflow  disturbance  and  controlled  by  the  optimal  tenth 
order  proportional  feedback  controller  or  by  the  third  order  propor- 
tional  feedback  control  law  designed  by  the  proposed  method.  The 
third  order  control  law  results  in  satisfactory  control  but  is,  of 
course,  not  as  good  as  the  optimal  control  law. 

The  above  examples  demonstrate  that  the  reduced  order  control 
laws  designed  by  the  proposed  method  give  excellent  results  when  applied 
to  the  control  of  the  high  order  model.  The  experimental  response  shown 
in  Figure  5.6  indicates  that  the  resulting  control  law  also  performs 
well  when  controlling  the  actual  pilot  plant  evaporator.  The  pro¬ 
portional  control  scheme  results  in  satisfactory  control  of  Wl,  W2,  and 
C2  after  two  20%  step  changes  in  feed  flow.  (The  proportional  feed¬ 
back  control  matrix  for  Figure  5.6  is  shown  in  Table  5.1).  This  exper¬ 
imental  response  represents  a  significant  improvement  over  previously 
reported  results  for  a  conventional  multiloop  controller  [4].  A  com¬ 
parison  of  control  by  the  proposed  method  with  control  by  this  multi¬ 
loop  control  scheme  is  presented  in  Chapter  Six,  for  proportional  plus 
integral  control. 

In  the  development  of  this  new  approach  for  designing  a  reduced 
order  control  law,  it  appeared  that  a  possible  application  might  be  for 
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TIME  IN  MINUTES  TIME  IN  MINUTES 

FIGURE  5.5  SIMULATED  COMPARISON  OF  PROPORTIONAL  FEEDBACK  CONTROL  OF  A  FEED  DISTURBANCE 
(10L(2018)/D,  +20%F/FB/100PT,  3RED10/3040A,  3040B) 
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FIGURE  5.6  EXPERIMENTAL  RESPONSE  OF  PROPORTIONAL  FEEDBACK  CONTROL  OF  A  FEED  DISTURBANCE 
(PP/D,  20%F/FB/ 3RED10/FB19) 


"implied"  control.  Suppose,  for  example,  there  is  a  key  state  variable 
xk->  which  cannot  be  measured,  but  must  be  controlled.  It  appeared  that 
if  the  optimal,  high  order  control  law  for  the  system  was  designed  with 
heavy  weighting  on  x^,  then  the  methods  of  this  chapter  could  be 
applied  to  obtain  a  reduced  order  control  law  which  was  not  a  function 
of  xk,  but  which  would  still  control  it  satisfactorily.  This  strategy 
was  attempted  in  the  evaporator  simulation  studies.  The  variable  to 
be  controlled  was  taken  as  C2  and  several  feedback  control  laws  were 
tried,  all  without  success.  It  is  possible  that,  with  further  study 
into  the  structure,  controllability  and/or  observability  of  the 
system,  a  satisfactory  "implied"  control  scheme  could  be  devised. 
However,  this  was  not  pursued  as  part  of  this  study. 

5.4  CONCLUSIONS 

A  method  has  been  proposed  for  deriving  reduced  order  control 
laws  which  eliminate  selected  state  variables  from  high  order  control 
laws.  In  effect,  an  incomplete  state  feedback  control  law  is  derived 
from  a  state  feedback  control  law  using  a  modal  analysis  of  the  high 
order  state  space  model.  Both  the  high  order  and  low  order  control 
laws  may  contain  feedforward,  integral,  and  setpoint  modes.  The 
feasibility  of  the  proposed  method  has  been  demonstrated  in  simulation 
and  experimental  studies  involving  a  pilot  scale  evaporator.  When 
applied  to  the  evaporator  model,  the  low  order  control  laws  designed 
using  the  proposed  method  gave  better  results  than  controllers 
designed  from  reduced  order  models  which  had  been  obtained  using 
model  reduction  techniques,  as  indicated  by  the  smaller  value  of  the 


performance  index. 
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CHAPTER  SIX 


MODEL  REDUCTION  AND  REDUCED  ORDER  CONTROL 
OF  A  PILOT  PLANT  EVAPORATOR 


ABSTRACT 


The  literature  dealing  with  the  structure  of  the  field  of 
reduced  order  control  law  design  is  reviewed.  The  most  promising 
methods  were  applied  to  the  experimental  control  of  a  pilot  plant 
evaporator.  The  control  law  which  gave  the  best  overall  control, 
of  the  methods  tried,  was  designed  by  the  elimination  of  selected 
state  variables  from  the  feedback  portion  of  the  high  order  control 
law  for  the  process  using  a  modal  analysis.  The  subject  of  model 
reduction  is  briefly  treated.  It  is  shown  that  most  of  the  modal 
methods  of  model  reduction  produce  one  of  two  basic  results  or  two 
modifications  to  one  of  these  basic  results. 
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6.1  INTRODUCTION 

This  chapter  considers,  in  more  detail,  the  general  problems 
of  model  reduction  and  reduced  order  control  law  design  which  have  been 
discussed  for  specific  topics  in  Chapters  Three,  Four,  and  Five. 

The  models  of  some  processes  are  of  such  high  order,  and 
their  high  order  multivariable  control  laws  so  complex,  that  they  are 
impractical  for  use  in  simulation  and  control.  The  simplification  of 
process  models,  and  the  design  of  control  laws  which  require  only  a 
small  subset  of  the  states  of  the  model  have  been  the  topic  of  many 
investigations.  These  simplified  models  and  control  laws  can  then  be 
used  to  study  and  control  the  process.  This  chapter  summarizes  the 
work  done  in  these  two  general  fields  as  well  as  applying  some  of 

the  methods  for  the  model  reduction  and  control  of  a  pilot  plant 
evaporator. 

Figure  6.1  shows  how  a  reduced  order  discrete-time  model,, 
and  a  reduced  order  control  law  can  be  derived  from  a  high  order 
continuous -time  model.  The  desired  reduced  order  model  and  control 
law  are  to  be  in  the  discrete- time  form  for  ease  in  computer  simula¬ 
tion  and  control  system  application  using  an  on-line  computer. 

The  model  reduction  paths  are  numbered  1  and  4  in  Figure 
6.1.  Most  of  the  literature  deals  with  the  continuous- time  model 
reduction  of  path  1  with  work  being  done  in  the  following  areas: 

a)  Reduction  of  a  state  space  model  using  a  modal  analysis  of 
the  system.  This  topic  is  discussed  later  in  this  chapter. 
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b)  Reduction  of  a  state  space  model  by  minimizing  some  measure  of 
difference  between  the  high  order  and  the  low  order  models. 

One  of  these  methods  involves  a  least  squares  analysis  which 
was  discussed  in  detail  in  Chapter  Four. 

c)  Reduction  of  a  transfer  function  model.  This  has  been  discussed 
using  different  approaches  by  Chen  and  Shieh  [4],  Gibilaro  and 
Lees  [20] ,  Sinha  and  Bereznai  [35],  Hsia  [23]  and  Fellows  et  al 
[18]. 

The  reduction  of  discrete- time  models  (path  4)  was  discussed 
in  Chapter  Three  using  a  modal  analysis  and  in  Chapter  Four  using  a 
least  squares  analysis.  One  form  of  the  least  squares  approach  is 
applied  to  the  evaporator  model  later  in  this  chapter. 

The  various  methods  of  calculating  a  reduced  order  control 
law  from  a  high  order  discrete- time  model  shown  in  Figure  6.1  are: 

a)  Model  reduction  followed  by  control  system  design  based  only 
on  the  reduced  order  model  (path  4  followed  by  path  7).  This 
has  been  considered  by  Nicholson  [31  -  33]  and  Andersop  tl] 
for  state  space  models  and  Bereznai  and  Sinha  [2]  and  Chen 
and  Shieh  [3]  for  transfer  function  models  and  is  also  con¬ 
sidered  in  this  chapter. 

b)  Model  reduction  with  the  added  restriction  that  the  optimal 
control  law  for  this  reduced  order  model  must  be  the  "best" 

suboptimal  control  law  for  the  original  high  order  system. 

This  involves  paths  4  and  7  at  the  same  time.  Rogers  and 

Sworder  [34]  utilize  this  approach. 


c)  Generation  of  the  reduced  order  control  system  directly  from 
the  high  order,  open-loop  model,  without  calculating  a  reduced 
order  model.  (Path  6)  This  approach  has  been  considered  by 
Dabke  [7,8],  Levine  and  Athans  [26],  Kosut  [25]  and  Davison 

et  al.  [9,10,15-17]. 

d)  Calculation  of  the  reduced  order  control  law  from  a  high  order 
control  system  designed  for  the  original  high  order  model. 
(Path  5  followed  by  path  8).  Kosut  [25]  requires  the  optimal 
control  law  for  two  simple  forms  of  his  controller.  One 
result  using  this  approach  was  derived  in  Chapter  Five  and  is 
compared  to  other  approaches  later  in  this  chapter. 

The  rest  of  this  chapter  is  organized  as  follows.  Section 
6.2  discusses  the  modal  approach  to  model  reduction  in  more  detail. 

It  also  summarizes  the  model  reduction  and  reduced  order  control  law 
design  techniques  which  are  applied  to  the  pilot  plant  evaporator. 
Section  6.3  discusses  the  model  reduction  results  and  Section  6.4  dis¬ 
cusses  the  application  of  the  reduced  order  control  laws  to  the  con¬ 
trol  of  the  pilot  plant  evaporator. 

6-2  MODEL  reduction  and  control  law  design  methods  to  be  applied  to 

THE  EVAPORATOR 

This  section  discusses  the  methods  of  model  reduction  and 
reduced  order  control  law  design  which  are  applied  to  the  evaporator 
(described  in  Appendix  A)  later  in  this  chapter.  The  general  form  of 
the  model  and  control  law  will  first  be  discussed.  These  are  shown 
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in  Figure  6.1a,  in  which  the  model  and  control  laws  have  been  added  to 
the  basic  structure  of  Figure  6.1.  The  nomenclature  in  Figure  6.1a 
include : 


x  —  n-dimensional  state  vector 

— 1  =  /“dimensional  subset  of  x  which  will  be  retained  in  the  low 
order  model  or  control  law 

2^2  =  (n“/)“ dimensional  subset  of  x 

— r  =  /“dimensional  state  vector  of  the  reduced  order  model 
u  =  (m+q) -dimensional  input  vector  which  can  be  partitioned  into 
an  m-dimensional  control  vector,  u,  and  a  q-dimensional 
disturbance  vector  d,  as 


A 

u 


d 


(6.11) 


X.  ~  p-dimensional  output  vector 

SP  qp 

1  -  a  constant  setpoint  vector  such  that  y"  is  the  desired  value 

of  y£,  i  »  i,  P 

A  A 

=  ’  b  b  £>  &  =  constant  coefficient  matrices  of  appropriate  dimen¬ 
sions  representing  the  process.  The  input  forcing  matrices 
A  A 

B  and  A  can  be  partitioned  to  agree  with  the  partitioning 
,  A 

of  u  as 


2  ~  5)  (6.12) 

t  =  g) 


j  -  counter  for  time  intervals  such  that  x(j)  denotes  x(t)  at 
t  =  jT,  where  T  is  the  discrete  time  interval. 


(6.13) 


. 
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FB  FB  FF  SP  I 

=1  3  =2  3  =  3  =  3  =  =  constant  coefficient  matrices  of  approp¬ 

riate  dimensions  representing  the  parameters  of  the  control 
law 

Ji  ~  integral  contribution  to  the  control  law  represented  for 
continuous- time  systems  as 


*■£ 


/  SP,  J 
(X  -  Z  )  dt 


(6.14) 


or  for  discrete- time  systems  as 

j 

±(j)  =  ^2  (Z(i)  "  ZSP) 

i=0 


(6.15) 


The  subscript,  R,  refers  to  a  reduced  order  vector  or  matrix 
and  the  matrices  BR  and  ^  can  be  partitioned  (as  were  ^  and  A)  into 


A 


-  <2*  2r> 


(6.16) 


4  *  2r) 


(6.17) 


6-2.1  Modal  Approach  to  Model  Reduction 

The  modal  methods  of  model  reduction  are  those  which  retain 
specific  modes  (or  eigenvalues)  of  the  original  high  order  model  in 
the  reduced  order  model.  These  methods  make  use  of  the  modal  analysis 
of  the  high  order  system  as  follows.  A  matrix  M  exists  such  that 


—  — 

X1 

“m.  m„1 

z , 

-1 

=  1  =2 

-1 

2_ 

Jh  \ 

i 

CM 

N| 

_ 1 

(6.18) 


which  will  transform  Equation  (6.1)  into  its  Jordan  canonical  form 
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and  Equation  (6.3)  into  its  Jordan  canonical  form 


(6.19) 


'zj/j  +  1) 

z.2(j  +  !) 


=  ]_  - 


0  a, 


2_ 


%(j) 


Zo(j) 


+ 


1 

fa 


A 

u 


(6.20) 


where  z  = 


n- dimensional  canonical  state  vector  partitioned  into  an 
/-dimensional  vector  z ^  and  an  (n  -  /  )rdimensional  vector 


~2 


J  = 


ii 


0 


a  = 


0  j2 


=  1  2 


2  22j 


=  V  A  M 


=  V  0  M 


(6.21) 


(6.22) 


^  = 


A 

Si 

A 

Sa 


*1 

% 


=  V  A 


(6.23) 


(6.24) 


V  = 


2l  =2 


V  V 
=3  =4 


=  M 


-1 


A  A 


(6.25) 


Matrices  G  and  8  can  be  partitioned  (as  were  B  and  ^)  into 


. 
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2  -  (£>  M)  (6.26) 


Even  though  there  have  been  many  papers  written  on  the  modal 
approach  to  model  reduction,  most  of  the  methods  produce  one  of  two 
basic  results  or  two  modified  forms  of  one  of  these  basic  results. 

Most  of  the  modal  methods  calculate  the  same  AR  or  <*>R.  They  differ 
only  in  the  approach  taken  and  in  the  resulting  BR  and  AR.  Most  of 
the  comments  to  follow  will  refer  specifically  to  the  continuous¬ 
time  reduction  methods  used  to  obtain  Equation  (6.6)  from  Equation 
(6.1). 


One  of  the  basic  modal  approaches  was  presented  in  different 
forms  by  Marshall  [27]  and  Chidambara  [5]  (Chidambara * s  method  marked 
C2).  The  equivalence  of  the  two  formulations  was  noted  by  Graham  [21]. 
They  assume  that  the  modes  which  do  not  contribute  significantly  to 
the  responses  of  the  states  of  interest,  reach  their  final  steady 
state  values  immediately  after  a  disturbance  enters  the  system.  The 
mam  advantage  >of  this  result  is  that  it  provides  a  reduced  order 
model  with  exact  steady  state  agreement  with  the  original  high  order 
model  for  step  disturbances.  Davison  [5,11],  however,  objects  to 
this  model  because  the  retained  modes  are  not  excited  in  the  same 
proportions  in  the  reduced  order  model  as  they  are  in  the  high  order 
model.  The  equivalence  of  these  results  is  shown  in  Appendix  C. 

The  second  basic  modal  approach  was  presented  in  slightly 
different  forms  by  Nicholson  [33]  and  Davison  [l3].  The  equivalence 
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of  their  formulations  is  shown  in  Appendix  C.  Their  resulting  reduced 
order  model  has  its  modes  excited  in  the  same  proportions  as  they  were 
in  the  high  order  model.  The  main  disadvantage  of  this  result  is  that 
a  steady  state  error  exists  between  the  high  order  and  the  low  order 
responses  for  sustained  disturbances.  Davison  .[11]  claims  that  this 
error  will  be  small  if  the  order  of  reduction  is  not  large. 

Several  different  formulations  have  been  presented  in  the 
literature  for  modifications  to  the  result  of  Davison  and  Nicholson 
to  provide  steady  state  agreement.  These,  however,  all  reduce  to  only 
two  different  results. 

Davison  [12]  presented  one  of  these  revisions;  each  state 
variable  of  the  reduced  order  model  is  adjusted  by  a  ratio  of  the 
desired  steady  state  to  the  steady  state  obtained  for  the  unmodified 
reduced  order  model.  For  systems  with  a  single  input,  Davison! s 
modified  result  can  be  arranged  into  the  standard  form  of  Equation 
(6.6).  His  modification  alters  both  AR  and  BR  but  retains  the  eigen¬ 
values.  However,  if  there  is  more  than  one  input,  this  arrangement 
is  not  possible.  Also,  his  modification  cannot  be  applied  to  models 
with  a  singular  A  or  AR  (unmodified  form  using  Davison's  basic  method 
mentioned  above)  since  A  ^  and  AR  ^  are  required.  However,  in  Davison's 
example,  both  A  and  AR  are  singular.  He  gets  around  this  problem  by 
changing  the  parameters  in  his  original  model  slightly,  so  that  A 
is  no  longer  singular  [14]. 

The  second  revision  to  Davison's  basic  method  has  been  pre¬ 
sented  by  Davison  [5,  11]  and  by  Fossard  [19].  This  result  was  also 
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obtained  by  Chidambara  [6]  (Chidambara ’ s  method  marked  Cl)  and  Graham 
[21],  not  as  revisions  to  Davison’s  basic  method,  but  in  independent 
derivations.  The  equivalence  of  these  results  is  derived  and  discussed 
in  Appendix  C.  This  revision  can  be  written  in  two  equivalent  forms 
as: 


or  as 


*Jl  =  =R  +  2r  (6.28) 

2£r  =  2^  +  Er  u  (6.29) 

\  =  ^r2r  +  SrH  +  Ir$  (6.30) 


where 


B° 


AM 

=R 


=  a  general  ^-dimensional  vector 
A 

=  B  using  Davison  s  basic  method 

— K 
A 

=  B  using  Marshall  s  basic  method 
—  K 


It  can  be  seen  from  Equation  (6.30)  that  this  result  reduces  to  Mar- 
shall’ s  result  for  the  case  of  u  =  0. 


Even  though  all  four  of  these  results  are  equivalent,  the 
form  presented  by  Davison  [5,11]  cannot  be  used  when  either  A  or 
A  is  singular. 


Graham  [21]  discusses  the  idea  of  applying  this  method 

with  a  non-zero  initial  vector  for  x  in  Equation  (6.30).  His  initial 

K 

value  vector  is 


xR(0+)=  Mx  xx(0+)+  V2  x2(Q+)}+  =r  H(°+)  (6-31> 


This  modification  provides  better  agreement  with  the  high  order 
model  earlier  in  the  response  but  with  the  disadvantage  that  there  is 
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a  large  error  at  the  start  of  the  trajectory.  If  the  model  were  to 

be  used  for  control  purposes,  the  initial  time  error  could  not  be 

accepted.  In  the  formulation  of  Equations  (6.28)  and  (6.29)  the 

vector  x*  ,  at  time  zero,  can  be  set  arbitrarily.  Thus,  in  the 

examples  in  Section  6.3.1,  to  ensure  that  x  (0)  =  x  (0)  the 

R  — 1  ’ 

initial  value  of  x^  will  be  specified,  using  Equation  (6.29),  to 
be 


A . 


=  2^(0)  -  Er  u(0) 


(6.32) 


The  modal  methods  which  are  applied  to  the  evaporator  model 
in  this  chapter  are  those  presented  by  Marshall  [27],  Davison  [13]  and 
Fossard  [19]  for  the  reduction  of  continuous-time  models  and  extended 
for  use  with  discrete-time  models  in  Chapter  Three.  Each  of  these 
methods  produce  the  same  A^  and  0^  as 


-1  ■  •  -  -1 . .  -1 


h  m  Ml  h  Si  -  -  a2  v4  v3  =  Ax  +  a2  m3  m3 


(6.33) 


and 


-1 


2r  =  Si  “i  «!  "  -  ±1  -  t2  lh  1  V3  =  +  g2  M3  ^  1  (6.34) 


-1 


These  equivalences  can  be  shown  using  relations  in  Appendix  D. 
Expressions  for  BR,  ER  and  for  each  of  these  methods  are  presented 
in  Table  6.1. 


In  the  application  of  these  methods,  x ^(Q)  is  calculated 
using  Equation  (6.32).  Also,  no  reference  will  be  made  to  indicate 


whether  the  continuous— time  model  or  the  discrete— time  model  was 
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identical  and  are  shown  by  Equations  (6.33)  and  (6.34). 


reduced,  since,  as  shown  in  Chapter  Three,  the  same  discrete-time  model 
is  obtained. 


6.2.2  Least  Squares  Approach  to  Model  Reduction 

The  application  of  least  squares  to  model  reduction  falls 
within  the  general  category  of  methods  which  minimize  a  function  of 
the  difference  between  the  response  of  the  high  order  model  and  that 
of  the  low  order  model.  Specifically,  least  squares  minimizes  the  sum 
of  the  squares  of  the  difference  between  the  high  order  model  response 
and  resulting  reduced  order  model  response  for  all  the  states  in  the 
reduced  order  model. 


In  this  approach  the  transition  matrix,  d>  is  calculated 

— R 

using  least  squares  while  is  calculated  to  give  agreement  between 
the  high  order  and  reduced  order  models  after  the  initial  transient 
dies  out.  The  high  order  model  and  a  sequence  of  uniformly  distrib¬ 
uted  random  numbers  are  used  to  calculate  the  data  for  use  with 
least  squares,  as  discussed  in  Chapter  Four.  The  results  are  sum¬ 
marized  here  for  the  reduction  of  discrete-time  models  with  at  least 
one  integrating  state  (that  is,  an  eigenvalue  of  0  equal  to  one). 


If  the  reduced  order  state  vector  x  is  partitioned  as 

R 


(6.35) 


corresponding  to  the  integrating,  Xj,  and  the  non-integrating,  x^, 
state  variables  of  x^5  then  0^  can  similarly  be  partitioned  as. 


■ 
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The  partition  of  d>R  corresponding  to  the  integrating  state  variables 
can  be  pre-specified  as 


(6.37) 


The  dimension  of  the  identity  matrix  in  Equation  (6.37)  is  the  same 

as  the  number  of  integrating  state  variables.  The  partition  of  ± 

— R 

corresponding  to  the  non-integrating  states  is  calculated  using 
least  squares  as 


"  (ST  1  f  W  (6.38) 

where 


wT  (1) 

z  = 

T 

4  <2> 

and  W  = 

wT  (2) 

1 

' — ✓ 

■  tyf 

_ I 

- 1 

H 

•  •  *15^ 

and  where  k  is  the  number  of  data  vectors.  The  jth  row  of  W  is 
defined  as 


-(j)  =  =1N  Vj)  +  =2  -2(j) 


where  has  been  partitioned  as 


(6.40) 


(6.41) 


corresponding  to  the  partitioning  of  xR  in  Equation  (6.35).  The 
elements  of  vectors  x^  and  x^  are  specified  as  random  numbers. 


After  has  been  obtained,  AR  is  calculated  to  ensure  that 
the  reduced  order  and  high  order  models  give  the  same  response  for 
large,  values  of  the  counter,  j.  Thus, 
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£ 


-R  (j  +  D  "  £R  2£r  (i) 


R  -R 


A 

u, 

1 


(6.42) 


where  ^ 
R 


A 

u. 

1 


1 
— R 


.th  n  A 

l  column  of  A 
— R 

.  th  _  A 

l  element  of  u 


for  large  values  of  the  counter  j  for  a  step  change  in 
u±.  (u^  =  0,  k  ±  i) 


The  values  of  the  elements  of  x^  (j  +  1)  and  x^  (j)  corresponding  to 
the  non-integrating  states  are  equivalent  and  are  the  steady  state 
values  of  these  states  for  the  non-zero  value  of  u. ,  (u^  =  0,  k  4  i) , 
as  calculated  from  the  high  order  model.  The  values  of  the  elements 
of  (j  4-  1)  corresponding  to  the  integrating  states  are  calculated 


as 


+  =  tx  +  t2  ±2(1)  +  £i  “i 


(6.43) 


where  =  i^  column  of  ZV.  ,  and  the  elements  of  x^  ( j )  corresponding 

to  the  integrating  states  are  specified  by  random  numbers  and  the 

elements  of  x^(j)  corresponding  to  the  non-integrating  states  are 

the  above  mentioned  steady  state  values  for  the  same  non-zero  value 
A 

of  u.  . 
l 


Equation  (6.42)  provides  for  steady  state  agreement  of  the 
non- integrating  states  and  agreement  of  the  rate  of  change  of  the 
integrating  states  after  the  transients  die  out. 

6.2.3  Approaches  Used  for  Reduced  Order  Control  Law  Design 

The  modal  analysis  and  the  least  squares  analysis  (using 


random  data)  were  used  to  design  reduced  order  control  laws  using  the 
following  paths  shown  in  Figure  6.1: 

a)  Path  4  followed  by  path  7  -  the  calculation  of  a  reduced 
order  model  and  the  design  of  the  reduced  order  control  law 
directly  from  this  model  using  standard  techniques. 

b)  Path  5  followed  by  path  8  -  the  calculation  of  the  high  order 
optimal  control  law  and  the  reduction  of  this  controller  to 
form  the  reduced  order  control  law. 

In  the  first  approach,  the  reduced  order  control  law  is 

calculated  from  a  reduced  order  model.  The  modal  approach  to  model 

reduction  used  in  this  step  was  that  of  Marshall  as  described  in 

Chapter  Three  and  by  Marshall  [27]  and  was  chosen  from  the  modal 

methods  discussed  in  Section  6.2.1,  since  it  was  the  only  method 

whose  result  could  be  arranged  in  the  standard  form  of  Equation  (6.8) 

and  which  guaranteed  steady  state  agreement  between  the  high  order 

and  the  reduced  order  models.  The  least  squares  approach  used  in 

this  step  is  that  outlined  in  Section  6.2.2.  It  calculates  <t> 

— R 

using  random  data  and  for  agreement  at  large  times.  Once  the 
required  reduced  order  model  is  obtained,  the  desired  reduced  order 
control  law  is  calculated  using  any  applicable  method.  In  this  work 
the  optimal  control  approach  outlined  by  Newell  [28]  was  used. 

In  the  second  approach,  the  reduced  order  control  law  is 
calculated  from  the  high  order  control  law  for  the  system.  A  modal 
analysis  and  a  least  squares  analysis  (with  random  data)  are  used 
to  eliminate  selected  states  from  the  high  order  control  law.  The 
modal  approach  is  that  which  is  discussed  in  detail  in  Chapter  Five 
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and  so  it  will  not  be  discussed  further  here.  The  least  squares  ap¬ 
proach  is  similar  to  that  applied  to  model  reduction  in  Chapter  Four. 
The  reduced  order  control  law  is  calculated  by  fitting  the  response  of 
the  reduced  order  controller  to  the  response  of  the  high  order  control¬ 
ler.  Appendix  E  shows  that  the  same  reduced  order  controllers  can  be 
obtained  from  two  different  (but  equivalent)  formulations,  as: 

a)  Calculating  the  control  matrices  so  that  the  control  vector, 

_u,  as  calculated  with  the  reduced  order  control  law,  agrees 
with  that  calculated  by  the  high  order  control  law. 

b)  Calculating  the  control  matrices  so  that  the  closed-loop 
response  of  the  high  order  system,  when  controlled  by  the 
reduced  order  control  law,  agrees  with  that  obtained  when 
controlled  by  the  high  order  control  law. 

The  resulting  reduced  order  control  matrices  are  calculated  from 

T  T  -1  T 

Kr  =  (Z  Z)  z1  W  (6.44) 

where  KR  =  the  control  matrices  to  be  calculated.  In  the  case  of 
proportional  feedback,  feedforward,  setpoint,  integral 
feedback  control. 


v_  =  rKFB  KFF  KSP  TC h 
=R  -R  ’  =R  ’  5  ^R 


T 

z  (1) 

T 

w  (1) 

z  = 

T 

z  (2) 

and  W  = 

T 

w  (2) 

• 

• 

• 

T 

_z  (k) 

• 

• 

wT  (k) 

(6.45) 


(6.46) 


. 

' 
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where  k  is  the  number  of  data  vectors.  The  vectors  _z(j )  which  make 
up  the  rows  of  Z,  for  the  case  of  proportional  feedback,  feedforward, 

setpoint,  integral  feedback  control,  represent  the  vectors  x  (j  )  , 

SP  P 

Z  and  s_(j  )  so  that 


z(j)  = 


£1 

d  (j) 
SP 

z 

s  (j) 


(6.47) 


The  vectors  w(j)  which  make  up  the  rows  of  W,  for  the  case  of  pro¬ 
portional  feedback,  feedforward,  setpoint,  and  integral  feedback  con 
trol,  are  calculated  as 


r  .  \  ...FB  ,  ,  . 

w(j)  =  (j  ) 


4-  VFB 

+  S2  —2 


FTi1 

(j)+K  d(j) 


+  KSP  ZSP  +  K1  s  (j  ) 


(6.48) 

The  values  of  ( j )  ,  Xg  (j  )  ,  d_(j  )  ,  Z  ^  and  s_(j  )  in  Equations  (6.47) 

and  (6.48)  are  specified  as  a  sequence  of  uniformly  distributed  ran¬ 
dom  numbers . 


In  the  application  of  the  result ,  if  a  control  law  which 
does  not  include  all  four  control  modes  discussed  in  the  above  analysis 
is  desired,  then  the  terms  in  Equations  (6.45),  (6.47)  and  (6.48) 
corresponding  to  the  undesired  control  matrices  are  simply  left  out 
of  the  analysis.  For  the  example  of  proportional  feedback  control 
only j  Equations  (6.45),  (6.47)  and  (6.48)  will  appear  as  1 


..  JO  -> 
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z(j)  =  ^(3) 

w(j)  =  k[B  Xx(j)  +  K™  x2(j) 
6.3  REDUCTION  OF  AN  EVAPORATOR  MODEL 


(6.49) 


The  model  reduction  methods  discussed  in  Sections  6.2.1 
and  6.2.2  have  been  applied  to  the  reduction  of  the  tenth  order  evap¬ 
orator  model  described  in  Appendix  A.  These  results  are  discussed 
in  this  section  .  The  various  reduced  order  models  calculated  are 
shown  in  Appendix  F.  Most  of  the  reduced  order  models  discussed  in 

this  section  are  third  order  models  in  which  the  reduced  order  state 

T 

vector  x.£  =  [Wl,  W2,  C2]  .  These  are  the  three  process  variables 
which  are  of  most  interest  since,  as  is  discussed  in  Section  6.4, 
they  are  the  process  variables  which  must  be  controlled. 


6.3.1  Modal  Methods 

The  modal  methods  applied  to  the  reduction  of  the  evapor¬ 
ator  model  are  those  presented  in  Table  6.1. 


Figure  6.2  compares  the  time  domain  response  to  a  +20%  step 
change  in  feedflow  using  various  third  order  open  loop  models  and  the 
tenth  order  open  loop  model.  It  shows  that  the  third  order  model 
calculated  using  Davison's  method  has  considerable  steady  state 
error.  (It  should  be  pointed  out,  however,  that  the  fifth  order 
model  calculated  using  Davison's  method  has  negligible  steady  state 
error.)  Figure  6.2  also  shows  that  the  third  order  models  calculated 
using  Marshall's  method  and  Fossard's  method  are  identical,  lead  the 
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tenth  order  model  slightly,  and  result  in  zero  steady  state  error. 

The  models  calculated  by  Marshall's  method  and  Fossard's  method  yield 
the  same  time  domain  response  since  the  input  is  constant  (+20%  step 
in  feedflow,  open  loop).  In  order  to  have  zero  initial  error  using 
Fossard  s  method,  the  initial  values  of  the  elements  of  x  were  set, 
using  Equation  (6.32)  to 


"Wl  (0)  " 

~  0.0014' 

Xj/  (°)  =  (0)  -  Er  u(0)  = 

W2(0) 

= 

0.0024 

_C2  (0)  _ 

--0.097  _ 

Figure  6.3  compares  the  response  of  the  third  order  model 
calculated  using  Fossard's  method  when  x^(0)  =  2£^(0)  and  when  xR(0 ) 
is  defined  by  Equation  (6.31).  The  response  with  x  (0)  x  (0) 
converges  to  the  high  order  model  response  faster  than  the  other 
response.  However,  it  has  a  large  initial  error  (at  zero  time) 
in  C2 ,  which  is  undesirable.  Also,  for  this  example,  the  response 
with  x.R(0)  =  xx(0)  is  satisfactory  since  the  error  is  not  large 
anywhere  along  the  response. 

Figure  6.4  demonstrates  that  the  agreement  with  the  tenth 

order  model  becomes  better  as  the  order  of  the  reduced  order  model 

is  increased,  as  would  be  expected.  The  response  of  the  fifth  order 
T 

model,  with  _xR  =  [Wl,  Cl,  HI,  W2,  C2],  agrees  better  with  the  response 
of  the  tenth  order  model  than  does  that  of  the  third  order  model. 

Figure  6.5  demonstrates  that  the  eigenvalues  to  be  retained 

in  the  reduced  order  model  must  be  carefully  selected  by  showing  the 

T 

response  of  two  fourth  order  models  with  *  =  [Wl,  HI,  W2 ,  C21.  When 

— R 
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TIME  IN  MINJTES  TIME  IN  MINLJTES 

FIGURE  6.4  COMPARISON  OF  DIFFERENT  ORDER  REDUCED  ORDER  MODELS  CALCULATED  BY  MARSHALL'S  APPROACH 
(10L(2018),  5L(2020),  3L(2021)/D,  +20%F/OL/2018,  2020,  2021) 
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the  four  eigenvalues  retained  in  0R  were  those  of  0  which  were  closest 
to  one,  poor  agreement  in  the  C2  response  resulted.  The  modal  matrix 
for  the  tenth  order  model  in  Appendix  A,  shows  that  none  of  these  four 
largest  eigenvalues  affect  HI  (x(5))  significantly.  However,  the 
fifth  largest  eigenvalue  does  affect  HI  significantly.  Thus,  when 
the  fifth  largest  eigenvalue  (0.7354)  replaces  the  fourth  largest 
eigenvalue  (0.9212),  much  better  results  were  obtained. 

6.3.2  Least  Squares  Method 

Fi-Shi’e  6.6  shows  the  comparison  of  the  least  squares  approach 
with  the  modal  approach  of  Marshall.  The  response  of  the  third  order 
model  obtained  by  the  least  squares  approach  agrees  slightly  better 
with  that  of  the  high  order  model  than  does  the  response  of  the  third 
order  model  obtained  by  the  modal  analysis.  Table  6.2  shows  that  the 
third  eigenvalue  of  0R  from  the  least  squares  model  is  larger  than  the 
corresponding  eigenvalue  of  0R  from  the  modal  model.  This  larger 
eigenvalue  causes  the  slower  response  and  the  better  agreement. 

6.4  REDUCED  ORDER  CONTROL  OF  A  PILOT  PLANT  EVAPORATOR 

The  reduced  order  control  laws  discussed  in  Section  6.2.3 
were  applied  to  the  pilot  plant  evaporator  described  in  Appendix  A. 

The  objective  was  the  lowest  order  control  law  which  would  still 
satisfactorily  stabilize  the  system,  The  result  was  a  third  order 
control  law  with  vector  3^,  containing  Wl,  W2  and  C2.  The  product 
concentration,  C2,  must  be  included  in  x_  since  regulation  of  C2  is 
the  primary  control  objective.  Wl  and  W2  must  be  included  since  these 
are  integrating  states  (corresponding  to  eigenvalues  of  0  equal  to 
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TABLE  6.2 


EIGENVALUES  OF  <£_  FOR  TWO  THIRD  ORDER 

■*R 


Eigenvalue 

Number 


Marshall '  s 
Modal  Analysis 


1 


1. 


2 


1. 


0.9600 


MODELS 


Least  Squares 
Analysis 

1. 

1. 

0.9662 
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one)  which  would  tend  to  exceed  physical  operating  limits  if  not  con¬ 
trolled.  The  remaining  state  variables  are  less  critical  and  their 
values  are  not  of  direct  concern. 

In  the  discussion  to  follow,  the  control  laws  will  be 
referred  to  by  an  abbreviated  code  for  clarity  and  brevity.  The 
third  order  control  laws  which  were  designed  by  reducing  the  order 
of  the  tenth  order  optimal  control  law,  are  referred  to  as  "3RED10" 
while  those  which  were  optimal  for  a  third  order  model  are  referred 
to  as  "30PT".  Following  this  designation  is  reference  to  the 
approach  used  in  the  model  or  the  control  law  reduction  as  either 
(M)  if  a  modal  analysis  is  used  or  "(LS)"  if  a  least  squares  analysis 
is  used.  Thus,  the  third  order  control  law  reduced  from  a  tenth 
order  controller  using  the  modal  analysis  is  referred  to  as  "3RED10(M)". 

The  "optimal"  control  laws  mentioned  above  were  calculated 
using  dynamic  programming  as  described  by  Newell  [28].  The  weighting 
rostrices  in  the  performance  index  used  with  dynamic  programming  are 
shown  in  Appendix  A.  It  should  be  pointed  out  here  that  different 
weighting  matrices  could  have  been  used  for  the  third  order  control 
laws  which  would  have  affected  the  results.  However,  to  be  able  to 
compare  the  results  it  was  felt  that  each  state  should  be  weighted 
the  same  amount  in  all  calculations. 

The  control  laws  were  applied  experimentally  for  the  con¬ 
trol  of  the  pilot  plant  evaporator  using  the  general  scheme  shown  in 
Figure  6.7.  The  pilot  plant  is  interfaced  to  the  IBM  1800  computer 
through  the  DDC  (direct  digital  control)  monitor  programs.  The 
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actual  multivariable  control  programs,  written  by  Newell  [29],  obtain 
the  process  measurements  from  the  DDC  loop  record  table,  calculate 
a  set  of  control  values,  and  then  change  the  setpoints  for  these 
control  loops  in  the  DDC  loop  record  table.  The  state  estimation 
step,  which  has  been  discussed  in  detail  by  Hamilton  [22],  is  not 
required  in  this  application  since  all  the  control  laws  used  are  a 
function  of  measurable  state  variables  only.  The  units  conversion 
steps  are  needed  since  the  control  laws  require  measurements  as  nor¬ 
malized  perturbation  variables  and  calculate  a  control  vector  in  this 
same  form.  A  more  complete  description  of  the  control  programs  is 
presented  by  Newell  and  Fisher  [29]  and  of  the  application  of  these 
programs  to  the  control  of  the  evaporator  by  Newell  [28]. 

For  comparison  with  the  multivariable  control  schemes  of 
this  chapter,  Figure  6.8  shows  the  response  of  the  evaporator  when 
upset  by  a  +20%  increase  in  feedflow  and  controlled  by  a  multiloop, 
proportional  plus  integral  control  scheme.  Each  liquid  hold  up  is 
controlled  by  the  outlet  flow  from  the  particular  effect,  while  the 
product  concentration  is  controlled  by  the  steam.  This  was  shown 
by  Newell  [28]  to  be  the  best  multiloop  configuration.  The  control¬ 
ler  constants  and  sample  time  used  in  this  run  are  shown  in  Table  6.3 
and  are  identical  to  the  values  used  by  Jacobson  [24]  for  his  "stand¬ 
ard"  case.  The  sample  time  of  64  seconds  is  the  same  as  that  used 
for  the  multivariable  controllers,  which  will  be  discussed  next.  The 
scales  in  Figure  6.8  are  the  same  as  those  in  the  figures  to  follow 
except  that  the  total  time  is  different  and  the  scale  of  the  levels 
is  shifted  vertically  by  two  inches. 
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TABLE  6.3 


CONTROL  LOOP  SETTINGS  FOR  MULTILOOP  RUNS 


Control 

Interval  Proportional  Integral  Constant 

Loop  (sec.)  Constant  (sec.) 


Bl 

W1 

64 

0.  75 

2048 

B2 

W2 

64 

2.5 

4096 

S 

C2 

64 

1.5 

2728 
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Various  multivariable  control  schemes  have  been  applied  to 
the  control  of  the  evaporator.  These  are  discussed  in  the  following 
sections.  For  the  proportional  feedback  and  the  proportional  plus 
integral  feedback,  control  schemes,  the  results  of  using  a  fifth 
order  control  law  is  also  presented.  This  control  law  is  optimal 
for  the  fifth  order  model  which  was  derived  by  Newell  [28].  This  will 
provide  a  comparison  between  the  work  presented  in  this  chapter  and 
the  previous  work  by  Newell  [28]. 

6-4-1  Proportional  Feedback  Plus  Feedforward  Control 

All  four  third  order  proportional  feedback  plus  feedforward 
(FB  +  FF)  control  laws  which  were  applied  to  the  evaporator,  con¬ 
trolled  it  well  with  only  slight  differences.  The  control  matrices 
used  are  shown  in  Table  6.4,  along  with  the  setpoint  control  matrices 
to  be  considered  later.  The  response  of  the  evaporator  when  con¬ 
trolled  by  the  FB  +  FF  control  law  designed  as 

a)  3RED10(M)  is  shown  in  Figure  (6.9) 

b)  3RED10(LS)  is  shown  in  Figure  (6.10) 

c)  30PT(M)  is  shown  in  Figure  (6.11) 

d)  30PT(LS)  is  shown  in  Figure  (6.12). 

These  responses  are  all  very  similar.  However,  the  con¬ 
trollers  designed  as  3RED10(M)  (Figure  6.9)  and  30PT(LS)  (Figure  6.12) 
provide  slightly  steadier  control  than  the  other  two,  with  less  oscil¬ 
latory  responses,  particularly  in  the  holdups,  Wl  and  W2,  and  in  the 
controls,  Bl,  B2  and  S. 

The  controller  calculated  as  3RED10(LS)  (Figure  6.10)  results 
in  a  larger  offset  in  the  three  controlled  variables,  Wl,  W2  and  C2. 


CONTROL  MATRICES  FOR  PROPORTIONAL  FEEDBACK 


135 


P-i 

co  Pi 
Wll 


Pm  Pi 
Mil 


S) 


TJ 
O 
•H  X 
CO  4-) 


mo  o 

1  1  00 

1 

^  O  CN 

o  o  cn 

rx  x®  in 

on  o  rx 

cm  h  in 

CN  CN  X® 

v®  mo  cn 

rH  xf  00 

in  On  O 

N  CN  cn 

CM  sf  O 

x®  in  o 

1  rH 

1 

rH  i — 1  rH 

1  1 

rH  1  rH 

1 

«H  |  rH 

1 

CN 

rx 

ix 

h  o  n 

no  in  co 

rH  ON  MO 

mo  x®  n 

CN  X®  mO 

cn  v®  cn 

O  ON  00 

n  n  n 

o  n  n 

o  cn  rx 

NT  LO  MO 

CM  NO  NO 

OH® 

1  1 

O  >H  ON 

l  i 

O  i — 1 

1 

O  rH  ON 

1 

00  ON  00 

H  ®  H 

Xf  XT  sf 

X®  MO  MO 

MO  00  CN 

rH  MO  Cn 

o  oo  on 

H  M  N 

X®  CN  i — 1 

Mn  n 

ON  N  O 

ON  f  o 

•  •  • 

•  •  • 

•  •  • 

•  •  • 

cn 

l  I  I 

00  CN  CN 

|  |  | 

lo  <r 

1  1  1 

<r  in 

l _ 1  L 

1 1 _ 

ii 

i 

1  1  1 

1 1 - 

7  *<f 

1 

cn  o 

rx  LO 

m  cn 

rx  ON 

00  CN  X® 

ON  i — 1  N 

co  m  cn 

cn  in  cn 

cn  h  h 

rH  ON  ON 

CN  O  H 

CN  O  CN 

H  O  O 

ON  CN  O 

IH  O  O 

rH  O  O 

<t  O  O 

NT  O  O 

x®  o  o 

x®  O  O 

O  O  O 

1  1  1 

o  o  o 

1 

o  o  o 

1  1  1 

o  o  o 

1  1  1 

ON  N  N 

cn  ix  v® 

cn  h  n 

rH  CN  CN 

cn  N  IN 

cn  oo  cn 

oo  cn  cn 

NO  CN  rH 

lO  H  CO 

x®  ix  <xj 

n  cn  ® 

n  CN  ON 

m  cm  ® 

CN  O  CN 

in  cn  on 

n  cn  on 

o  o  o 

1 

o  o  o 

1 

o  o  o 

i 

o  o  o 

1 

m  oo 

LT)  CN 

cn  cn 

oo  cn 

00  H  N 

'j  m  cm 

oo  cn  oo 

cn  «n  x® 

cn  oo  on 

i — l  cn 

cn  oo  on 

X®  00  On 

CN  ON  ON 

cn  >3-  oo 

CN  On  ON 

CN  ON  On 

L  H°0,. 

cn  o  o 

..HO° 

,  .  rH  °  O 

1 

in  m  oo  r- 

1  1  m  n 

-*1  cn  cn 

l 

O  00  CN 

CN  N  on 

cn  cn  in 

cn  co  fs 

oo  in 

cn  MO  rH 

on  x®  cn 

MO  rH  00 

xt  O0  ON 

m  H  cm 

H  >®  On 

cn  in  on 

1 

rH  1 — I  rH 

1 

rH 

1 

1 — 1 

1 

cn 

MO 

cn 

cn 

n 

H  O  O 

o  <®  o 

h  o  in 

mo  x®  n 

cn  <r  mo 

ON  CN  CN 

o  o  00 

n  n  n 

o  cn  rx 

CN  CN  00 

Xf  NO  NO 

CN  MO  xO 

OH® 

1 

O  rH  ON 

1  1 

O  rH  ON 

1  1 

O  H  On 

1  1 

N  00  00 

O  MO  00 

xr  xr  cn 

X®  MO  MO 

MO  00  CN 

rH  o  m 

O  00  ON 

r — 1  00  rx 

X®  CN  i— 1 

IT)  MO  00 

ON  N  O 

ON  H  O 

CM  xfst 

00  CN  CN 

, ,  x®  n  x® 

,,x®  m  xt 

1 - 1  L 

J  i 

| 

CO 

s 

P 

/—N 

CO 

o 

O 

S 

H 

rH 

1 - 1 

Nw' 

Q 

P 

H 

H 

9 

0 

PM 

PM 

P5 

PM 

O 

O 

cn 

cn 

cn 

cn 

r~ 

1 

1 - 

1 

*<r 

n 

rx 

*sf 

MO 

O 

CN 

MO 

rH 

MO 

CN 

• 

• 

• 

• 

• 

'  • 

n 

1 — 1 

CN 

n 

rH 

CN 

iH 

1 

1 — 1 

1 - 1 

» — 1 

1 

rH 

rH 

IX 

X® 

CN 

O 

cn 

in 

CN 

On 

'sf 

n 

CN 

n 

00 

o 

CN 

CN 

00 

o 

CN 

00 

• 

• 

• 

• 

• 

• 

o 

1 - 1 

1 

ON 

o 

1 

1 - 1 

1 

ON 

On 

rx 

ON 

cn 

cn 

rx 

00 

00 

ON 

ON 

<f 

n 

n 

IX 

ON 

MO 

• 

• 

• 

• 

• 

• 

00 

CN 

CN 

00 

CN 

CN 

1 

1 _ 

| 

cn 

H 

•H 

CO 

Pi 

o 

CO 

•H 

H 

4-1 

Xr 

3 

O 

O 

r — 1 

CM 

Q 

O 

W 

3 

Pi 

PM 

cn 

• 

Sm 

Pm 

3 

O 

Pm  pi 
Mil 

O 

<4H 

H 

X 

CD 

PM 

a) 

X 

+J 

CO  Pi 

WII 

•H 

i — 1 

PM 

X 

X 

CM 

4-1 

4-1 

CO 

•H 

•H 

£ 

£ 

CD 

X 

t>N 

T) 

r— 1 

CD 

X 

Pi 

cn 

1 — 1 

O 

3 

3 

X) 

cn 

O 

£  cn 

0) 

•H 

•H 

cn 

X 

3 

p: 

o 

o 

•frl  t — 1 

cn 

•H 

X  3 

•H 

X 

£  o 

PQ 

4-1  X 

Pm  pim 

cO  O 

P4ll  Pm  pi 

X  3 
4-4  X 

cn 

X 

•H 

CD 

X  QJ 

X 

£ 

3  QJ 

H 

CO  H 

■K 

. 


136 


NIW/B"!  NI  liVdMOld 


j  no  ni  luruvaiui^i 

OLE  oax 

|  - - - +; - vs - 1 — 


O 


XT  0T-  B 

TOAT1  3d  NI  N3ND3 


CEH  -NI  NI  13 All 


TEC  IN  MINUTES  TIME  IN  MINUTES 

FIGURE  6.9  EXPERIMENTAL  PROPORTIONAL  FEEDBACK,  FEEDFORWARD  CONTROL  (3RED10(M)) 

(PP/D,  20%F/FB  +  FF/ 3RED10 (M) /FF21) 
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TIME  IN  MINUTES  TIME  IN  MINUTES 

FIGURE  6.12  EXPERIMENTAL  PROPORTIONAL  FEEDBACK,  FEEDFORWARD  CONTROL  (30PT(LS)) 

(PP/D,  20%F/FB  +  FF/30PT(LS)/FF25) 
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With  perfect  feedforward  control,  there  should  be  no  steady  state  off¬ 
set.  This  offset  results  since  the  control  law  is  calculated  with  a 
least  squares  analysis,  with  no  adjustment  for  steady  states  agree¬ 
ment  between  the  system  controlled  optimally  and  suboptimally .  This 

FB 

offset  could  be  removed  by  designing  alone  using  least  squares, 

FF 

and  then  designing  KR  for  zero  offsets  in  the  three  controlled  var¬ 
iables,  as  was  discussed  by  Newell  et  al.  [30], 

FF 

Three  of  the  four  matrices  shown  in  Table  6.4  are  very 

similar,  the  only  different  one  is  that  designed  as  3RED10(LS). 

FF 

However,  of  these  three  similar  matrices,  the  corresponding 
matrices  are  not  similar,  with  some  of  the  larger  elements  differ** 
ing  by  as  much  as  a  factor  of  three. 


6.4.2  Proportional  Feedback  Plus  Setpoint  Control 

The  proportional  feedback  plus  setpoint  (FB  +  SP)  control 
laws  were  applied  to  the  evaporator  for  two  setpoint  changes  in  C2. 
The  first  changed  C2  from  a  normalized  perturbated  value  of  0  to 
+0.1,  while  the  second  returned  C2  to  0 .  The  form  of  the  control 
law  being  used  is 


u 


„FB  SP  SP 

Sr  -i"+  Sr  z 


(6.51) 


SP 

Thus,  when  _y_  =  0^,  as  it  was  for  the  second  step  (0.1  to  0),  the 

run  which  resulted  was  the  same  as  a  proportional  feedback  run  with 
the  system  initially  at  a  non-zero  state. 


The  response  of  the  evaporator  when  controlled  by  the 


controllers  designed  as 
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a)  3RED10(M)  is  shown  in  Figure  (6.13) 

b)  3RED10(LS)  is  shown  in  Figure  (6.14) 

c)  30PT(M)  is  shown  in  Figure  (6.15) 

d)  30PT(LS)  is  shown  in  Figure  (6.16) 

All  four  control  systems  obtained  the  desired  +10%  change  in 
C2 .  However,  the  controller  designed  as  3RED10(M)  responded  with 
about  half  the  overshoot  (10%)  compared  with  16%  -  19%  for  the  response 
produced  by  the  other  three  control  laws.  Also,  the  response  of  the 
process  variables  is  less  oscillatory  when  controlled  by  the  control 
law  designed  as  3RED10 (M) . 

There  is  a  larger  difference  between  the  control  produced 
by  the  four  results  in  the  second  part  of  the  runs,  when  K^B  is  con¬ 
trolling  the  process  with  a  non-zero  initial  condition  in  C2 .  All 
four  controllers  eventually  return  the  process  to  its  expected  steady 
state.  Once  again,  however,  the  controller  designed  as  3RED10(M), 
(Figure  6.13),  controls  the  system  in  the  smoothest  manner  without 
excessive  oscillation  in  any  of  the  process  variables.  The  control 
calculated  as  3RED10(LS)  produced  a  most  undesirable  response  with 

very  large  oscillations.  This  response  is  not  surprising.  The 
FB 

which  is  being  used  alone  to  provide  proportional  feedback  control, 
was  designed  along  with  a  setpoint  control  matrix  to  provide  for 
changes  in  the  setpoint  vector.  A  better  response  would  result  if 
had  been  designed  using  least  squares  and  then  ^  was  calculated 
to  provide  for  setpoint  changes. 

The  responses  for  the  non-zero  initial  C2  value  of  both  con¬ 
trollers  designed  using  the  third  order  models  are  very  similar,  as 
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FIGURE  6.13  EXPERIMENTAL  PROPORTIONAL  FEEDBACK,  SETPOINT  CONTROL  (3RED10(M)) 

(PP/SP,  10%C2/FB  +  SP/ 3RED10 (M) / SP22) 
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shown  in  Figures  6.15  and  6.16.  Both  allow  the  levels  to  depart  by 
a  large  amount  from  their  initial  (and  desired)  values.  Both 
responses  are  more  oscillatoty  than  that  produced  by  the  control 
law  designed  as  3RED10 (M) . 


The  control  matrices  in  these  FB  +  SP  runs  are  those  shown 
SP 

in  Table  6.4.  The  two  matrices,  that  were  optimal  for  the  two 


.SP 


reduced  order  models,  are  very  similar.  The  other  two  matrices 
are  considerably  different,  with  some  large  elements  differing  by 


SP 


a  factor  of  three.  Each  of  these  matrices  are  used  with  the  same 
FB  FF 

as  was  used  with  the  except  for  the  controller  designed  as 
3RED10 (LS) . 


6.4.3  Proportional  Plus  Integral  Feedback  Control 

Only  three  reduced  order  proportional  plus  integral  (FB  +  I) 
control  laws  were  applied  to  the  control  of  the  pilot  plant,  and  of 
these,  only  one  resulted  in  satisfactory  control.  The  FB  +  I  control 
matrices  that  were  calculated  are  shown  in  Table  6.5. 


The  only  reduced  order  FB  +  I  control  law  which  successfully 
controlled  the  evaporator  was  designed  as  3RED10 (M) .  This  run  is  shown 
in  Figure  6.17.  The  controller  returned  each  of  the  controlled  vari¬ 
ables  Wl,  W2,  and  C2  to  their  desired  steady  state  values. 

The  response  of  the  evaporator  being  controlled  by  the  con¬ 
trol  law  designed  as  30PT(M)  and  30PT(LS)  is  shown  in  Figures  6.18 
and  6.19,  respectively.  For  each  of  these  runs,  only  the  +20%  change 
in  feedflow  was  introduced  since  these  runs  were  tending  to  go  un- 
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FIGURE  6.17  EXPERIMENTAL  PROPORTIONAL  PLUS  INTEGRAL  FEEDBACK  CONTROL  (3RED10(M>) 

(PP/D,  20%F/FB  +  1/ 3RED10  (M) /FB20) 
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FIGURE  6.18  EXPERIMENTAL  PROPORTIONAL  PLUS  INTEGRAL  FEEDBACK  CONTROL  (30PT(M)) 

(PP/D,  +2G%F/FB  +  1/ 30PT (M) /FB29 ) 
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FIGURE  6.19  EXPERIMENTAL  PROPORTIONAL  PLUS  INTEGRAL  FEEDBACK  CONTROL  (30PT(LS)) 

(PP/D,  +20%F/FB  +  I/30PT(LS)/FB31) 
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stable  as  shown  by  the  oscillations  in  all  the  controlled  and  control 
variables.  The  control  matrices  for  these  two  approaches  are  nearly 
identical,  as  shown  in  Table  6.5. 

The  FB  +  I  control  law  designed  as  3RED10(LS)  was  not  applied 
to  the  pilot  plant  since  its  response  was  unstable  in  the  simulation. 

A  comparison  of  these  matrices  with  the  other  matrices  in  Table  6.5 
shows  the  elements  in  the  result  designed  as  3RED10(LS)  are  consider¬ 
ably  larger  than  the  corresponding  elements  in  the  other  matrices. 

These  larger  elements  cause  the  oscillatory  unstable  response.  A 
comparison  of  the  matrices  designed  as  3RED10(LS)  with  the  FB  +  I 

matrices  of  the  tenth  order  control  law  in  Table  A. 5  is  also  inter- 

FB 

are  nearly 
FB 

the  same  as  the  corresponding  columns  of  K  (columns  3,  7  and  8). 

Thus,  the  approach  3RED10(LS)  has  produced  a  result  which  is  essen¬ 
tially  the  same  as  the  tenth  otder  control  law  with  the  undesired 
FB 

columns  of  K  set  to  zero.  Considering  this,  the  unstable  simulated 
response  is  not  surprising. 

For  purpose  of  comparison,  the  control  of  the  evaporator 

by  a  fifth  order  FB  +  I  control  law  which  was  optimal  for  a  fifth 

order  model  is  shown  in  Figure  6.20.  For  this  control  law, 

T 

x^  =  [Wl,  Cl,  HI,  W2,  C2].  In  the  application  of  this  control  law 
the  state  estimation  procedure  used  by  Newell  [28]  provided  an  esti¬ 
mate  of  the  unmeasurable  state  variable,  Cl.  The  response  shown  in 
Figure  6.20  is  slightly  more  oscillatory,  but  with  smaller  maximum 
deviations  in  C2  and  Wl  than  the  responses  shown  in  Figure  6.17 


esting.  is  almost  the  same  as  K'*', 
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(control  law  designed  as  3RED10(M)).  Thus,  the  control  provided  by 
the  third  order  control  law  is  not  much  poorer  than  that  provided  by 
this  fifth  order  controller.  The  fifth  order  FB  +  I  control  matrices 
used  for  this  comparison  are  shown  in  Table  6.6 

6.4.4  Proportional  Feedback  Control 

To  complete  the  comparison  of  the  control  laws  obtained  by 

a  modal  analysis,  the  proportional  feedback  (FB)  controllers  were 

applied  as  calculated  by  3RED10(M)  and  30PT(M).  These  are  shown  in 

Figures  6.21  and  6.22,  respectively.  The  control  matrices  used  here 
FB 

are  the  shown  in  Table  6.4.  The  response  shown  in  Figure  6.21 
results  in  the  expected  offsets  after  the  first  change  in  feed 
flow  rate.  These  offsets  were  eliminated,  as  expected,  after  the 
feed  flow  was  returned  to  its  original  value.  All  the  controlled 
and  control  variables  react  very  smoothly,  without  oscillations. 

The  response  shown  in  Figure  6.22  is  more  oscillatory  than  that  in 
Figure  6.21.  However,  after  the  first  step  in  feed  flow  there  is  a 
smaller  offset  than  that  shown  in  Figure  6.21.  This  is  because  this 
proportional  controller  (30PT(M))  has  generally  higher  gains  and 
so  has  overcorrected  causing  the  oscillatory  response  and  the 
smaller  offset. 

The  feedback  controller  designed  as  3RED10(M)  was  also  used 
to  recover  the  two  oscillatory  FB  +  I  control  runs  shown  in  Figures 
6.18  and  6.19.  These  "recovery"  runs  are  shown  in  Figures  6.23  and 
6.24.  In  each  case,  after  the  oscillatory  FB  +  I  run  was  stopped  and 
this  FB  controller  started,  the  evaporator  responses  settled  out 
immediately.  This  demonstrates  the  use  of  the  3RED10(M)  FB  controller 
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TABLE  6.6 


CONTROL  MATRICES  FOR  FIFTH  ORDER 


PROPORTIONAL  PLUS  INTEGRAL  FEEDBACK  CONTROL 


8.209 

4.542 

4.238 


-1.239 

0.3714 

1.166 


-3.640 

0.5527 

-0.05915 


0.1391 

-1.284 

12.26 


-15.44 

9.070 

14.22 


1.277  0.02747  -1.432 

0.7931  -0.2981  0.8948 

0.6475  1.936  1.306 
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FIGURE  6.21  EXPERIMENTAL  PROPORTIONAL  FEEDBACK  CONTROL  (3RED10(M)) 

(PP/D,  20%F/FB/ 3RED10 (M) /FB19) 
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FIGURE  6.22  EXPERIMENTAL  PROPORTIONAL  FEEDBACK  CONTROL  (30PT<M) ) 

(PP/D,  20%F/FB/30PT(M)/FB33) 
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in  controlling  a  process  at  an  initial  non-zero  state.  The  response 
of  the  evaporator  to  the  -20%  feedflow  change  in  Figures  6.23  and 
6.24  is  similar  to  that  for  the  -20%  feedflow  change  in  Figure  6.21. 

The  response  of  the  evaporator  controlled  with  the  fifth 
order  controller  as  50PT  is  shown  in  Figure  6.25.  It  shows  smaller 
maximum  deviations  than  the  responses  controlled  by  the  third  order 
controller  with  smooth  controller  action  in  Bl,  B2  and  S.  Thus,  the 
third  order  FB  control  law  provides  only  slightly  poorer  control  than 
a  fifth  order  controller.  The  fifth  order  FB  matrix  used  for  this 
comparison  is  shown  in  Table  6.7. 

6.4.5  Summary  of  Experimental  Results 

The  design  approach  which  produced  the  most  consistent, 
stable,  smooth  acting  reduced  order  controllers  of  all  the  types  con¬ 
sidered  was  that  which  obtained  the  third  order  control  law  by  reduc¬ 
ing  the  order  of  the  tenth  order  controller  using  the  modal  analysis> 
(3RED10(M)).  This  is  especially  evident  for  the  case  of  proportional 
plus  integral  feedback  control  where  it  is  the  only  reduced  order 
controller  to  produce  a  stable  experimental  response.  This  design 
approach  also  produced  control  laws  which  resulted  in  nearly  as  good 
control  as  did  the  fifth  order  control  laws  designed  from  a  fifth 
order  model  and  applied  to  the  evaporator  by  Newell  [28], 

Both  the  third  order  models,  which  were  used  to  calculate 
the  control  laws  as  30PT,  were  very  similar.  This  accounts  for  the 
similarity  of  the  resulting  reduced  order  controllers  and  of  the 
resulting  controlled  responses.  The  simulated  comparison  of  pro- 
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FIGURE  6.25  EXPERIMENTAL  FIFTH  ORDER  PROPORTIONAL  FEEDBACK  CONTROL 
(PP/D,  20%F/FB/50PT/FB17) 


TABLE  6.  7 


CONTROL  MATRIX  FOR  FIFTH  ORDER 


PROPORTIONAL  FEEDBACK  CONTROL 


6.508 

-1.229 

-3.235 

-0.08719 

-13.14 

3.825 

0.3688 

0.6885 

-1.378 

9.708 

3.058 

1.056 

0.1784 

9.750 

11.56 

} 
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portional  feedback  control  obtained  from  these  third  order  models,  was 
also  almost  identical. 


The  approach  taken  in  reducing  the  high  order  control  laws 

with  a  least  squares  analysis  was  unsatisfactory  as  is  shown  by  the 

FB 

results  using  this  method.  Also,  a  different  must  be  used  in  each 

of  the  four  types  of  control  considered.  This  is  undesirable.  If 

a  process  were  being  controlled  by  a  feedback  feedforward  scheme  and 

a  setpoint  change  must  be  entered,  the  entire  control  law  must  be 

changed,  the  new  setpoint  entered  and  the  feedback  feedforward  control 

law  put  back  after  the  new  setpoint  has  been  reached.  This  control 

FB 

law  swapping  would  not  be  necessary  if  the  same  could  be  used  for 

all  control  types.  Thus,  an  alternate  design  approach  would  be  to 
FB 

design  using  least  squares  and  to  design  the  other  matrices  sep- 

FF 

arately.  The  matrix  can  be  designed  to  compensate  for  measurable 

SP 

disturbances  as  was  discussed  by  Newell  et  al.  [30],  to  provide 

for  setpoint  control,  and  to  eliminate  offsets  caused  by  unmeas¬ 
urable  disturbances. 


The  similarity  between  the  simulated  and  experimental  closed- 
loop  responses  is  shown  in  Figure  6.26,  which  reproduces  the  experi¬ 
mental  run  of  Figure  6.17  overplotted  by  the  simulated  response  for  a 
single  +20%  step  change  in  feedflow.  The  one  noticeable  difference 
is  the  time  delay  in  the  experimental  results  that  is  not  adequately 
described  by  the  process  model. 
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FIGURE  6.26  COMPARISON  OF  SIMULATED  AND  EXPERIMENTAL  APPLICATION  OF  3RED10 (M) 

(PP,  10L(2018)/D,  20%F/FB  +  I/3RED10 (M) /FB20 ,  3023B) 


6.5  CONCLUSIONS 


A  unifying  treatment  for  the  methods  of  designing  reduced 
order  control  laws  to  yield  incomplete  state  feedback  controllers  is 
presented.  A  modal  analysis,  used  to  eliminate  selected  state  vari¬ 
ables  from  a  high  order  control  law  produced  a  reduced  order  control 
law  which  controlled  the  pilot  plant  evaporator  better  than  the  other 
methods  tried.  This  reduced  order  controller  also  compared  favorably 
with  the  control  using  a  fifth  order  control  law. 

The  topic  of  model  reduction  is  discussed.  Most  of  the 
modal  methods  which  have  appeared  in  the  literature  can  be  shown  to 
be  one  of  two  basic  results  or  one  of  two  revisions  to  one  of  these 


basic  results. 


CHAPTER  SEVEN 


RECOMMENDATIONS  AND  CONCLUSIONS 

7.1  RECOMMENDED  FUTURE  WORK 

As  a  result  of  this  study  of  model  reduction  techniques  and 
different  approaches  to  the  design  of  reduced  order  control  laws, 
several  areas  for  future  work  were  identified. 

7.1.1  GEMS  COPE 

For  the  simulations  in  this  thesis,  GEMSCOPE  was  found  to  be 
a  convenient  package  of  programs.  All  the  control  law  designs  and 
simulated  runs  were  calculated  using  GEMSCOPE.  However,  there  are 
many  design  aids  and  techniques  that  have  been  developed . that  are 
not  included  as  part  of  GEMSCOPE.  Some  of  these  are  discussed  in 
Chapter  Two  and  include: 

1.  A  general  noutine  to  calculate  a  linear  state  space  model  from 
a  nonlinear  model  by  a  numerical  linearization  procedure 

2.  Provision  for  studying  noisy  processes 

3.  Provision  for  designing  multivariable  control  laws  by  methods 
other  than  discrete  dynamic  programming 

To  make  GEMSCOPE  useful  in  a  wider  range  of  applications,  it  must 
be  expanded  to  include  these  simulation  and  design  techniques. 

7.2.2  Least  Squares  Model  Reduction  Using  Random  Data 

The  approach  to  model  reduction  using  random  data  with  least 
squares  resulted  in  acceptable  reduced  order  continuous-time  and  dis¬ 
crete-time  models  of  the  evaporator.  However,  some  difficulty  was 
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experienced  in  the  reduction  of  the  twentieth  order  state  space  model 
of  the  department's  distillation  column.  A  summary  of  the  work  done 
in  this  area  is  described  by  the  author  [14].  A  worthwhile  future  pro¬ 
ject  would  be  to  determine  the  source  of  the  trouble  described  in  this 
reference. 

7.1.3  Least  Squares  Control  Law  Reduction  Using  Random  Data 

The  use  of  least  squares  with  random  data  to  calculate  ah 
reduced  order  control  law  from  a  high  order  control  law  gave  poor 
results.  Suggested  work  in  this  area  is  to  design  the  reduced  order 
feedback  control  matrix  using  least  squares  plus  random  data,  and 
then  design  the  feedforward,  integral  and  setpoint  matrices  to  account 
for  measurable  disturbances,  unmeasurable  disturbances  and  setpoint 
changes.  This  approach  for  the  design  of  feedforward  control,  given 
a  feedback  matrix,  is  considered  by  Newell  et  al  [12]. 

7.1.4  "Implied  Control" 

An  important  area  of  control  theory  is  to  control  a  variable 
which  cannot  be  measured  on-line,  and  which  requires  a  lengthy  labor¬ 
atory  analysis.  One  approach  to  solving  this  problem  could  be  the 
"implied  control"  technique  outlined  in  Chapter  Five.  Perhaps,  with 
a  more  detailed  examination  of  the  systems  controllability,  observa- 

s 

bility  and  structural  characteristics  a  feasible  "implied  control" 
scheme  designed  by  the  method  presented  in  Chapter  Five,  can  be 


devised . 
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7.1.5  Controller  Design  Using  a  Model  with  Input  Derivation  Terms 

One  of  the  forms  of  the  reduced  order  model  calculated  using 
a  modal  analysis  contains  an  input  derivative  in  its  state  equation 
as 

x=Ax+Bu+  Eu  (7.1) 

Marshall  and  Nicholson  [10]  present  an  interesting  approach  to  obtain 
a  simple  control  law  from  this  model.  However,  since  this  control  law 
is  calculated  from  a  continuous- time  model,  it  must  be  applied  as  a 
continuous  (analog)  controller.  Thus,  it  is  not  useful  in  controlling 
the  process  using  a  digital  computer,  unless  the  sampling  interval  is 
quite  small. 

By  defining  a  new  variable,  x*  Equation  (7.1)  can  be  rewritten 
as 

x  =  A  x  +  B  u  (7.2) 

x  =  x  +  E  u  (7.3) 

and  these  can  be  written  in  the  discrete- time  form  of 

x(I  +  1)  =  i  x(j)  +  a  u(j)  (7.4) 

x(j)  =  x(j)  +  1  u.(j)  (7.5) 

An  interesting  project  would  be  to  extend  the  work  of  Marshall 
and  Nicholson  [10]  to  obtain  a  discrete-time  control  law  from  Equa¬ 
tions  (7.4)  and  (7.5)*  which  would  include  the  effect  of  E  u(j )  term 
in  Equation  (7.5). 
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7.1.6  Output  Feedback 

Most  of  the  approaches  to  output  feedback  control  [l  -  9,  13] 
are  for  continuous- time  models,  whose  resulting  control  laws  are  also 
continuous- time.  An  extension  of  some  of  these  approaches  to  design 
discrete-time  controllers  that  could  be  applied  to  the  control  of  a 
process  using  a  digital  computer  would  be  useful.  A  direct  compari¬ 
son  of  the  approaches  used  in  this  thesis  with  these  other  output 
feedback  approaches  could  then  be  made. 

7.1.7  Choice  of  x-^  in  Reduced  Order  Control  Law 

In  the  derivation  of  the  reduced  order  control  law  it  was 
assumed  that  x-^  could  satisfactorily  stabilize  the  system.  As  was 
pointed  out  in  Chapter  Five,  this  is  a  reasonable  assumption  for  the 
evaporator  since  x-^  had  been  used  in  previous  multi-loop  control 
studies  [  1 1  ] .  A  better  theoretical  base  for  choosing  x^  would  be 
useful . 

7.2  CONCLUSIONS 

The  GEneral  Multipurpose  Simulation  and  COntrol  PackagE 
(GEMSCOPE)  which  was  assembled  for  use  in  simulation  and  control  sys- 
tem  design  studies,  proved  to  be  a  convenient  package  of  programs. 

All  the  simulations  and  control  laws  used  in  this  thesis  were  calculated 
using  GEMSCOPE. 

The  general  problem  of  obtaining  a  reduced  order  discrete- time 
model  from  a  high  order  continuous- time  model  was  considered  in  detail. 
The  required  discrete- time  model  can  be  calculated  by  two  alternate 


. 
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approaches.  These  approaches  were  shown  to  be  equivalent  for  several 
modal  approaches  to  model  reduction,  but  were  not  equivalent  when  a 
least  squares  model  reduction  approach  was  used. 

The  approach  to  model  reduction  using  least  squares  with 
random  data  was  more  satisfactory  than  the  more  conventional  trajectory 
fitting  approach. 

In  general,  selection  of  the  "best"  approach  to  model  reduction 
depends  upon  the  end  use  of  the  resulting  reduced  order  model.  Both 
the  least  squares  with  random  data  approach  and  the  modal  approach 
resulted  in  reduced  order  models  which  satisfactorily  approximated 
the  response  of  the  high  order  model.  Thus,  they  were  satisfactory 
in  open  loop  comparisons.  Furthermore,  proportional  feedback,  pro¬ 
portional  feedback  plus  feedforward  and  proportional  feedback  plus 
setpoint  control  laws  designed  from  these  reduced  order  models  gave 
adequate  control  of  the  pilot  plant  evaporator.  However,  not  all  of 
the  reduced  order  models  were  good  enough  for  use  in  the  design  of 
proportional  plus  integral  feedback  controllers. 

The  reduced  order  control  law  which  controlled  the  pilot 
plant  evaporator  the  best  was  that  designed  by  the  elimination  of 
selected  states  from  the  high  order  control  law.  The  controllers 
calculated  by  this  approach  were  as  good  as,  or  better  than  the 
control  laws  calculated  by  other  approaches  when  applied  to  the  con¬ 
trol  of  the  pilot  plant  evaporator.  Furthermore,  this  control  law 
reduction  approach  was  straightforward  to  use,  requiring  only  a  high 
order  control  law  and  a  modal  analysis  of  the  system. 
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NOMENCLATURE 


Eq uation 
where 
first 

Alphabetic  used 

A  Continuous-time  state  matrix  (n  x  n)  (1.1) 

B  Continuous-time  input  (controls)  matrix  (n  x  m)  (1.1) 

A 

I  =  5)  (n  x  (m  +  q))  (2.1) 

C  State  to  output  matrix  (p  x  n)  (1.2) 

D  Continuous-time  input  (disturbance)  matrix  (n  x  q)  (1.1) 

d_  Disturbance  vector  (q  x  1)  (1.1) 

E^  Input  to  Output  matrix  ( /  x  m)  (3.34) 

E*  Input  plus  delayed  vector  to  output  matrix  (2.2) 

F  State  feedback  matrix  used  in  decoupling  (m  x  n)  (2.6) 

G  Continuous-time  canonical  input  (controls)  matrix 

(n  x  m)  (5.22) 

G  =  (G,  H)  (n  x  (m  +  q))  (3.38) 

H  Continuous-time  canonical  input  (disturbance)  matrix 

(n  x  q)  (5.22) 

I  Identity  matrix  (2.5) 

J  Continuous-time  canonical  state  matrix  (n  x  n)  (3.38) 

J  Quadratic  performance  criterion  (2.11) 

j  Time  counter  such  that  2£(j)  denotes  2^(0  at  t  =  jT  (1.3) 

K  Control  matrix  (2.8) 

k  Number  of  data  vectors  in  least  squares  solution  (4.8) 

£  Dimension  of  the  reduced  order  state  vector 

M  Matrix  to  transform  system  into  its  canonical 

form  (n  x  n)  (3.4) 

m  Dimension  of  the  control  vector 
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Nomenclature  (Continued) 


m 

Integer 

(2.4) 

N 

Matrix  used  in  decoupling 

(2.6) 

N 

Total  number  of  time  intervals  in  the  quadratic 
performance  criterion 

(2.11) 

n 

Dimension  of  the  state  vector 

P 

Dimension  of  the  output  vector 

9 

State  weighting  matrix  used  in  the  quadratic  per¬ 
formance  criterion 

(2.11) 

q 

Dimension  of  the  disturbance  vector 

R 

Control 

weighting  matrix  used  in  the  quadratic 

performance  criterion 

(2.11) 

s 

Final  state  weighting  matrix  used  in  the  quadratic 
performance  criterion 

(2.11) 

s 

Integral 

contribution  to  the  control  law 

(6.14) 

s 

Laplace 

transform  variable 

(2.7) 

T 

Discrete 

time  interval 

(2.4) 

t 

Time 

(5.20) 

u 

Control 

vector  (m  x  1) 

(1.1) 

A 

u 

— 

u 

((m  +  q)  x  1) 

(2.1) 

d 

V 

=  M-1 

(3.8) 

V 

Reference  vector  used  in  decoupling 

(2.6) 

W 

Transfer 

matrix  used  in  decoupling 

(2.7) 

w 

Matrix  used  in  the  least  squares  solution 

(4.6) 

w 

Row  of  W 

used  in  the  least  squares  solution 

(4.8) 

w 

Input 

plus  delayed  vector 

(2.2) 

X 

State 

vector  (n  x  1) 

(1.1) 

y 

Output  vector  (p  x  1) 

(1.2) 

r  .  ■ 
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Nomenclature  (Continued) 


Z 

z 

z 


Matrix  used  in  the  least  squares  solution 
Canonical  state  vector  (n  x  1) 

Augmented  vector  used  in  least  squares  control  law 
reduction 


(4.6) 

(3.4) 

(6.46) 


Greek 

a 

0 

A 

A 

A 

| 

A 

o 

r 

v 

e 

1 

d> 


Discrete-time  canonical  state  matrix  (n  x  n) 

Time  weighting  factor  in  the  quadratic  performance 
criterion 

Discrete-time  input  (controls)  matrix  (n  x  m) 

=  (A,  0)  (n  x  (m  +  q)) 

Discrete-time  canonical  input  (controls)  matrix 
(n  x  m) 

=  (5 ,  2)  (n  x  (m  +  q) ) 

Convergence  criterion  in  least  squares  model  reduc¬ 
tion 

Discrete- time  canonical  input  (disturbance)  matrix 
(n  x  q) 

Discrete-time  input  (disturbance)  matrix  (n  x  q) 
Vector  used  in  Davison’s  Modal  reduction 
Discrete-time  state  matrix  (n  x  n) 


(3.6) 

(2.11) 

(1.3) 

(2.3) 

(5.9) 

(3.6) 

(4.31) 

(5.9) 

(1.3) 
(3.20) 
(1.3) 


Subs  cripts 

1,2, 3,4  Refer  to  partitions  of  a  matrix 

I  Partition  corresponding  to  the  integrating  states 

i,j,k  Elements  of  a  vector  or  matrix 

N  Partition  corresponding  to  the  non-integrating 

states 


(3.2) 

(4.14) 


(4.14) 


R 


Refers  to  reduced  order 


(3.3) 
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Nomenclature  (Continued) 


Superscripts 


FB 

Proportional  feedback 

(2.8) 

FF 

Feedforward 

(2.8) 

I 

Integral  feedback 

(2.8) 

SP 

Setpoint 

(2.8) 

T 

Transpose  of  a  vector  or  matrix 

(2.11) 

u 

Intermediate  matrix 

(5.17) 

Nomenclature  for  Computer  Graphs 

Each  computer  graph  is  identified  by  a  string  of  characters. 
These  are  in  the  general  form  of 

(Model  or  Pilot  Plant/Run  Conditions /Control  Type/ 

Source  of  Control  Law/Run  Number) 


The  codes  used  are: 

Model  or  Pilot  Plant: 

PP  Pilot  Plant 

nL (m)  Linear  model  (L)  of  order  n  calculated  in  run  number  m 

Run  Conditions : 

Disturbance 
Set  point 

Initial  condition  of  variable  xx 
Positive  or  negative  change  (no  sign  means  both) 

Step  size  as  a  percentage  of  steady  state 
Process  variable  disturbed 


D 

SP 

xx(0) 

+  >  “ 

X  % 
XX 


Nomenclature  (Continued) 


Type  of  Control  : 


OL 

Open  loop  (no  control) 

FB 

Proportional  feedback  control 

FB  +  I 

Proportional  plus  integral  feedback  control 

FB  +  FF 

Proportional  feedback  plus  feedforward  control 

FB  +  SP 

Proportional  feedback  plus  setpoint  control 

Source  of  Control  Law: 


nOPT 

th 

Optimal  control  law  for  an  n  order  model 

3RED10 

Third  order  control  law  reduced  from  the  tenth  order 
control  law 

DDC 

Multiloop  control 

Run  Number: 

The  simulated  or  experimental  run  in  which  the  response 


was  calculated. 
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APPENDIX  A 


PILOT  PLANT  EVAPORATOR  AND  ITS  MODELS 

A. 1  THE  EVAPORATOR 

The  pilot  plant  evaporator  used  in  this  work  is  a  double 
effect  unit  with  the  two  effects  operating  in  sequence.  It  was  built 
by  Andre  [1].  The  major  pieces  of  process  equipment  are  shown  in  the 
schematic  diagram  in  Figure  A.l.  The  control  loops  shown  in  Figure  A.l 
represent  the  multiloop  control  scheme  applied  to  the  evaporator  in 
previous  studies  [1,6,10,11,13]. 

The  first  effect  has  natural  circulation  through  its  18  inch 
long,  3/4  inch  O.D.  tubes.  It  is  heated  by  process  steam.  The  second 
effect  is  a  long  tube  vertical  unit  which  was  run  in  its  forced  cir¬ 
culation  mode.  It  has  three,  six  foot  long,  one  inch  O.D.  tubes.  It 
is  operated  at  a  lower  pressure  than  the  first  effect  and  is  heated  by 
the  vapour  produced  in  the  first  effect. 

The  evaporator  is  fully  instrumented  and  can  be  controlled 
by  either  Foxboro  electronic  controllers  or  under  Direct  Digital 
Control  (DDC)  from  an  IBM  1800  Data  Acquisition  and  Control  Com¬ 
puter  operating  under  MPX.  Multiloop  DDC  can  be  applied  directly 
using  the  computer  system  control  package  and  advanced  control  schemes 
by  user  provided  programs  and  a  set  of  system  programs  to  interface 
between  the  user  and  system  control  programs. 
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FIGURE  A. 1  PILOT  PLANT  DOUBLE  EFFECT  EVAPORATOR 


A.  2  THE  EVAPORATOR  MODEL 
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The  evaporator  model  was  first  described  by  Andre  [1].  Newell 
[11]  presents  a  complete  outline  of  a  tenth  order  model  in  its  nonlinear 
form.  The  model  used  in  this  work  is  a  linearized  form  of  Newell's  non¬ 
linear  model,  calculated  using  a  numerical  linearization  procedure. 

The  linearized  model  is  in  the  continuous -time  form  of 

x  =  A  x  +  B  _u  +  D  d_  (A.l) 

£  =  C  x  (A. 2) 

or  in  the  discrete-time  form  of  Equation  (A. 2)  and 

x(j  +1)  =  <£  x(j)  +  A  u(j)  +  0  4(j)  (A. 3) 

The  elements  of  the  vectors  x,  _u,  d_,  are  defined  as  normalized 
perturbation  variables  as 

W1  -  W1 

ss 

X3  W1 

ss 

where  W1  is  the  normal  steady  state  value  of  Wl.  The  vectors  x, 
ss  — * 

u,  d_,  and  are  defined  as  follows 

Normal  Steady 

State  Vector,  x:  State  Value 


TS 

Steam  Temperature  (x^) 

244°F 

TW1 

First 

effect  tube  wall  temperature  (X2) 

227°F 

Wl 

First 

effect  holdup  (x^) 

45.5  lb. 

Cl 

First 

effect  concentration  (x^) 

4.59%  glycol 

HI 

First 

effect  enthalpy  (x^) 

189.2  BTU/lb 

, 
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Normal  Steady 

State  Vector,  x:  (continued)  State  Value 


TW2 

Second  effect  tube  wall  temperature  (x. ) 

6 

181°F 

W2 

Second  effect  holdup  (x.,) 

41.5  lb. 

C2 

Second  effect  concentration  (xQ) 

O 

10.11%  glycol 

H2 

Second  effect  enthalpy  (xQ) 

y 

134.1  BTU/lb 

TW3 

Condenser  tube  wall  temperature  (x^q) 

108°F 

Control  Vector,  u: 

S: 

Steam  flow  (u^) 

2.  lb /min 

Bl 

First  effect  bottoms  flow  (u^) 

3.485  lb /min 

B2 

Second  effect  bottoms  flow  (u^) 

1.581  lb /min 

Disturbance  Vector,  d: 

F 

Feed  flow  (d^) 

5 .  lb  /  min 

CF 

Feed  concentration  (d^) 

3.2%  glycol 

HF 

Feed  enthalpy  (d^) 

156.9  BTU/lb 

Output  Vector, 

X  =  (Wl,  W2,  C2)T 

The  coefficient  matrices  of  the  continuous -time  model  are 
shown  in  Table  A.l  and  those  of  the  discrete-time  model,  with  a  64  second 
time  base,  are  shown  in  Table  A. 2.  Table  A. 3  shows  the  system  modal 
matrix,  M,  along  with  the  corresponding  eigenvalues  of  both  A  and  d>. 

The  general  open-loop  model  contains  two  integrating  states 
(with  eigenvalues  of  0  of  unity) .  A  closed  loop  model  is  also  used 
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which  changes  the  unity  eigenvalues  of  to  lie  between  zero  and  one. 
This  is  calculated,  with  a  discrete-time  form  of  the  evaporator  model 
which  has  an  eight  second  time  base,  and  the  relations 
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Bl  =  0.9203  W1 
B2  =  2.768  W2 


(A. 4) 


to  get  the  model  shown  in  Table  A. 4.  This  model  has  an  eight  second 
time  base. 


The  tenth  order  optimal  control  matrices  for  the  model  in 
Tables  A.l  and  A. 2  have  been  calculated  using  dynamic  programming  as 
is  outlined  by  Newell  [11]  and  using  the  following  weighting  matrices. 
For  proportional  feedback,  feedforward  and  setpoint  control 

Q  =  diag  [0,  0,  10,  0,  0,  0,  10,  100,  0,  0]  (A. 5) 

R=  diag  [0.05,  0.05,  0.05]  (A. 6) 

and  for  proportional  plus  integral  feedback  control,  R  is  the  same  as 
Equation  (A. 6)  while  Q  is  given  by 

g  =  diag  [0,  0,  10,  0,  0,  0,  10,  100,  0,  0,  0.5,  0.5,  1.0]  (A. 7) 

The  controller  matrices  for  this  tenth  order  model  are  given  in 
Table  A.  5. 


A  fifth  order  model  is  also  used.  It  is  in  the  general  form 
of  Equations  (A.l)  and  (A. 2)  with  a  state  vector  of 

x  =  [Wl,  Cl,  HI,  W2,  C2]T 
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and  coefficient  matrices  as  shown  in  Table  A. 6.  This  model  is  the  basic 
fifth  order  model  used  by  Newell  [11],  However,  the  coefficient  mat¬ 
rices  are  different  since  it  was  linearized  about  a  different  steady 
state. 
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APPENDIX  B 


CALCULATION  OF  A  REDUCED  ORDER  DISCRETE-TIME 
MODEL  FROM  A  HIGH  ORDER  CONTINUOUS-TIME 
MODEL  BY  TWO  APPROACHES 

A  reduced  order  discrete-time  model  can  be  calculated  from 
a  high  order  continuous-time  model  by  two  approaches,  as  is  discussed 
in  Chapters  Three  and  Four  and  shown  in  Figures  3.1  and  4.3.  Using 
the  modal  analysis  the  same  result  is  calculated  while  using  the 
least  squares  analysis  different  results  are  obtained.  These  conclu¬ 
sions  are  shown  in  this  Appendix.  In  this  Appendix,  the  elements  in 
the  control  and  the  disturbance  vectors  will  both  be  included  in 
one  vector,  called  u_. 

B.l  EQUIVALENCE  OF  RESULTS  BY  MARSHALL'S  ANALYSIS 

B.1.1  <t>  Derived  from  Path  I  (<£>  (I)) 

-ft =R 

By  definition,  (I)  =  e=R  T 

By  Equation  (3.42)  A^  =  M^ Jj  M.^1 
Thus,  using  the  theorem  in  Appendix  D.l 

gR(I)  =  M2  e=lT  M1_1  (B.l) 

B.l. 2  Derived  from  Path  II  (0_(II)) 

-ft =R 

By  path  II,  Equation  (3.13) 

<t>  (II)  =  <t>  -  <6  V  V 

=Rv  ;  =1  =2  =4  =3 

Using  ^  and  ^  from  Appendix  D.3 
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£r(II)  =  Mx  e=lT  (V±  -  l2  V4  1  V3) 


Using  Equation  (D.19) 


0R(II)  =  M1  e=lT  M1  1  (B . 2) 


Thus,  0  by  Marshall’s  analysis  is  the  same  using  the  two 
— R 

approaches  of  Figure  3.1. 


B.1.3  Derived  from  Path  I  (A,,  (I)) 

— K 


r T  A  (T-r) 

By  definition,  A- (I)  =  I  e  B -dr 

'R  Jo  "R 

Using  Equations  (3.42)  and  (3.43)  and  the  theorem  in  Appen¬ 
dix  D.l 

51  -  =2  1a'1  i2'X  (13  h  +  1a  hi  dT 

(B .  3) 


B.1.4  A_.  Derived  from  Path  II  (4^(11)) 

=R -R 

By  Path  II,  Equation  (3.14)  gives 

VID  =  i1+t2  V^1  (I  -  g2)_1  J2 

Using  the  relations  for  A^  and  ^  ^rom  Appendices  D.6  and  D.3, 
and  «2  anci  ^  ^rom  Appendix  D.7, 
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Ar(II)  =  (M1S  M2)  |  e  V  B  dr 


S 


T  J(T-r) 
0 


-1 


<f  -  jn 


+  (M1,  M2)  e  ^?2j  V4  1  1 1  -  e  2  |  (V,,  VA)  \  Me  V  B  dr 


'i'i’  =4' 


j 


T  J(T-t) 

0 


=  (M1, 


h  y  r=Bi 


=  3  =4j  1=2 


dr 


JT 


+  e 


%  Si  +  U  S3>  h  =2  +  =4  M4)  \  e 


£ 


T  J  (T-r ) 


V  B  dr 


This  can  be  rearranged  using  relations  in  Appendix  D.4  to  get 

Jl(T-r) 


=R 


J  (T-r) 

=  V  B  +  e  V  B 

J  (T-r)  1  1  J  (T-r) 

r  h  2i  +  h  h 


dr 


'  ilT  -l 
e  =2  24 


J  T\  ‘  f T  J  (T-r) 

e  1  l  e  (V3  B1  +  V4  B2)  dr 


fT  7J  (T-r)  /* T 

=  M,  l  e  V  B  dr  +  M  I 

=1  Jo  =1  =1  =1  Jo 


T  J,  (T-r) 


dr  +  M,  l  e  1  V2  B2  dr  +  L  +  L2 


where 


ilT  -1  /  il7  \  fT 
=  Si  e  hh  {I'*  )  )0 


J2T  \  1  rT  (T-r) 


<23  il  +  =4  I2>  dr 


=2 


-  (*  *  •*  (■  - 


J  T  /  J  T\  ~\fT  J  (T-r) 

and  L„  =  1 1  +  e  II  -  e  I  lie  (V  B1  +  |2)_  dr 


Expressions  for  V  and  in  Equations  (D.16)  and  (D.17)  give 


r  T  J.CT-r) 

AR(II)  =  M  \  e  M1  (I  -  M2  V3)  |1  dr 


(T  J  (T-r) 

+  Si  jo  <-«r  s2  v  h  ir  +  h  +  h 


T  J,  (T-  r) 


■  h  -  Si 


=1  1  =2  (=3  =1  +  =4  i2)  dr  +  =1  +  =2 


C  T  J-,  (T-  r) 
where  L  =  M  I  e  M  B  dr 

Jo 


Also,  using  the  expression  for  M2  in  Equation  (D.18)  gives 
^(11)  =  L3  +  L4  +  L5  +  L±  +  L2 


T  J  (T-r) 


where  L.  = 
=4 


fT  k 

i e 


Mi  1  A2  V4  1  J2  1  (V  B  +  V4  B2)  dr 


and  =  l  JM„  e 


J±  (T-  r) 


ii  h  k1  h1  %  2i  +  k  S2)  dT 


Now  consider 


=  3  + 


J  (T-r)  _ 

1  -1  -1  ,  / 

M  e  tL 

Si  -  ±2  U  k  %  2l  +  k  S2> 

=  A_(I)  by  Equation  (B.3) 
—  K 


'  i 
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Thus  AR(II)  =  AR(I)  +  L5  +  Lx  +  L2 


Performing  the  integration  on  terms  and  using  the 
result  of  Appendix  D.9  and  on  term  using  the  result  of  Appendix 
gives 


k  ■ 


k  - 


h  “ 


J  T 

-Si  e=1  h  U1  k'1  Si  +  24  12> 


-M2  II  +  e 


(i  +  e^T  (i  -  j  ^  ^  k'1  %  8 

■Sl  -  e=lT\  =2  U1  k'1  <23  8l  +  24  12> 


1  +  24  S2> 


Therefore , 


J,  T 


J  T  -1  =1“  -1 

,  -1  =1  V  V  -M  e  V, 

Ar(II)  =  Ar(I)  +  -K±  V2  V4  1  +  K±  e  =2  =4  =1  =2  =4 


(i  -  ^  -  S; 


j  t  v  j  t  L 

-  S2  u  -  e=Z  ]  -  s2  6=2  i2_1  <23  2i  +  24  s2) 


^(II)  =  4r(i)  -  %  v2  v4  1  +  m2)  J2  1  (V3  B1  +  v4  b2) 


and  using  the  relation  for  in  Equation  (D.17) 


^(H)  -  ^(1) 

/ 

Thus,  the  two  approaches  using  Marshall’s  analysis  give 


.10 


identical  reduced  order  discrete-time  models. 
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By 2  EQUIVALENCE  OF  RESULTS  BY  DAVISON'S  ANALYSIS 

B.2.1  p  Derived  from  Path  I  (d>  (I)) 

— K.  =K 

Davison’s  result  for  A^  using  Equation  (3.44)  is 

Ad  =  A..  +  A„  M„  M  _1 
—  K  — _L  —  Z  —  j  —X 

The  expressions  for  A^  and  A^  in  Appendix  D.2  give 

Ar  ■  Si  [2i  +  =2  S3  «1_1]  +  Sz  i2  [=3  +  U  S3  Si'1] 

Equations  (D.5)  and  (D. 7)  allow  A  to  be  written  as 

— 

ts.  =  Si  ii  Si"1  (b.4) 

Then,  by  the  definition  of  4>  and  the  theorem  of  Appendix  D.l, 

ArT  J,T  _1 

<f>R(I)  =  e  =  M  e  Mx  (B.5) 


B.2.2  0  Derived  from  Path  II  (0  (II)) 

— K  —w 


=R 


</>  from  path  II,  (Equation  (3.22))  is 

— K 


iR(II)  =  +  £  2  =3  =1 


-1 


Using  the  expressions  for  0^  and  0>2  in  Appendix  D.3,  this 
can  be  written  as 

J  T  J  T 

0R(II)  =  Mx  e=1  (Vx  +  V2  M3  M^1)  +  M2  e=2  (V  +  V4  M3  M^”^) 

Making  use  of  Equations  D.5  and  D.7,  this  is  written  as 


J  T 

0R(II)  =  M3  e=1  M1"1 


' i  dJ  fX  baa  Z.C  f.noiJtu  i  lo  *au  v  Le  [ 
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which  is  identical  to  Equation  (B.5) 


B.2.3  Derived  from  Path  (I)  (A  (I)) 


A-CD  by  definition  is 
—  K 


f  T  A_n  (T-  T) 

^  4  e 


1r  dT 


Using  Equations  (B.4),  (3.45)  and  the  theorem  of  Appendix  D.l 


'T  J  (T-t) 

Ad(I)  =  M,  \  e  (V1  B1  +  V2  B2)  dr 


-R 


-Si  j 


0 


B.2.4  AR  Derived  from  Path  (II)  (AR(II)) 
Equation  (3.25)  gives 


4r(II)  =  Mx  (Vx  Ax  +  V2  A2) 


Using  the  definition  of  A  in  Appendix  (D.6) 


$ 


T  J  (T-r) 


-1 


Ar(II)  =  (V15  V2)  ^  Me  M  *  B  dr 


T  J(T-r) 


-1 


=  M1  (Vx  M1  +  l2  M3,  Vx  M2  +  V2  M4)  ^  e  M  ~  B  dr 


Equations  (D.5)  and  (D.6)  give 


4r(ii)  = 


Mi  Q>  2)  ^ 


T  J(T-r) 

e  M  B  dr 

0 


so  that 


Sr 


(II)  =  M1  \  e  (Vx  Bx 


+  V2  B2)  dr 


(B.6) 


which  is  identical  to  Equation  (B.6). 


<*a  d  +  ss  xs 
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Thus,  Davison's  analysis  produces  identical  reduced  order 
discrete-time  models  by  the  two  approaches; 

B. 3  EQUIVALENCE  OF  RESULTS  BY  FOSSARD'S  ANALYSIS 

Fossard's  analysis  produces  the  same  A  ,  B  ,  and  &  as 
does  Davison's  analysis.  The  equivalence  of  this  and  by  the 
two  approaches  in  Figure  3.1  was  shown  in  Appendix  B.2.  The  equivalence 
of  the  E  matrix  by  the  two  approaches  must  now  be  shown. 

B.3.1  E  Derived  from  Path  I  (E  (I)) 

— R  -R 

From  Equation  (3.48) 

Ir(i)  ■  -M2  i2_1  s2  (B-7) 

B.3.2  E  Derived  from  Path  II  (E  (II)) 

— R  — R 

From  Equation  (3.35),  using  Equations  (D.14)  and  (D.15)  and 
Appendix  D.6, 

Er(II)  =  M2  (I  -  a2)_1  |2 

<x3’  V  4 

J2T\_1  rT  J(T-r) 

-  e  J  (Y3>  Y4)  Ml  e  V  B  dr 

J  2T\1  rT  J(T-r) 

=  M2  / I  -  e  J  (V3  Mx  +  V4  M3,  V3  M2  +  V4  M4)  I  e  G  dr 

Making  use  of  Equations  (D.7)  and  (D.8) 


;■ 


iX*iS-eX>.  ( 
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Ir(ii)  - 


M2  (l 

Mifl 


-  e 


-  e 


J0T  ' 

-1 

no 

II M 

Vw' 

i 

=2T  ' 

f  r 

/  JO 

'T  J(T-r) 

e  G  dr 

>0 


=2  dr 


The  integral  in  Appendix  D.9  is  used  to  obtain 

-1 


!r(ii) 


=2  =2 


Thus , 


1r(ii)  ■  -Ml  ii _1  S2 


which  is  identical  to  Equation  (B.7)  above.  Thus,  Fossard's  analysis 
produces  identical  reduced  order  discrete-time  models  by  the  two 
approaches  of  Figure  3.1. 

B . 4  DIFFERENCE  OF  RESULTS  BY  LEAST  SQUARES 

If  the  same  data,  (trajectory  or  random  numbers)  are  used  for 
both  the  continuous-time  reduction  and  the  dis crete- time  reduction  of 
the  two  paths  in  Figure  4.3,  the  resulting  reduced  order  models  will, 
in  general,  be  different.  This  can  be  shown  as  follows 

Equation  (4.25)  shows  that  A^  =  (Z^  Z)  Z^  W 

where  W  is  calculated  from  the  continuous-time  model.  Each  row  of 
=c 

W  is 
=c 

*  T  T  T  T  T 
-1  =  -1  =1  +  —2  =2 

so  that 

T  T 

W  =  Z  A  +  X„  a; 

=c  =  =1  =z  =2 


216 


where  X 2  is  the  matrix  made  up  of  the  data  in  the  3£9(j)  vectors  as 


h  ~ 


4(i) 

4(2) 


4(k> 


Thus , 


T  T  -1  T  T  T 

4  =  (z  z)  z1  (Z  +  X2  A2) 


(ZT  z)  1  ZT  (Z,  x2)( 


=  K 


4 

£ 


where 


Thus , 


K  =  (f  Z)  1  ZT  (Z,  X  ) 


=  ^  =J 


Similarly,  for  the  discrete-time  reduction.  Equation  (4.7)  is 


T  -IT 

tR  =  (Z  Z)  1  z1  wD 


where  Z  contains  the  same  data  as  for  the  continuous  case  above  and 
is  calculated  from  the  discrete-time  model.  Each  row  of  WD  is 


T  T  T  T  T 

x-j^j  +1)  =  Xp  (i )  +  —2  )  ^2 


Md  -  2  4  +  tl 


so  that 


nS  TS  '<5  *  it 
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Now  we  can  arrange  <£R  like  we  did  to  give 

$R  -  t2>  i1 

This  can  be  expanded  using  the  following  relation  of  and 
a  and  relations  in  Appendix  D.3  (a  =  exp  (JT)) 

<t>  =  M  q  V 


($v  i2)  =  (M^  M2)  a  V  =  (M1#M2)  exp  (JT)  V 


Thus , 


£r  =  M2)  exp  (JT)  V  K 


(B.  8) 


The  expression  for  A  can  similarly  be  expanded  into 

=  lt 


^  =  (MrM2)  J  V  KJ 


The  disctete-time  model  corresponding  to  this  A  is 

— R 


£r  =  exp  (ArT)  =  exp  j  (M^  M2)  JTV  K 

This  expression  only  equals  the  Expression  (B.8)  if  V  K1  (M^,  I^)  =  I 
as  is  shown  in  Appendix  D.12. 

Expand  this  expression  using 

K  =  (ZT  Z)'1  ZT  (Z,  X,)  =  (I,  (ZT  Z)"1  ZT  X,)  =  (I,  K,) 


h  =2  \  /  =1 

=  r  %.  «2>  ■  s\g}  (“i-  ife)  =  ^v3  vjy  Ml  k2  m2 


Thus , 


This  is  only  equal  to  the  identity  matrix  if 


^2  =1  =  =3  and  I2  =2  =  =4  (B.9) 

T 

can  be  eliminated  from  Equations  (B.9)  to  give 

M  M  M  =  M 

=3=1  =2  =4 

which  can  be  rearranged,  using  Equation  (D.7)  to  give 

1^2  +  =4=4  =  2  (b.io) 

However,  Equation  (D.8)  says  that 

V  M  4-  V  M  =  I 

=3=2  =4=4  = 

so  that  Equation  (B.10)  and  thus.  Equation  (B.9)  cannot  hold.  Thus, 

T 

VK  (M^ ,  li^)  4  I  and  different  discrete-time  models  are  calculated 
by  the  two  paths  shown  in  Figure  4.3. 


APPENDIX  C 


MODEL  REDUCTION  -  EQUIVALENCE  OF  FORMS 


Many  of  the  modal  reduction  results  presented  in  the  lit¬ 
erature  are  equivalent.  This  Appendix  will  show  these  equivalences. 
In  this  Appendix,  the  control  vector  and  disturbance  vector  will  be 
combined  as  one  input  vector,  _u.  The  high  order  system  will  be 
represented  as 


(SMI; 


while  the  low  order  model  will  be  represented  as 


k  =  h  £1  +  2r  h 

C.l  EQUIVALENCE  OF  NICHOLSON'S  AND  DAVISON'S  RESULTS 


(C.l) 


(C.2) 


(C .  3) 


(C.4) 


Both  Nicholson  [12]  and  Davison  [5]  present  the  same  reduced 

order  model  in  slightly  different  forms.  The  difference  in  their 

results  is  in  the  representation  of  B  . 

— R 


Davison  represents  B  as 

—K 


il  -  h  h 


(C  .5) 
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where  ^  B 


(C.  6) 


Nicholson  represents  B  as 

— R 


Sr  -  Si  2a  S 


(C.  7) 


where  M  is  a  normalized  form  of  M,  ,Q.  is  the  first  /  rows  of  M  ^  as 
=1  =1  -A 


=  <2l’  X2> 


(C .  8) 


_  £  ^ 

and  is  equivalent  to  with  the  i  row  multiplied  by 


m.  . 

li 


m.  . 

li 


(i  =  1,  2 . /) 


The  scalars  m. .  and  m. .  are  elements  of  matrices  M..  and  , 
n  n  =1  =1’ 

respectively. 

D  N 

In  order  to  show  that  B^  =  B^,  Equation  (C.7)  must  be  expanded 

in  terms  of  M..  and  0A . 

=1  -A 

M  is  the  normalized  form  of  M  .  That  is,  each  column  of 

is  obtained  from  by  some  normalizing  multiplication  factor.  Let 

th 

this  factor  be  N_.  for  the  j  column  so  that 


m. .  =  N.  m. . 
ij  3  ij 


(C.9) 


Thus,  there  is  an  N^.  for  each  column  of  M^.  Letting 


N  =  diag  (N1,  N2  ...  N^) 


(C.10) 


matrix  M  can  be  written  as 


Mx  =  Mi  N 


(C.ll) 


5  8 1 


' 


th 


Similarly,  Q^  is  obtained  from  Q^  by  multiplying  the  i  row  by 


m../m. ..  Defining 
11  11  e 


S  =  diagl 


mll  m22 


mll  m22 


m 


'll 


m 


'll 


(C.12) 


Matrix  Q can  be  written  as 


2a  -  S  2a 


(C. 13) 


Now,  Equation  (C.7)  can  be  written  as 


,N 


Br  =  N)  (S  gA)  B 


=  (N  S)  (Vr  V2)  B 

=  M1  (N  S)  G  (C.  14) 


m. 


m„ 


‘11  “‘22 

Now,  N  S  =  diag  [  $  N2  3^,  .  . .  N^  -±JL 


m 


11 


m 


22 


m 


'll 


=  I,  using  Equation  (C.9) 

N  D 

Thus,  B  =  M  G  =  B  as  desired. 

— K  — J_  — 1  — K 

C . 2  EQUIVALENCE  OF  RESULTS  OF  DAVISON,  CHIDAMBARA,  FOSSARD  AND  GRAHAM 

As  mentioned  in  Chapter  Six,  the  same  reduced  order  model 
is  presented  in  the  literature  in  four  places  in  slightly  revised 
forms.  The  result  was  first  presented  by  Chidambara  [3]  as  his  method 
Cl  in  April,  1967,  and  also  by  Davison  [3]  as  a  revision  to  his  orig¬ 


inal  method  in  the  same  correspondence.  Davison  presented  the  result 
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again  in  [4].  The  equivalence  of  these  two  was  stated  by  Chidambara 
[2]  in  a  further  correspondence  in  December,  1967.  Graham  [9]  pre¬ 
sented  a  third  form  in  July,  1968  by  a  different  analysis.  Although 
Graham  referred  to  the  correspondence  by  Chidambara  and  Davison  [3,4] 
he  didn't  recognize  that  his  result  was  identical  to  theirs.  Finally, 
Fossard  [7]  presented  his  form  in  April,  1970  as  an  alternative  result 
to  Davison's  for  use  when  A  and/or  A^  have  zero  eigenvalues.  Although 
some  of  the  equivalences  have  been  stated  before,  they  all  will  be 
shown  here  due  to  the  simplicity  of  the  proofs.  Some  basic  relations 
will  be  used  in  this  work.  These  are  found  in  Appendix  D. 

C.2.1  Chidambara 's  Method  (Cl) 

Chidambara  presented  his  method,  referred  to  as  Cl,  as 


From  Equations  (C.16)  and  (C.17), 


4  =  ii 

+  2i 

u 

(C.15) 

-l  =  4  -i 

+  «2 

—2 

(C.16) 

Z  2  =  F  u  = 

-  k 

CM 

Oil 

rH 

1 

(C.17) 

(C.17), 

=  M  “1  x  _  - 

M  _1 

M_  z„ 

=1 

-1 

=1 

=2  —2 

-1 

-1 

-1 

(C.18) 

-Si 

+ 

Si 

M2  J2  G2  u 

Thus,  from  Equation  (C.16),  using  Equations  (C.15),  (C.17),  and  (C.18), 


we  can  write, 


. 
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k  =  Si  ix  +  m2  k 


=  M1  [jx  +  Gx  u J  +  M2  (-  J2  1  G2  u) 


Jx  (Mx  1  x-l  +  Mx  1  M2  J2  1  G2  u)  +  M1  G1  u  -  M2  J2  1  g2  u 

=i  =1  =1-1  -1  +  =1  (=1  +  il  =1_1  =2  =2_1  =2^  -  "  =2  =2_1  =2  A 


(C. 19) 


C.2.2  Davison's  Modified  Method  (DM) 

Davison  presented  his  revision  to  his  earlier  method  as 

k^=  (A^  +  A2  M3  1)  u 

{  ^  SR  -  [a-1  b]  J  « 

where  J^A  ^  bJ  are  the  first  /  rows  of  A  ^  B 


-i =  H  +< 


(C. 20) 
(C. 21) 


First,  using  Equations  (D.l)  and  (D.2),  A^  can  be  written  as 

A  =  A.  +  A_  M„  M  _1 
— K  — 1  —A  —  j  —1 

=  (M1  J1  V1  +  M2  J2  V3)  +  (Mx  J1  V2  +  M2  J2  V4)  M3  M1 
=  “i  ii  +  v2  m3  m^1)  +  m2  j2  (V3  +  v4  m3  M^1) 

Now,  using  Equations  (D.7)  and  (D.19) 


Sr  -  k  h  -  =2  h1 S3>  +  S2  %  +  h  ar1) 


-1 


-  2?l  ii  ^i 


-l 


(C . 22) 


Now,  from  Equations  (C.20),  (C.21),  and  (C.22) 


- 
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x,  = 


H  + 1  -r’1  Sr  '  [s'1  5^  ]  A 

*i  +  Si  Si  -  +  |=p_1  h  -  [s'1  s]  ii 
Si  -  |=r_1  Sr  -  IV1  ay  »  +  Sx  Si  “  +  j^"1  s 
-  [s'1  S^Ji 

iSi+Si  Si  -  ii  Si'1  4  1  Sr  -  [s'1  l]  J 


■  Si  h  Si 


"  Si  h  Si 


*  Si  Si  Si 


+  {hl  h  -  [s'1  l\L  }  A 


(C. 23) 


Now,  using  Appendix  D.ll,  consider  the  term, 


Sr  *  [S'1  l\ 

=  (si  Si'1  Si'1)  %  Si)  -  h  h1  Si  -  «2  s2_1  g. 


=1  =1  1  Si  "  Ml  Si  1  ^1  "  -2  -2  1  -2 


-  h  =2_1  S2 


(C.24) 


Thus,  the  reduced  order  model  can  be  written  from  Equation  (C.23)  as 


k  =  Si  Si  Si"1  «1  +  Si  (gi  +  Jj  Sx-1  S2  J2_1  g2)  »  -  m2  j 


2'1  §2  A 


(C. 25) 


Equation  (C.25)  is  equivalent  to  Equation  (C.19)  and  so 
Chidambara's  Cl  method  and  Davison’s  DM  method  are  equi valent. 

C.2.3  Fossard's  Revision 

Fossard  presented  the  revision  to  Davison's  method  as 


' 


->%  r>  :  'J 
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—£  =  ^-1  +  =2  =3  =1  ^  — H  +  -1  -i  — 


=  x^  -  ^2  I2  ^  -2  — 


(C.26) 

(C. 27) 


In  Section  C.2.2,  Equation  (C.24)  showed  that 


-  m2  j2  1  g2 


(C. 28) 


so  that  the  correction  applied  by  Fossard  is  identical  to  that  applied 
by  Davison,  and  so  this  result  is  equivalent  to  the  previous  two 
results,  since  Equations  (C.26)  and  (C.27)  can  be  written  as 


+  A2  M3 


+  Si  §1  H. 


which  is  identical  to  Equations  (C.20)  and  (C.21). 


(C.29) 

(C. 30) 


C.2.4  Graham’s  Result 

Graham  presented  his  result  as 

x  =  Mx  J  M1"1  x±  +  (Bx  -  A2  V4_1  J2_1  G2)  u  -  M2  J^1  g2  U 

(C. 31) 


This  is  already  in  the  general  form  of  Equations  (C.19)  and  (C.25) 
except  for  the  representation  of  B^.  from  Equation  (C.19)  can  be 

arranged  into  the  form  as  presented  above.  Thus,  from  Equation  (C.19) 
using  Equation  (D.10), 
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h  ■  Si  %  +  k  sr1  S2  S2_1  S2> 

=  Mi  Si  +  «!  ii  (-  v2  v^1)  j2-x  g2 
■  Ml  %.  Si  +  S2  12)  -  Mx  Jx  V2  v^1  3-1  G 


Now,  using  Equations  (D.2)  and  (D.9) 


2r=  (s  -  M2  v3)  bx  +  m2  v2  |2  -  (a2  -  m2  j2  ff4)  v4_1  j.,'1  g2 


Si  -  M2  S3  Si  +  Mi  l2  S2  -  t2  V*  k'1  £2  +  M2  S2 


-1  ,  -1 


Si  -  k  k  k  S2  -  M2  h  Si  +  Ml  h  I2  +  S2  h  li  +  «2  v4 


Si  -  =2  S4’1  S2"1  §2  +  (Mi  l2  +  M2  V  S2 


Finally,  using  Equation  (D.10) 


=R  =1  =2  =4  =2  =2 


which  is  the  form  of  BR  above  in  Equation  (C.31) 


Thus,  Graham's  result  can  also  be  written  as 


-1 


k  ’  Ml  k  Ml’X  2l  +  +  J2  M^1  M2  J2_1  G2)  u  -  M2  J=-L  G2  u 


-1 


(C . 32) 


which  is  equivalent  to  the  other  three  results  in  Appendix  C.2, 


C. 3  EQUIVALENCE  OF  RESULTS  OF  MARSHALL  AND  CHIDAMBARA 


The  results  presented  by  Marshall  and  by  Chidambara  as  method 


C2  are  identical.  Marshall  presents  A  as 

— R 


while  Chidambara  presents  it  as 


Ar  -  +  A2M3M1-1) 

These  two  forms  of  AR  are  shown  to  be  equivalent  in  Appendix  C.2.2. 

The  form  of  B  presented  by  Marshall  is 

Br  =  (B1  -  A2V4_1  J2_1  (V3|1  +  V4B2))  (C. 33) 

while  that  presented  by  Chidambara  is 

Br  =  Bi  +  (A2M4  -  A2M3M1"1  M  )  F  (C. 34) 

where  F  =  -J2  1  G2  (C.35) 

Equations  (C.34)  and  (C.35)  can  be  combined,  and  rearranged  as  follows 

2r  -  Si  -  h  %  -  m1  sy  ii1  S2 

-  !i  -  a2  %  -  m3m1'1  m2)  I'1  (V3  Sl  +  v4b2) 

Using  Equations  (D.7)  and  (D.8) 

l^1  =  %  ~  =  3=1"" 1  M2) 

so  that 

Sr  "  Si  -  =2=4_1  k1  +  V412) 

which  is  identical  to  Equation  (C.33).  Thus,  Marshall’s  result  is 
identical  to  Chidambara’s  C2  result. 


APPENDIX  D 


MATHEMATICAL  EQUIVALENCES 


D . 1  THEOREM 


Gantmacher  [8,  p.  98]  presents  a  theorem  that  is  used  in 
this  thesis.  It  states: 

If  A  and  B  are  similar  and  W  transforms  A  into  B, 

B  =  W*1  A  W, 

then  the  matrices  f(A)  and  f(B)  are  also  similar  and  W  trans¬ 
forms  f (A)  into  f  (B) , 

f (B)  =  W_1  f (A)  W 


For  the  continuous- time  system  of  Equation  (3.36),  from 
Equation  (3.39), 


From  the  above  theorem,  and  the  definition  of  d>  as  <£  =  exp  (A  T),  it 
follows  that, 

exp  (JT)  =  M"1  exp  (AT)  M  =  M'1  j  M  =  a 


Thus,  the  same  transformation  matrix,  M,  will  transform  both  the  dis¬ 
crete-time  and  the  continuous- time  systems  into  their  corresponding 
Jordan  canonical  forms,  and  M  can  be  calculated  using  either  the  dis¬ 
crete-time  system  or  the  continuous-time  system. 
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D .  2  PARTITIONS  OF  A 


Starting  with  A  =  M  J  V,  it  follows  that 


A  A  rM_  J.  V,  +  J0  V  M.  J,  Vn  +  V, 

=1  =2  =1  =1  =1  =2  =2  =3  =1  =1  =2  =2  =2  =4 


=3  U/  \hhh  +  hhh  hi xh+hiiU. 


Consequently, 


(D.l) 


(D.2) 


D.  3  PARTITIONS  OF  <*> 


Following  from  A  =  M  J  V,  using  the  theorem  in  section  D.l, 


JT 


<f>  =  M  e=  V.  Thus, 


d>  d> 

=  1  =2 


d>  d> 

=3  =4 


£iT 


=1 

S3 


£iT 


il  =2 

h  +  h 


£2t 

h 

£iT 

M  e 

h 

+  h 

£2T 

e 

>11 

£2t 

h 

iiT 

h  e 

h 

+  M, 
=4 

J=2T 

e 

>11 

Consequently, 


J,  T  J2T  .  .  JT  /  V 

ii  =  Si  e=  Si  +  fc  e=  =3  =  <Si*  e=  (  vj 


(D.3) 


J-i T  i2T  ^T  /  V2 

=2  =  =1  6  l2  +  S2  e  V4  =  (Mp  Ife)  e  ( 


(D.4) 


D.4  EXPANSION  OF  VM  =  I 


-1  . 


By  definition  of  V  =  M  ,  it  follows  that  VM  =  I  and 
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(h  «1  +  h  %3  Zitk  +  h  Ma\  /l  2  \ 

\h  Si  +  h  h  h  +  =4  S  J  \2  i  / 

which  implies 


h  ft  +  h  ft  =  i 

=1  =2  +  =2  =4  =  2 

ft  ft  +  ft  ft  =  2 
ft  ft  +  ft  ft  =  i 

D . 5  EXPANSION  OF  MV  =  I 

In  a  manner  similar  to  Appendix  D.4,  MV  =  I  and 
/«!  h  +  «2  h  +  h  =4 

\=3  Zi  +  £4  S3  J2  +  »4  =4 

and 

M  V  +  M  V  =  I 
=1  =1  =2  =3  = 

Si  h  ±  h  U  =  2 

M-  V,  +  M.  V  -  0 
=3  =1  =4  =3  = 

ft  ft  +  ft  ft  =  i 


(D.5) 
(D.  6) 
(D.7) 
(D.8) 


(D.9) 
(D. 10) 
(D. 11) 
(D. 12) 


D . 6  PARTITIONS  OF  A 


D. 


1, 


Starting  with  A  = 
it  follows  that 


A(T-t) 

e= 


B  dr  and  the  theorem  in  Appendix 


A* 


e  0 

J.(T-r) 
0  e  1 


and  therefore 


Si  ■  $0 


J  (T-r) 
e  0 


0  e 


j2(t-  r) 


V  V 
=  1  =2 


V  V 
=3  =4_ 


12 


dr  (D. 13) 


D . 7  DEFINITIONS  OF  c*2  AND 


JT  =2^ 

Since  a  =  e=  ,  then  «2  =  e 


(D.14) 


Since  £  =  V  A,  then  6^  =  (V^,  V^)  A 


(D- 15) 


D.  8  RELATIONS  FOR  Vp  V  ,  M2,  M. 


-1 


-1 


From  Equation  (D.9),  (I  -  V3) 


(D. 16) 


-1 


From  Equation  (D.10).  V_  =  -  M.  ML  V, 

=2  =1  =2  =4 


(D- 17) 


Rearranging  Equation  (D.2),  can  be  calculated  as 


=  A0  V,'1  J  “1  -  m  J,  V.  V,  1  J0  1 
=2  =2  =4  =2  =1  =1  =2  =4  =2 


(D. 18) 


Eliminating  M2  from  Equations  (D.9)  and  (D.10)  gives 


' 


■ 


bsvI;  (0  3)  ut >z  <  U)  c  &  m  , 
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M^1  “  (VL  -  V2  V4_1  V3)  (D.19) 

D . 9  INTEGRAL  EVALUATION  NO.  1 


The  following  integral  can  be  evaluated  as  shown  here 


D.10  INTEGRAL  EVALUATION  NO.  2 
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Starting  with  the  equality  A  =  M  J  V; 


A_1  B  =  M  J_1  V  B 


X  h' 

"j^1  2 

r—f 

>ii 

h 

r 

1 

Slf 

^1? 

_ 1 

2  i2'\ 

_h 

=4_ 

Now,  j^A  bJ  is  the  first  /  rows  of  4  ^  S  as 


[A'1  s],  ■  fc) 


ii'1  2\/Si 


S  =2  V\g2 


=  M1  J1"1  Gl  +  H2  J  _1  G2 


D. 12  CONDITIONS  FOR  WHICH  EXP  (A  B  C)  =  A  EXP  (B)  C 


Consider  the  power  series  expansion  of  exp  (A  B  C)  as 

(A  B  C)  (A  B  C)  (A  B  C)  (A  B  C)  (A  B  C) 
exp  (ABC)  =  I  +  ABC  +  - yT~ -  +  - - —  +  .  .  . 

if  C  =  A  ■*"  this  can  be  written  as 

2  3 

ABC  ABC 

exp  (A  B  C)  =  A  (I)  C  +  A  B  C  +  z—p -  +  — ~y  ~  +  .  .  . 


PQ  II  CQ  II 
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=  A  exp  (B)  C 


Thus,  exp  (A  B  C)  =  A  exp  (B)  C  if  A  =  C 


APPENDIX  E 


LEAST  SQUARES  APPROACH  TO  CONTROL  LAW  DESIGN 

E.l  MINIMIZATION  WITH  RESPECT  TO  THE  CONTROL  LAW 


The  desired  form  of  the  reduced  order  feedback  control  law 


is 


„FB 

i  =  Ir  *1 


or 


— 

fb 

U-, 

K_, 

1 

•  R1 

n 

.2 

• 

• 

FB 

u 

L 

m 

— R 

M  _ 

_  m  _ 

(E.l) 


Consider  only  the  first  of  these  equations,  that  for  the  scalar  u^  as 

(E.2) 


„FB 

U1  =  \  ^ 


*1 


(*?) 


(E.3) 


Rewrite  Equation  (E.3)  for  the  k  data  points  used  in  the  least  squares 
solution 


or 


—  — 
u^  (1) 

"  t  ~ 
2^(1) 

(2) 

• 

ii 

T 

x£(2) 

• 

• 

u^  (k) 

• 

•  rp 

(k) 

(E.  4) 


(E.  5) 


We  want 
between 


:o  calc 


ft) 


in  order  to  minimize  the  sum  of  the  squares 


:he  left  side  and  the  right  side  of  Equation  (E.5)  as 
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Minimize  S  to  give 


-2 


T 

Z  w 


+  2 


T 

Z  Z 


(E.7) 


Thus , 


=  (ZT  Z)  1  ZT  Wj_ 


(E.8) 


The  steps  from  Equations  (E.2)  to  (E.8)  can  be  done  for 
each  of  the  equations  represented  in  Equation  (E.l)  and  the  results 
will  be 


(E.9) 


etc 


Thus,  Equations  (E.8)  and  (E.9)  can  be  augmented  together  to  give 
one  matrix  equation  for  the  resulting  control  law,  as 


(*)'  ■  (v 


kJb 


k£B  \  =  (|T  Z)  1  ZT  (wlf  w2,  ...  wj 

=  (zT  z)-1  zT  w 


m 


(E. 10) 
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If  the  desired  reduced  order  control  law  is  in  the  form  of 


u 


„FB 

£r 


,  „FB  ,  ,  SP  SP 
+  d  + 


=  R  - 


=  R 


+  =R 


(E. 11) 


it  can  be  rewritten  as 


(E. 12) 


where 


and 


(E.13) 


(E. 14) 


and  the  procedure  is  continued  as  above. 


In  the  evaluation  of  Equation  (E.10),  the  elements  of  W  are 
evaluated  using  the  high  order,  optimal  control  law,  so  that  the  rows 
of  W  are  as 


or 


and  W  becomes 


(E. 15) 


(E. 16) 


(E.17) 


f 


I 
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where 


T 

x'U) 

T 


(E. 18) 


Now,  Equation  (E.10)  can  be  written  as 


(E. 19) 


E. 2  MINIMIZATION  WITH  RESPECT  TO  THE  HIGH  ORDER  CLOSED-LOOP  MATRIX 


In  this  second  formulation,  the  desired  resulting  reduced 
order  control  law  is  still  as  in  Equation  (E.l)  but  the  minimization 
is  done  with  respect  to  the  system’s  closed  loop  matrix.  The  closed 
loop  form  is 


FB 


x(j  +1)  =  £•[.  Xj_(j)  +  ^2  -2^  +  =  =R  — 1^ 


■  (si +  £  SrB)  % 


(j )  +  ti  —2^) 


where 


&-1V  £-21 


—12 *  *  £-22 

•  • 

•  • 

•  • 

^ln’  ^2n 


(E. 20) 


(E. 21) 
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Also,  A  = 


'n 


Consider  only  the  first  equation  of  (E.20),  for  scalar  x^(j  +1)  as, 


(j  +1)  =  C_2  ^(j)  +  ±21  -2^^ 


xJ(J  )  £±  +  x^Cj  )  ^21 


(E. 22) 


where 


%  ■  *ii  * 


£■ ) 


Writing  Equation  (E-22)  for  k  data  points  gives 


x1(l)  =  4(1)  C±  +  x^1)  ^21 


xx(2)  =  x^2)  +  4(2)  4i 


x1(k)  =  4(k)  c1  +  4(k)  4 


\ (E. 23) 


Lump  Equation  (E.23)  to  form  one  matrix  equation  as 


-i  =  2  £i  +  h  *2i 


(E. 24) 


FB 

We  want  to  calculate  K  by  minimizing  the  sum  of  the  squares  of  the 

— R 

difference  between  the  left  side  and  right  side  of  Equation  (E.24), 


s2  -  (w x  -  z  c  1  -  z2  ^21)  (w^  -  Z  Cj  -  z2  021) 
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■  4  <»i  -  2  -  h  4i>  -  ^  ST  i  "  2  5-1  -h  —2i^ 


-21  §2  I  £i  ?2 


with  respect  to  as 


%  '  ‘  jI  -i  '  sT  (-i  '  2  5i  -  l2  +  ST  2  ^  +  f  h  <t£  “  5 


T  -1  T  T 

c_i  =  (z1  z)  1  z  (wx  -  z2  ^21) 


(E.25) 


The  steps  from  Equations  (E.22)  to  (E.25)  can  be  done  for  each  of  the 
equations  in  (E.20)  to  give 


T  -1  T  T 

c2  =  (Z  z)  x  z  e»2  -  z2  ^22) 


T  -1  T  T 

C  =  (Z  Z)  1  Z1  (w3  -  Z2  023) 


T  —1  T  T 

Z  (Wj  -  Z2  *2„> 


-\ 


>(E.26) 


Thus,  Equations  (E.26)  and  (E.25)  can  be  augmented  together  to  form 
the  resulting  suboptimal  closed  loop  matrix  as 


(C1.  c2 . Cn)  =  (ZT  Z)  1  zT 


&L»  —2*  ’  *■*  — n^  =2  ^—21  22  ’  *  * ’— 2n^ 


T  -1  I  T  T 

C  =  (Z Z)  <  z  w  -  z  z2 


sl! 


(E.27) 


*' 
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Also 


}  (E.28) 


J 


Using  Equations  (E.27)  and  (E.29), 


si  *  (s?)1  s 


AT  =  (zT  z)  1  |zT  w  -  zT  z2  4>2( 


FB  T 

which  can  be  solved  for  ,  if  (A  A)  exists,  as 

=  j(=T  2)-1  f  l2  =2}  "  =l|  =  (=T  =)_1 


(E. 30) 


T  -1 

The  matrix  (A  A)  will  usually  exist.  One  case  where  it  will  not 
exist  is  if  A  has  a  column  of  zeros,  which  is  unlikely  since  then 
there  is  an  input  which  does  not  affect  the  system. 


In  the  calculation  of  Equation  (E.30),  the  elements  of  W 
are  evaluated  using  the  high  order  optimal  closed  loop  system,  so 
that  the  rows  of  W  are  as 

x(j  +  1)  =  (J  +  A  KFB)  x(j) 
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or 


W  =  (Z,  Z_) 


‘  ,  ,  A  JB,T 

(^1  +  A  ^  ) 


,  ,  .  *  FB.T 

(^2  +  A  |2  ) 


T  T 

W  =  Z  (^i  +  A  Kf)  +  Z2  (^2  +  A  K2B) 


(E.31) 


Now  Equation  (E.30)  can  be  written  as 


T  -1 
(z1  Z)  x 


ZT  Z  (^  +  A  KBB)  +  ZT  Z2  (^2  +  A  Kf) 


T  T 

-  1  S2  *2 


T  I  T  — 1 

-  £-l  (  A  (A  A) 


£  +  ( K™')  AT  +  <ZT  Z)-1  ZT  Z2  (k2FB)T  AT  -  * 


T  /  .  /aT  aN-1 

A  (A  A) 


T  -1 
+  (Z  Z)  1 


(E. 32) 


which  is  identical  to  Equation  (E.19)  formed  by  the  simpler  approach 
of  a  least  squares  fit  of  Equation  (E.l). 


Figure  E.l  shows  a  time  domain  response  comparing  the  high 
order  evaporator  model  optimally  controlled  with  the  high  order  model 
controlled  by  reduced  order  controllers  calculated  with  each  of  these 
two  approaches.  The  reduced  order  controllers  give  identical  results. 


j  nil  nt  •3amvd3dwn 


NIW/Bl  NI  31VfcWQU  NIW/B1  NI  KJU  MV315 


TIME  IN  MINUTES  ,  TIME  IN  MINUTES 

FIGURE  E.l  SIMULATED  COMPARISON  OF  THIRD  ORDER  CONTROL  LAWS  OBTAINED  BY  LEAST  SQUARES  REDUCTION 
(10 (2018) /D,  +20%F/FB/ 10QPT ,  3RED10(LS),  3RED10 (LS) /3040A,  3034,  3050) 


APPENDIX  F 


REDUCED  ORDER  MODELS 

The  reduced  order  models  used  in  this  thesis  are  shown  in 
this  Appendix.  Table  F.l  is  an  index.  The  models  are  presented  in 
Tables  F.2  to  F.16. 
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TABLE  F.l 

INDEX  TO  REDUCED  ORDER  MODELS 


Table 

No. 

Reduction 

Method* 

Model 

Order 

C  or 

p  ** 

Run 

No.  *** 

Comments- 

F-2 

MM 

5 

C 

2020 

F-3 

MM 

4 

C 

2059A 

Keeping  eigenvalues  1,2 

F-4 

MM 

4 

c 

2060A 

Keeping  eigenvalues  1,2 

F-5 

MM 

3 

c 

2021 

F-6 

MM 

3 

D 

2050 

F-7 

DM 

3 

D 

2065 

F-8 

FM 

3 

D 

2094 

F-9 

RLS1 

3 

D 

2119 

Start  sequence  at  19653 

F-10 

RLS1 

3 

D 

2137 

Start  sequence  at  9639 

F-ll 

RLS1 

3 

D 

2136 

Start  sequence  at  1 

F-12 

RLS1 

3 

C 

2141 

Start  sequence  at  19653 

F-13 

DISC 

3 

D 

2141 

F-14 

RLS1 

3 

D 

2113 

Closed-loop  model 

F-15 

RLS2 

3 

D 

2108 

Closed-loop  model 

F-16 

TFLS 

3 

D 

2144 

Closed-loop  model 

*  MM 
DM 
FM 
RLS1 


=  Marshall's  Modal  Approach 
=  Davison's  Modal  Approach 
=  Fossard's  Modal  Approach  ^ 

=  0  by  random  least  squares,  &  for  steady  state. 
—  R  —  w 


SR 


A 


0  and  A  by  random  least  squares. 
— R  — R 


A 


RLS2 

TFLS  =  0R  by  trajectory  fitting  least  squares,  for  steady  state 
DISC 


=R  ^  . y  - - 43 - n - ’  =R 

Discrete  version  of  model  in  Table  F.12. 


**  C  =  Continuous -time  model,  D  =  discrete-time  model 


***  Simulated  run  in  which  the  model  was  calculated. 


■ 

i *»l  snlsJtt 


FIFTH  ORDER  CONTINUOUS-TIME  (MARSHALL’S)  MODEL  (2020) 
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